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Abstract 

The approximate symmetries of Quantum ChromoDynamics in the infinite heavy quark {Q = c, b) 
mass limit (rriQ oo) and in the chiral limit for the light quarks (m, — » 0, q — u, d, s) can be 
used together to build up an effective chiral lagrangian for heavy and light mesons describing strong 
interactions among effective meson fields as well as their couplings to electromagnetic and weak 
currents, including the relevant symmetry breaking terms. The effective theory includes heavy {Qqj 
mesons of both negative and positive parity, light pseudoscalars, as well as light vector mesons. We 
summarize the estimates for the parameters entering the effective lagrangian and discuss in particular 
some phenomenologically important couplings, such as qb^Btv The hyperfine splitting of heavy 
mesons is discussed in detail. The effective lagrangian allows for the possibility to describe consistently 
weak couplings of heavy [B, D) to light (tt, p, K* , etc.) mesons. The method has however its own 
limitations, due to the requirement that the light meson momenta should be small, and we discuss 
how such limitations can be circumvented through reasonable ansatz on the form factors. Flavour 
conserving (e. g. B* ~> B"f) and flavour changing (e. g. B — > K* 7) radiative decays provide another 
fleld of applications of effective lagrangians; they are discussed together with their phenomenological 
implications. Finally we analyze effective lagrangians describing heavy charmonium- like (QQ) mesons 
and their strong and electromagnetic interactions. The role of approximate heavy quark symmetries 
for this case and the phenomenological tests of these models are also discussed. 
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1 Introduction 



There is a general agreement at the present time that quantum chromodynamics (QCD) is the correct 
theory of strong interactions. Although QCD is simple and elegant in its formulation, the derivation of its 
physical predictions presents however arduous difhculties because of long distance QCD effects that are 
essentially non perturbative. Related to them is, for example, the most prominent expected implication 
of QCD, color confinement. 

Inevitably, QCD effects enter any calculation of processes involving hadrons, such as electroweak 
transitions between hadronic states. Predictions for such transitions and their comparison with data are 
essential to complete the program of determining the parameters of the standard electroweak model. The 
main source of uncertainty for such predictions is our inability to calculate the relevant non perturbative 
QCD effects. 

The theoretical progress in the field has gone through various directions, including lattice simulations 
and the use of sum rules, but one framework has emerged as basic to advance our understanding, namely 
the one employing approximate symmetries, broken explicitly or spontaneously, or both ways. 

The empirical pattern of quark masses, that are widely different, is the essential logical guide to the 
formulation of the symmetries that have been introduced. Historically, different roads were followed, 
some symmetries being already known and investigated even before the notion of quark was established. 
The first important development was isotopic spin, vastly used already in the physics of nuclei, suggested 
by the approximate equality of proton and neutron mass. In the quark language it is the closeness of 
the masses of the up and down quark that induces isotopic spin symmetry. The strange quark being 
much heavier, the extension of the SU{2) isospin symmetry to SU{3), to include the strange quark, then 
necessarily implied dealing with stronger symmetry breaking effects. Later on it was realized that there 
is a typical energy scale of hadronic phenomena, such that it is the relative magnitude of the symmetry 
breaking mass parameters, as compared to such a scale, which suggests the degree of accuracy of the 
symmetry predictions. 

From this point of view the magnitude of the SU (2) breaking was generally expected to be related to 
the ratio of the up and down quark mass difference to the hadronic scale (plus the effects of electromag- 
netism, which breaks isospin as well). Both masses are now known to be very small in comparison to the 
scale, which suggests a larger symmetry, SU{2) x SU{2), the light quark chiral symmetry, exactly valid 
in QCD in the limit when both the up and down quark have zero mass. 

Spontaneous symmetry breaking takes place and breaks the chiral symmetry into isospin, thus ex- 
plaining the better experimental viability of isospin in strong phenomena as compared to chiral symmetry. 
Historically, the progress went the other way around, with chiral symmetry proposed before the quark 
mass values were roughly known. Basic to this progress was the interpretation of the pion as the Goldstone 
boson of the spontaneous symmetry breaking. 

One can attempt to treat the strange quark as a massless quark in some first approximation, ready 
to deal subsequently with substantial deviations from the symmetry. The approximate chiral symmetry 
is then extended to chiral SU{S) x SU{3). 

Within such a frame, current algebra provides for a number of useful results. The other useful 
approach is a systematical lagrangian expansion, known as chiral perturbation theory. In this approach 
the symmetry is used to provide for a catalogue raisonne of the terms appearing in the chiral expansion. 
In this way one can determine which phenomenological inputs are needed to fix at a given chiral order 
the full expansion and decide how to determine them from experiment 

At the opposite side with respect to the hadronic scale of QCD, are the heavier quark masses, i.e. 
those of the beauty b and charm c quark. In the limit of infinite masses (mc, mi, — > oo), three phenomena 
appear. 

The first one consists in the fact that the resulting effective lagrangian exhibits a SU{2) heavy flavour 
symmetry; this symmetry applies to quantities that remain finite in the limit mq —^ oo and arises because, 
in such extreme limit, the exact value of the heavy quark mass plays no role in its interaction with the 
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light sector. For finite quark masses the heavy flavour symmetry is broken, and the breaking can be 
relevant especially in the charm sector, since the c quark is substantially lighter than the b quark. 

The second phenomenon is a heavy quark velocity superselection rule, which is due to the fact that the 
strong interactions of the heavy quark, in the rng — s- oo limit, do not change its velocity vq that always 
remain equal to the heavy meson velocity (only weak and electromagnetic interactions can change vq). As 
a consequence of the velocity superselection rule, the effective lagrangian describing strong interactions of 
the heavy quarks should be written as a sum of terms that are diagonal in the velocity dependent heavy 
quark field operators. 

The last phenomenon appearing in the limit of infinite quark mass is the decoupling of the gluon from 
the quark spin; in other words the effective lagrangian is invariant under heavy quark spin transformations 
and has, therefore, a further SU{2) spin symmetry. In conclusion, the complete symmetry of the effective 
lagrangian is a SU{2Nf) of fiavor {Nf is the number of heavy fiavours) and spin for each value of the 
heavy quark velocity. The resulting effective theory is nowadays known as Heavy Quark Effective Theory 
(HQET) (see |3| ^, |^, ^, ^ ||). In the physical world the symmetry is broken explicitly because of 
the finite heavy quark masses. Symmetry breaking terms are expected to be particularly important for 
charm quark and they can be systematically added to the lagrangian of HQET and parameterize, order 
by order, the deviations from the heavy mass limit. 

One of the first and most important applications of the heavy quark symmetries has been the study 
of the semileptonic decay of i? — > D*W. In the infinite quark mass limit this process is described by one 
form factor whose normalization is fixed at the kincmatical point where the two heavy quarks have the 
same velocity. The velocity of the heavy quark in this limit is, as we have already stressed, the velocity 
of the meson. 

To illustrate in more detail the usefulness of the heavy quark symmetry one can consider the analogy 
between the determination of the Vus, element of the CKM matrix from the semileptonic decay K — > ireD 
and the possible determination of the element Vcb from the semileptonic decay of B into D*. For the 
if to TT decay, a non renormalization theorem says that corrections to the SU{3) normalization of the 
form factor at the symmetry point, i.e. zero momentum transfer, vanish at first order in the difference 
between the strange quark mass and the nonstrange quark mass. For the heavy transition the symmetry 
is the heavy quark symmetry, which is valid for very large quark masses; in this limit some relevant form 
factors are renormalized only at second order in the symmetry breaking parameter (the inverse quark 
mass) at the relevant (zero velocity) symmetry point 

The last example shows the usefulness of the heavy quark symmetry not only to provide us with exact 
relations valid in the terra firma of the exact limit, but also as a platform for studying corrections away 
from the limit. 

In some kineniatical regions, which, at the same time, are not very far from the heavy quark limit and 
from the chiral limit for the light particles, one can try to use simultaneously both the heavy quark and 
the chiral approach in the two distinct sectors and, as we have already discussed, the most economical 
way to do this consists in using phenomenological lagrangians. In other words, chiral SU{3) x SU{3) 
symmetry can be used together with the spin-flavour heavy quark symmetry of HQET and the velocity 
superselection rule to build up an effective lagrangian whose basic flclds are heavy and light meson 
operators. This approach has been proposed in a number of papers jl^, |ll|, |l2|, ^ ^4, 15, 1^, [ij, |l^, |l^ 
and the purpose of this report is to review this method and its applications to the interactions among 
heavy and light mesons. 

The first results we describe are in the field of the strong interactions and concern the properties of 
the effective fields describing the heavy mesons as well as their couplings to the pseudoscalar octet of 
the Goldstone bosons. We also discuss the introduction in the lagrangian of the light vector resonances 
p, K* etc and the inclusion of positive parity heavy meson states. Applications of these ideas to heavy 
baryons containing one or more heavy quarks have been also studied |2^, ^ |2^, ^ |2^, but they will 
not be reviewed here because the experimental situation concerning heavy baryons is still poor and there 
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are therefore too few constraints on the parameters of the resulting effective theory. 

The chiral lagrangian approach has the advantage of allowing for a perturbative theory including 
not only tree level contributions, but also loop calculations. Such calculations, at present, can only give 
an order of magnitude estimate of the effects, because we do not have yet sufficient experimental data 
to fix the arbitrary coefficients in the counterterms of the effective lagrangian. Nevertheless they offer 
a clue to the size of the loop effects and can be extremely useful in reconciling data with theoretical 
expectations based on tree level calculations. For pedagogical purposes we shall present two explicit and 
detailed examples of these calculations; the first one is the evaluation of the loop effects in the hypcrfine 
— Mb mass splitting. The second example we shall show is given by the chiral loop effects to the 
ratio JdsI Idi where and fo^ are the D and Ds meson leptonic decay constants. Other examples 
of computed chiral loop effects are given by the corrections to the strong coupling constants go* Dm to 
semileptonic form factors and to B and D meson radiative and rare decays: they will be also reviewed and, 
whenever possible, the results of the lagrangian approach will be compared to other existing theoretical 
methods. 

The main pitfall of the effective lagrangian approach is the abundance of coupling constants and 
parameters appearing in the lagrangian. Even if one works at the lowest order in the light meson 
derivatives and in the l/niq expansion, one has to fix several couplings from data. A typical example 
is the already mentioned D*Dtt coupling constant, whose experimental determination is still missing. 
In absence of experimental inputs, one may rely on theoretical information coming, for example, from 
QCD sum rules (for a review of this subject see ||2^), or potential models |^ (for a review see 
psj ) or, when available, on the results obtained by Lattice QCD (for a review of the D and B meson 
phenomenology on the lattice see ||29|| ). An alternative is provided by the use of information coming 
not only from strong interactions, but also from weak and electromagnetic interactions among mesons. 
Actually the application of the chiral lagrangian to these processes offers the possibility not only to exploit 
experimental data to constrain the effective lagrangian, but also to relate different processes using the 
symmetries. This is the second main issue to be discussed in the present report. We shall see that two 
methods can be used to perform the task: the first one uses the chiral and heavy flavour symmetries to 
relate different weak and electromagnetic transitions by establishing scaling relations among them. The 
second method makes use of the chiral lagrangian to compute the different amplitudes. In both cases, 
however, some additional hypothesis on the g^-behaviour of the form factors must be made and we shall 
discuss the different scenarios as well as their comparison with the data. 

The third and final topic discussed in this paper is the application of the ideas of HQET to mesons 
made up by two heavy quarks (heavy quarkonium) . The effective quark theory resulting from the toq 
oo limit satisfies, as in the previous case, the velocity superselection rule and the spin symmetry, but not 
the heavy flavour symmetry. As a matter of fact, the non relativistic kinetic energy term of the effective 
QCD lagrangian, which is flavour dependent, cannot be neglected since it acts as an infrared regulator. 
Therefore the chiral effective lagrangian for light mesons and heavy quarkonium-like mesons does not 
possess the SU(2) heavy flavour symmetry; nevertheless, because of the spin and chiral symmetries, it 
allows for a number of relations among different strong and electromagnetic decay amplitudes of heavy 
quarkonia states: they will be discussed and compared with the data whenever they are available. 

In our opinion the chiral lagrangian approach to the interactions of the heavy mesons is a predictive 
method to relate a large amount of processes and decay rates of these states. We hope to convince the 
reader by this work that the chiral lagrangian method for heavy hadrons is a promising way to describe 
this most fascinating physics. 

We conclude this introduction with a brief summary of the subsequent sections. In section ^ we review 
the symmetries of the approach, we construct the effective chiral lagrangian for heavy mesons and we 
discuss the inclusion in the effective lagrangian of the light vector mesons and the positive parity heavy 
meson resonances. In section || we discuss some problems related to the strong interactions effective 
lagrangian: the strong coupling constant gB*B-n and its possible determinations; the one loop calculation 
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of the B* — B hyperfine splitting and the strong decays of positive parity states. In section ^ after a brief 
review of the i? — > Z?, D* semileptonic transition, we discuss the effective weak V—A current and the chiral 
corrections to the ratio foj fo- Semileptonic heavy mesons decays into a final state containing one light 
meson are discussed in section ^ where we also consider the constraints put on the q^-behaviour of the 
form factors by different theoretical approaches and by some weak non leptonic decay rates, most notably 
B J/ipK* . In section ^ we consider radiative heavy meson decays and we discuss the predictions arising 
from the chiral lagrangian approach. Sections |^ and |^ are devoted to heavy mesons containing two heavy 
quarks: we write down an effective lagrangian describing their interactions and we use it to relate different 
decay processes of these states. In particular we also discuss processes characterized by the breaking of 
the symmetries of the effective theory: spin and chiral symmetry. Finally three appendices conclude the 
work: the first contains a list of Feynman rules used to compute the amplitudes; in the second, some 
integrals encountered in the loop calculations are listed; the last appendix contains the formalism for 
higher angular momentum quarkonium states. 

2 Heavy quark and chiral symmetry 
2.1 Heavy Quark Effective Theory 

The HQET describes processes where a heavy quark interacts via soft gluons with the light degrees of 
freedom. The heavy scale in this case is clearly mg, the heavy quark mass, and the other physical scale for 
the processes of interest here is Aqcd- The identification of the heavy degrees of freedom to be removed 
requires some care: we do not want to integrate out completely the heavy quarks, being interested in 
decays of heavy hadrons and therefore in matrix elements with heavy quarks on the external legs. As we 
will see, the so called small component of the heavy quark spinor field, describing fluctuations around 
the mass shell, has to be eliminated. 

We indicate by z;^ the velocity of the hadron containing the heavy quark Q. This is almost on shell 
and its momentum pq can be written, introducing a residual momentum k of the order of Aqcd, as 

PQ = mqv + k . (I) 

We now extract the dominant part rnqv of the heavy quark momentum defining a new field Qv 

Qv{x) = exp{imqvx)Q{x) = hy{x) + Hy{x) . (2) 

The field is the large component field, satisfying the constraint ^hy = h^: if the quark Q is exactly on 
shell, it is the only term present in (H). the small component field, is of the order l/mq and satisfies 
■^Hy = — i?!,: it is integrated out when deriving the HQET effective lagrangian. 

The non local effective lagrangian is derived by integrating out the heavy fields in the QCD generating 
functional, as done in ]30| . At tree level one has simply to solve the equation of motion for Hy and 
substitute the result in the QCD lagrangian. The equation of motion is 



(2toq + ivD)H, 



I 



i0h. 



(3) 



V ~ 



2 



'V 



where, 



(4) 



with Ta the generators of SU{'i)c and cjs = \/47raI the strong coupling constant. We get 



= K{ivD)K + hyip^^^—^ (- ^——p:] i0hy 

2 \2mq + ivD) 



(5) 
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where a sum over velocities is understood. By using the foUowing identity: 
1+^ 1-/ l+yi . 1+/ 

^^7m^— 7;^^— = — - - ^o'm'')^— ' (6) 
we can write 

Ceff = hviivD)K + KiD^ig^'' -v^v" -laf^")- ^——iD^h, ■ (7) 

zniQ + ivD 

The expansion of this lagrangian in l/mg gives an infinite series of local terms. The leading one is 

C = K,{ivD)K (8) 

which, being mass independent, clearly exhibits the heavy-flavour symmetry. Moreover, since there are 
no Dirac matrices in (|^), the heavy quark spin is not affected by the interaction of the quarks with gluons 
and therefore the lagrangian has a SU (2)-spin symmetry. 

These symmetries are lost if we keep the next terms in the 1 /toq expansion 



C^-^K{iDfK + -^K(T^liG''^K + 0{l/ml) , (9) 
2mQ AniQ ^ 



where we have used the equation of motion, ivDhy = ©(I/toq), to get rid of the term hy{vD)'^hy. 
The first term is the kinetic energy arising from the off-shell motion of the heavy quark, the second one 
describes the chromomagnetic interaction of the heavy quark spin with the gluon field. 

The last step in building up the effective lagrangian is the inclusion of QCD radiative corrections. In 
and (^) the Wilson coefficients are taken at the matching scale mg, i.e. the scale at which the heavy 
degrees of freedom are integrated out. The evolution down to a scale /i < toq introduces logarithmic 
corrections. Details and references can be found in § |. We shaU summarize here only some results. 

The inclusion of quantum loop corrections due to hard gluon exchanges modifies the coefficients in 
the lagrangian (||) and (||), giving: 

C^K{ivD)hy + ^K{iDfK + ^Kao,pG''PK . (10) 
2mQ AniQ 

The tree level matching gives ai{mQ) = 02 (mg) — 1: in the leading logarithm approximation one finds 
at the scale /z < mg 



-9/{33-2Nf) 



(11) 



where Nf is the number of active quark flavours in the range between /i and mg. Notice that ai(/i) = 1: 
this is a consequence of the so-called reparametrization invariance pl| |, which relates the term in ai to 
the leading one. Such an invariance arises from the fact that the decomposition of the heavy quark 
momentum is not unique. The transformation 



V 



k k-q, (12) 

where • g = to satisfy the constraint = 1, is another possible decomposition and it has to give 
rise to the same physical observables: only the heavy quark momentum is a well defined quantity. The 
consequences of this invariance have been studied in ref. The main results are as follows. First of 

all the velocity and the derivative iD should appear only in the combination 

v,+^^ (13) 
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where M is the mass of the field under consideration (in this case toq). Second one has to modify the 
fields in the velocity representation, that is 

(l>^{x) = eyip{iMvx)(j){x) . (14) 

The scalar fields do not require any change, but a vector fields P^^ , at the order f/M should appear in 
the combination 

pr = n*"-«''^^ • (15) 

This is because the field P* should satisfy the constraint v^^P*^ = also after reparametrization. The 
field P^, as well as the scalar field, have a very simple transformation law under the reparametrization 
(|l|). They pick up a phase factor 

Pr->exp(*g.x)P^^ . (16) 

Invariant terms under reparametrization are then easily constructed. In particular one finds ai(fi) = 1 
from these constraints. 

We want now to implement the symmetries discussed before in the spectrum of physical states, in 
particular the pseudoscalar D and B meson states and the corresponding vector resonances D* and B* . 
The wave function of a heavy meson has to be independent of flavour and spin of the heavy quark: 
therefore it can be characterized by the total angular momentum se of the light degrees of freedom. To 
each value of corresponds a degenerate doublet of states with spin J = S£ ± 1/2. The mesons P and 
P* form the spin-symmetry doublet corresponding to Si = 1/2. 

The negative parity spin doublet (P, P*) can be represented by a 4 x 4 Dirac-type matrix H, with 
one spinor index for the heavy quark and the other for the light degrees of freedom. Such wave functions 
transform under a Lorentz transformation A as 

H ^ D{A)HD{Ay^ (17) 

where -D(A) is the usual 4x4 representation of the Lorentz group. Under a heavy quark spin transfor- 
mation S belonging to SU{2) one has: 

H ^ SH , (18) 

where S satisfies S] — to preserve the constraint = H. 
A matrix representation of current use is: 

2 

H = joHho ■ (20) 



H = '-^[p;r^Pi,] (19) 



Here v is the heavy meson velocity, v^P*^ = and AI^ = Mp = Mp- (we shall use also the notation 
Mh = M). Moreover jiH = -H^ ^ H, H^^ -^H = H. 

P*^ and P are annihilation operators normalized as follows: 

(0|P|Qg(0-)) = (21) 

{0\P*''\Qq{l-)) = e^V^ ■ (22) 

The general formalism for higher spin states is given in |Q. Here we will consider only the extension 
to P-waves of the system Qq. The heavy quark effective theory predicts two distinct multiplets, one 
containing a O"*" and a 1+ degenerate state, and the other one a l"*" and a 2+ state. In matrix notations, 
analogous to the ones used for the negative parity states, they are described by 

S=l{l+MD'^l^.l5-Do] (23) 
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and 



T^ = -(l+/) 



D2 Iv 



with the fohowing conditions: 



(24) 



(25) 



These two multiplets have S£ = 1/2 and S£ = 3/2 respectively, where s^, the angular momentum of the 
light degrees of freedom, is conserved together with the spin sq in the infinite quark mass limit because 

J = St + SQ. 



2.2 Chiral symmetry 

From the point of view of HQET it is natural to divide quarks into two classes by comparing their 
lagrangian mass with Aqcd- The u and d quarks belong definitely to the light quark class, to„, md <^ 
Aqcd- The situation for the strange quark is not so clear, but it is usually considered to belong to the light 
quark class, though non negligible mass corrections are expected. If we take the limit niu, rud, rUg 0, 
the QCD lagrangian for these three quarks possesses a SU{3)l SU{3)b. <Xi U{l)v symmetry which is 
spontaneously broken down to SU{3)v <8i U{l)v- The lightest pseudoscalar particles of the octect tt, K, 
K, T] are then identified with the Goldstone bosons corresponding to the broken generators. Of course, 
due to the explicit symmetry breaking given by the quark mass term, the mesons acquire a mass. 

As it is well known, the interactions among Goldstone bosons can be described by the chiral perturba- 
tion theory | |3^ , that is a low momentum expansion in momenta and meson masses. Chiral perturbation 
theory describes the Goldstone bosons in terms of a 3 x 3 matrix S(a;) e SU{3) transforming under 
SU{3)l<?) SU{S)b. as 



The meson octect is introduced via the exponential representation 



S — exp 



/ 



/ = 132 MeV 



(26) 



(27) 



where is a 3 x 3 hermitian, traceless matrix: 



M = 



V 



TT 

K- 



(28) 



3^ / 



To the lowest order in the momenta and in the massless quark limit, the most general invariant lagrangian 
is given by 



f2 



(29) 



where the constant /^/8 has been chosen such as to get a canonical kinetic term for the mesonic fields 
appearing inside the matrix M.. 
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Higher order terms in the momentum expansion are suppressed by powers of p/A^, where p is the 
typical momentum scale of the process and is the chiral symmetry breaking scale, which is evaluated 
to be of the order of 1 GeV. As we have already noticed, chiral symmetry is not an exact symmetry of 
QCD, being explicitly broken by the quark mass term 

^ qarhabqb (30) 

a—u,d,s 

where rh is the hght mass matrix: 

(31) 




























The expression ( |30| ) transforms as the representation (3l, 3ij) ® (3i,, 3_r). We can take into account this 
breaking, at the first order in the quark masses, by adding to the chiral lagrangian a term transforming 
exactly in the same way. This contribution can be written in the form 

AoTr (mS + Y.'^rh) . (32) 

The Goldstone bosons receive a contribution to their square mass from this term. This is the reason for 
treating formally the quark masses as second order terms in the momentum expansion. Then, the tree 
diagrams generated by (|2|) and ( ^ reproduce the same results of the soft pion theorems. Corrections to 
the leading terms come from higher derivative or mass terms and from loop diagrams. It is also important 
to stress that chiral perturbation theory is renormalizable at any fixed order in the momentum expansion. 

The interactions of the Goldstone fields with matter fields such as baryons, heavy mesons or light 
vector mesons (p, w), can be described by using the theory of non linear representations as discussed 
in the classical paper by Callan, Coleman, Wess and Zumino (CCWZ) |Q. The key ingredient in this 
theory is the coset field ^(cc), which is defined on the coset space SU{?>)l ® SU{3)r/ SU{3)v- In this 
context, the ^(x) field is simply related to T,{x) by the relation 

nx)^ei^) ■ (33) 

The transformation properties of ^(x) under chiral transformations (that is transformations of SU{S) l ® 
SU{i)R) are 

ax)^9Lmu\x)^U{xnx)g^j^ . (34) 

The matrix U{x) belongs to the SU{Z)v unbroken subgroup and it is defined by the previous equation. 
As a consequence, U{x) is generally a complicated non-linear function of the coset field S,{x) itself, 
and, as such, space-time dependent. The matter fields have definite transformation properties under the 
unbroken SU{i)v group. For instance, a heavy meson made up by a heavy quark Q and a light antiquark 
qa (a = w, d, s), transforms, under a chiral transformation, according to the representation 3 of S'C/(3)y, 
that is [Ha « Qqa) 

Ha -> HbUlix) (35) 

where U is the same matrix appearing in eq. (p4|). In view of the locality properties of the transformation 
U{x), one needs covariant derivatives or gauge fields, in order to be able to construct invariant derivative 
couplings. This is provided by the vector current 

V,. = l{ed^^ + ^d,^^) (36) 



transforming under the chiral transformation of eq. ( |34D as 

V^^UV^U^ + Ud^U^ . (37) 
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It is also possible to introduce an axial current, transforming as the adjoint representation of SU{3) 



V 



A^ = l{^^d^^-id^^^) (38) 

with 

A^^UA^U^ . (39) 
2.3 A chiral lagrangian for heavy mesons 

The effective lagrangian for the strong interactions of heavy mesons with light pseudoscalars must satisfy 
Lorentz and C, P, T invariance. Furthermore, at the leading order in the l/M expansion (M is the heavy 
meson mass), and in the massless quark limit, we shall require flavour and spin symmetry in the heavy 
meson sector, and chiral SU{3)l ^ SU{3)r invariance in the light one. The most general lagrangian is 
then 

C = i< HbV^D^haHa > +ig < Hb-ff.-f^A'i^Ha > + 

+ ^-^d^llabdA (40) 

where = 9^ + and < . . . > means trace over the 4x4 matrices. In ( pO| ) a sum over heavy meson 
velocities is understood. The first term in the lagrangian contains the kinetic term for the heavy mesons 
giving the P and P* propagators, 

(41) 



2v ■ k 

and 

2v ■ k 



(42) 



respectively. The interactions among heavy and light mesons are obtained by expanding the field ^(x) = 
exp(iA^(a;)//) and taking the traces. In the first term there are interactions among the heavy mesons 
and an even number of pions coming from the expansion of the vector current . The interactions with 
an odd number of pions originate from the second term. As an example, the first term in the expansion 
of the axial current gives 

« ^9^M + . . . (43) 

The last term in (^0|) is the non-linear lagrangian discussed in the previous section, describing the light 
meson self-interactions. Corrections to this lagrangian originate from higher terms in the l/M expansion 
and from chiral symmetry breaking. Let us start with the last issue. We proceed as in the previous 
section by considering, at the first order, breaking terms transforming as (3l,3_r) © (3i,,3ij) under the 
chiral group. The most general expression is 

CxB = Xoirhab^ba + mabi^la) + Ai < HaHb{£,'m^ + ^^m^^)ba > 

+ X[ < HaHa{m^ + m^^)bb > 

+ 1X3 < Hb-r^jsA'^^i^rh^ + £}m£}),aHa > + 

+ iX'3<HaHa-if,-i^A''^a[imi + (}m£})dc> ■ (44) 

Here we have neglected terms contributing to processes with more than one pion. Notice that the 
coefficients A3 and A3 should be of order 1/A,^ because they multiply operators of dimension five. In 
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principle there are other dimension five operators (see |p5|| ), which however contribute only to the order 
X 1/M (neglecting again interaction terms with more than one pion field). The Ai and A'^ terms 
give rise to a shift in the heavy meson propagators. For instance, in the case of the strange heavy mesons, 
they produce the shift v ■ k ^ v ■ k — 5 with S = Md^ — Md — Mb^ — Mb- 

Let us now discuss the 1/M corrections (see refs. ^). First, one has to take into account 

the constraints coming from the reparametrization invariance that tie together different orders in the 
expansion. In the present formalism one can define H fields transforming by the simple phase factor 
exp(igx) under the transformation (n2|) 



H' = H+^^D^[-I^^,H] . (45) 
In fact, it is easily seen that, neglecting terms proportional to the form of the free wave equation (con- 



tributing to the next order in the expansion) is equivalent to use the equation (13) for the four- velocity 
and the equation ( |l5|) for the vector field in the definition of H. This substitution modifies the zeroth 
order lagrangian in the way described in js^. We shall not report here the expressions because all the 
extra terms involve at least two derivatives and they contribute only at the order 1/M^. Finally we have 
0{1/M) terms which are invariant under four- velocity reparametrization, and therefore they appear with 
arbitrary coefficients. In this discussion an important role is played by the time-reversal invariance. In 
our case we have to require invariance under the following transformations 

M{x) -> ~M{-xp) 

A''{x)^A^{-xp) (46) 

where xp and vp are the parity reflections of x and v, that and Vp = v^. Also 

T^T-'^i; . (47) 

Taking into account this constraint and neglecting higher derivative terms (which contribute to the order 
1/AP), one finds [|| 

A2 - 

+ ifj<i,i,A^,^H,Ha> . (48) 

By writing A2 — —MA/2 — —M{Mp. — Mp)/2, one sees that the effect of the corresponding operator is 
to shift the P and P* propagators to 

2{vk + lA) ^^^^ 

and 

- 2(..fc-iA) 

respectively. The couplings 171 and g2 renormalize the coupling g appearing in equation ( ^o| ) in different 
way for the P*P*and P*P couplings. More precisely one finds 

g gp'P' ^g + ^{gi+ 52) 

5 ^ .9p*p = .9+ -^(.91 - .92) • (51) 
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2.4 Light vector resonances 

We want now to introduce in the previous effective lagrangian, equation (^o|), the light vector resonances, 
/o, K* , etc. We shall make the hypothesis that they can be treated as light degrees of freedom. Therefore 
they could be introduced as matter fields by using the CCWZ formalism [Q. However we prefer here 
to make use of the hidden gauge symmetry approach , as done in jl^ (see also ^ ) . The two 
methods are completely equivalent, but the second one is easier to deal with. The main idea lies in the 
observation that any non-linear ct- model based on the quotient space G/ H , where G is the symmetry 
group and H the unbroken subgroup, is equivalent to a linear model with enlarged symmetry G x Hiocai, 
where Hiocai is a local symmetry group isomorphic to the unbroken group H. In the linear model the 
fields have values in the group G, rather than in G/H as in the non-linear formulation. The extra degrees 
of freedom can be gauged away by taking advantage of the local invariance related to Hiocai- In the 
unitary gauge one recovers the CCWZ formulation. However, the explicit appearance, in the formalism, 
of a local invariance group gives room for the introduction of gauge fields with values in the Lie algebra 
of if, which will be interpreted as the light vector mesons. Again, one can show that in the unitary 
gauge these fields correspond to vector matter fields of the CCWZ formulation |^ . Originally |^ it was 
proposed that the p meson was the dynamical gauge boson of the hidden local symmetry Hiocai = SU{2)v 
in the SU{2)l ® SU{2)fj/SU{2)v nonlinear chiral lagrangian. The extension to SU{3) is straightforward 
pof and incorporates the p, K*, K* and (f> — uj mesons. 

Let us briefly describe the procedure. It consists in using two new SU{3) matrix-valued fields L and 
R to build up S 

Y. = LR'' . (52) 

The chiral lagrangian in ( |29| ) is then invariant under the group SU{3)l ® SU{3)ii (8) SU(3)h 

L^QLLh^x), R^gnRh^x) (53) 

where h G SU{3)h is a local gauge transformation. The local symmetry associated to the group SU{3)h 
is called hidden because the field S belongs to the singlet representation. It should be noticed that this 
description is equivalent to the previous one by the gauge fixing L = R'^ — ^, which can be reached 
through a gauge transformation of SU{3)h- With the fields L and R we can construct two currents 

= i {L^d^L + R^d^R) (54) 

- i (Lta^L - R^d^R) (55) 

which are singlets under SU{3)l <8) SU{3)b. and transform as 

^ hV^.h'' + hd^^h^ (56) 
^ hA^h^ (57) 

under the local group SU{3)h- In the unitary gauge L = R^^ = ^ they reduce to the V and A previously 



introduced in (36) and ( 



In this notation, the transformation (35) for Ha reads 

Ha - H.hlix) (58) 
and the covariant derivative is defined as 

D^H ={d^ + V^)H . (59) 
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The octet of vector resonances (p, etc.) is introduced as the gauge muhiplet associated to the group 
SU{3)h- We put 



where p is a hermitian 3x3 matrices analogous to the one defined in equation (g8|) . This field transforms 
under the full symmetry group as 

hp^h^ + hd^h^ . (61) 

The vector particles acquire a common mass through the breaking of SU{3)l <8) SU{3)u SU{3)h to 
SU{3)v- In fact, 8 out of the 16 Goldstone bosons coming from the breaking are the light pseudoscalar 
mesons, whereas the other 8 are absorbed by the p field. 

We can now build a lagrangian describing the interactions of heavy mesons with low momentum vector 
resonances, respecting chiral and heavy quark symmetries. The new terms we have to add to (^0|) are : 

-Cp = ^ < F,,{p)F^-^ip) > [< (A.r > +« < (Vp - >] 

+ if3<HbV^{V^-p^)^^Ha>+iX<Hba''-'F^,{p)baHa> , (62) 

where F^„{p) = dfj,p^ - d^Pf_, + [p^,Pi/]. 

In the first line in (|6|) there is the kinetic term for the light vector resonances. The second term gives 
back the non linear cr-model lagrangian (p9|), as it can be seen by using the identity 

{A^A^ = -i(9^Ea^St) (63) 

plus interactions among pions and p-like particles. The value of the parameters a and gy can be fixed by 
considering the electromagnetic couplings |3^ . In this way one can see that the first KSRF relation |^ 
is automatically satisfied 

9p 9p-K7rf , (64) 
with gp is the p — 7 mixing parameter, and gpT„T — o,gv l"^- Furthermore, from the second KSRF relation, 

Wp = gl^^f (65) 

and extracting the p mass from (|& 

ml^\aglf , (66) 
we see that a = 2, and that 

9v = -^~5.8 . (67) 

The terms proportional to j3 and A give the couplings of the light vector mesons with the heavy states, 
like PPp, PP*p etc. 

As usual in (|6^) we have considered the lowest derivative terms. Explicit symmetry breaking terms 
can be introduced as in ( 44 ) and ( ^ ) . 

As shown in | |4^ , the hidden symmetry approach has an interesting limit in which an additional 
symmetry appear, the so-called vector symmetry. In this limit the vector meson octet is massless and the 
chiral symmetry is realized in an unbroken way: the longitudinal components of the vector mesons are the 
chiral partners of the pions. A chiral lagrangian for heavy mesons incorporating both heavy quark and 
vector symmetries has been written down in p3| : having an additional symmetry, there is a reduction of 
the number of effective coupling constants, and in the exact symmetry limit only one unknown coupling 
constant appears. However large symmetry breaking effects are expected and corrections to the vector 
limit can be sizeable. 
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2.5 The chiral lagrangian for the positive parity states 

In the sequel we shall use also the chiral lagrangian for the positive parity states introduced in section 



2.1. This lagrangian, containing the fields Sa and T^' as well as their interactions with the Goldstone 



bosons and the fields Ha, has been derived in refs. ||3^, |18| : 

— Ckin + CiTT + Cs + Cd (68) 

Ckin = i < Sbiv ■ D)baSa > +i <Tl^iv ■ D)baf^a > 

- 5ms < SaSa > -Smr < T^%a > (69) 
£1, = ik< T^ixi^At^a > +i~k < Sbi^j^A'iJa > (70) 
Cs = ih < Sbi^ibA-laHa > +ih < T^A^bal^Sa > +h.c. (71) 

Cd = Tl^lXl5{D^A^)baHa > 



<T^^lX-/5iD^A^)baHa> ■ (72) 



In (|69|) Srris — Mjjo ~ Mp — Mn^ — Mp, SrriT — Mp,^ — Alp = — Mp . A mixing term between the 
S and Tf^ field is absent at the leading order. Indeed, saturating the fi index of with or gives 
a vanishing result, and derivative terms are forbidden by the reparametrization invariance [^2| , [l8| . We 
can also introduce the couplings of the vector meson light resonances to the positive and negative parity 
states as follows 

£4 = Csp + Ctp + a (73) 
Csp = < 5f,t;^(V^ - Pp)baSa > +t~Xi < Sba''" Fp,{p)baSa > (74) 

Ctp = if32 < T>^(V^ - Pp)baTaX > +«A2 < T^<J^''F^,{p)bafaX > (75) 

C = tC < SaHbJpiV^ - p^ba > +IP < SaHba^" FxAp)ba > 

+ zCl < HaT,^l^{V^ - P^ba > +ipi < HaT^^j''Fp,{p)ba > ■ (76) 

We shall see in the sequel that some information on the coupling constants g, p., A, and ^ can be obtained 
by the analysis of the semileptonic decays 

H PWi, H P*m (77) 

and from the radiative decay 

P* ^ P-f . (78) 

As discussed in the next sections g and A have been also evaluated by potential models and QCD sum 
rules. 
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3 Strong interactions 



In the limit of exact chiral, heavy flavour, and spin symmetries, the low-energy interaction among two 
heavy mesons and light pseudoscalars is governed by the lagrangian (|40|). The coupling constant 
describing the coupling of the heavy mesons to the pseudoscalar Goldstone bosons, is one of the funda- 
mental parameters of the effective lagrangian. As we shall see, via chiral loops, it enters into a variety of 
corrections to both the chiral and the spin symmetry limit of many quantities of interest. 

For the time being we limit the discussion to the strong interaction among the lowest lying, negative 
parity states, (0~) and P°* (1^), contained in the multiplet H'^. Later on we shall discuss strong 
interactions involving excited states. 

The terms containing one light pseudoscalar are readily obtained from the lagrangian (|40|). They 
read: 

9s 



C = 



f 



/ 



,d^'MP^ + h.. 



(79) 



The interaction term PPtt is forbidden by parity; the direct P* Ptt transition is not allowed in the B 
system because of lack of phase space. On the other hand, this transition occurs for D mesons. From eq. 
( |79| ) one obtains the partial widths: 



T{D* 



r(P>*+ - 

The decay D 
Dtt channels 

rtot{D*+) < 



(see table 
131 KeV 



T{D* 



D\+) 



9 



■,2 



g 



is also forbidden by the phase space 



(80) 



The D*^ decay is dominated by the 
There is an experimental upper bound on the total D*^ width: 
By combining this bound with the measured branching ratios of D*^ 
reported in table |l|, one obtains the following upper limit for g: 

2 




< 0.5 



(81) 



Table 1: Experimental D* branching ratios (%) 
Decay mode Branching ratio 
63.6 ±2.8 
36.4 ±2.8 
68.1 ±1.3 
30.8 ±0.8 
1.4±0.8 

Also the radiative partial widths r(Z?*'^ D'^l) and T{D*^ D^l) depend, via chiral loops, on 



D*^ - 




D*° - 




D*+ - 




D*+ - 




D*+ - 


P>+7 



the g coupling constant ^ ^ . This dependence will be discussed in section |6T . 

A list of measurable quantities which depend on the g coupling constant either directly, or via chiral 
loop corrections, includes: the rate for B D{D*)'kIv, the form factors for the weak transitions between 
heavy and hght pseudoscalars, the chiral corrections to the ratios BbJDu+, to the Isgur-Wise 

function S^{v ■ v'), to the double ratio {fBjfB°)/{fDjfD+)i and to several mass splittings in the P^jPa 
system. The discussion of these observables will be presented below. 
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3.1 Theoretical estimates of g 

In this section we review some theoretical estimate of the strong couphng constant g defined in eq. Eol. 



3.1.1 Constituent quark models 

In the constituent quark model, one finds g — 1 pH |. As a matter of fact the axial- vector current 
associated to the lagrangian of eq. 



reads: 



Jo fl 



(82) 



where {A = 1, ...8) are SU{3) generators and the dots stand for terms containing light pseudoscalar 
fields. The matrix element of the combination j^^ + i between the D*~ and the D'^ states can be 
easily evaluated. By working in the D*~ rest frame and by selecting the longitudinal helicity, one obtains: 

{D*-\{Jl+^J|),=,\D^)=-g (83) 

On the other hand, if one identifies, without further renormalization, the partially conserved axial currents 
of eq. (^) with the corresponding currents of QCD: 

(j5^)^^''=9«(r^).b7M75g6 (84) 
one can evaluate the same matrix element within the non-relativistic constituent quark model, obtaining: 

{D*-U+^Ji),^,\D'^) = -l . (85) 
The comparison between the eq. (pSl) and (|85|) leads to: 



5=1 



(86) 



A similar argument provides gA — 5/3 for the nucleon, to be compared with the experimental result 
gA — 1.25 (analogous result was previously obtained in the constituent quark model |Q). The authors 
in ref. find a slightly different value: g ~ 0.8, obtained in a calculation considering mock mesons (see 
references therein). A similar value g ~ 0.8 is obtained in using PC AC (see also |^ and iQ). 

In ref. it has been suggested that a departure from the naive constituent quark model might arise 
as a consequence of the relativistic motion of the light antiquark q inside the heavy meson. The model 
adopted in p7| is based on a constituent quark picture of the hadrons; the strong interaction between 
the quarks is described by a QCD inspired potential and the relativistic effects due to the kinematics 
are included by considering as wave equation the Salpeter equation [ p6| (for more details see ||57|). 

In this model one finds |27|: 



1 



dk 
2^ 



E„ 



1 - 



3{Eg + TUq 



(87) 



where Eo 



+ is the light quark energy, and ip is the wave function. By considering the 



non-relativistic limit {Eq ~ m 



3> k) one obtains g — 1, because of the normalization condition 



(27r)3 



This reproduces the constituent quark model result of eq. (p6|). 

Let us now take in ( ^7|) the limit of very small light quark masses (we note that there is no restriction 
to the values of ruq in the Salpeter equation and niq = is an acceptable value). In this case, we obtain: 



(89) 
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It is worth to stress that the strong reduction of the value of g from the naive non relativistic quark 
constituent model value g = 1 (eq. |8^) to the result (|8^) has a simple explanation in the effect of the 
relativistic kinematics taken into account by the Salpeter equation. Similar results have been obtained 
in [||. 

Including finite mass effects (m„ — rud — 38 AleV; rag = 115 MeV , nic = 1452 MeV, mb = 
4890 MeV are used in this fit) one obtains the numerical results: 

g = 0.40 {D case) (90) 

g = 0.39 {B case) . (91) 

3.1.2 QCD sum rules 



The coupling constant g has also been determined within the QCD sum rule approach |9[ |2||o[ |6l) Q. 
The starting point of this approach is the QCD correlation function: 

A^{P,q)^i J dx < Tr-{q)\T{V^,ix)j5{0)\0 > e-"^''' = Aq^ + BP^ (92) 

where, considering the case of the B system, — u'jfib, js = ibjsd, P — qi + q2, q — qi — q2 and A, B 
are scalar functions of qf, , q^. 

Both A and B satisfy dispersion relations and are computed, according to the QCD sum rules method, 
in two ways: either by means of the operator product expansion (OPE), or by writing a dispersion relation 
and saturating the associated spectral function by physical hadronic states. 

The OPE can be performed in the soft-pion limit 0, for large Euclidean momenta {qf = g| — *■ 

—oo). The various contributions come from the expansion of the heavy quark propagator and of the 
vector current . This leads to a combination of matrix elements of local operators bilinear in the light 
quark fields, taken between the vacuum and the pion state. On the other hand, when considering the 
dispersion relation for the correlator of eq. (^) , the constant g enters via the contribution of the B and 
B* poles to the spectral density, through the S-matrix element: 

<Tr-{q) B°{q2)\B*-{que)>= gB-Bne^-Qt. (93) 

From the lagrangian in eq. ([79|), one immediately finds: 

2AfB , , 

gB'Bw ~~^9 ■ ^^'^> 

The QCD sum rule approach allows to estimate directly the strong amplitude of eq. (|9^), characterized 
by the coupling constant gs* Btt- This includes the full dependence on the heavy quark mass mb, not only 
its asymptotic, large rrib, behaviour. 

On the other hand, by retaining only the leading terms in the limit mb oo, one obtains the following 
numerical results from the sum rule: 

g = 0.040 ± 0.005 GeV^ . (95) 

where F parametrizes the leading term in the decay constants fs and /s» : 

/b=/b. . (96) 

^For a complete list of earlier references, see ref. ]6o| . 
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F can be computed by QCD sum rules. For example, for cu = 0.625 GeV (w is the binding energy of the 
meson, finite in the large mass limit), and the continuum threshold parameter yo in the range 1.1 — 1.4 
GeV, and neglecting QCD corrections the result is 

F = 0.30±0.05GeV3/2 . (97) 

By including radiative corrections one finds higher values (around 0.4-0.5 GeV'^^^) that are compatible 
with the results obtained by lattice QCD; 

F = 0.55 ±0.07 [63] 
F = 0.61 ±0.08 H 

F = 0.49 ±0.05 |6| . (98) 



Since one has neglected in ( |95| ) radiative corrections, a safer value for F is given in eq. (|97|), which 
is also the value we shall use in the subsequent sections. From eqs. ( |95| ) and (^^, one would obtain: 
g = 0.44 ± 0.16. 

An independent estimate of g can be obtained by expanding the correlator of eq. ( p^ near the 
light-cone in terms of non-local operators whose matrix elements define pion wave functions of increasing 
twist (this method is called light-cone sum rules). In this way, an infinite series of matrix elements 
of local operators is effectively replaced by a universal, non-perturbative, wave-function whose high- 
energy asymptotic behaviour is dictated by the approximate conformal invariance of QCD. By using this 
technique, in the following result has been obtained: g = 0.32 ± 0.02. 

Our best estimate for g, based on the analyses of both QCD sum rules ^ and relativistic quark 
model is 

g « 0.38 (99) 
with an uncertainty that we estimate around ±20%. This is the value we shall use in the next sections. 



In section S.l we will show that also the results from radiative D* decays are compatible with (M 



3.2 Chiral corrections to g 

Due to the exact chiral symmetry of the interaction terms in eq. (^o|), the coupling constant g does not 
depend on the light flavour species. Chiral breaking effects can be accounted for by adding breaking 
terms to the symmetric lagrangian. The chiral breaking parameters are the light quark masses, and the 
lowest approximation consists in keeping all the terms of the first order in the quark mass matrix. 

On the other hand, in a given process, corrections to the chiral limit can arise in two ways: either 
via chiral loops, with mesons propagating with their physical, non-vanishing mass, or via counterterms 
which affect the considered quantity at tree-level. The latter corrections exhibit an analytic dependence 
on the quark masses and are typically unknown, being related to new independent parameters of the 
chiral lagrangian. On the contrary, the former terms contain a non-analytic dependence on the quark 
masses which is calculable via a loop computation. The loop corrections in turn depend explicitly on 
an arbitrary renormalization point /i^ (e.g. the t'Hooft mass of the dimensional regularization) . This 
dependence is cancelled by the fjl^ dependence of the counterterms. 

Although the overall result is given by the sum of these two separate contribution, it is current practice 
to estimate roughly the chiral corrections by neglecting the analytic dependence and by fixing to about 1 
GeV the renormalization scale /i in the loop computation. The adopted point of view is that the overall 
effect of adding the counterterm consists in replacing /x^ in the loop corrections with the physical scale 
relevant to the problem at hand, A^. Possible finite terms in the counterterm are supposed to be small 
compared to the large chiral logarithms due to the formal enhancement of the non-analytic terms as 
m^logmj over the analytic ones. In view of these uncertainties the results of this method are more an 
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indication of the size of the corrections than a true quantitative calculation, since there are examples e.g. 
in kaon physics where the finite counterterms are not negligible ( [|66| ). 

With this philosophy in mind, the chiral corrections to the g coupling constant have been evaluated 
in ref. [|67| [6^ , |36| . Neglecting the u and d quark masses in comparison to the strange quark mass and by 
using the Gell-Mann-Okubo formula to express in terms of m\ (m^ = 4/3m^), the leading one-loop 
logarithmic corrections can be expressed in terms of 



1 



X 



16^2 f2 

-0.125 in 



1 GeV) 



The one- loop coupling constant, ffe//, is given by: 



Jeff = 9 



1 - (1 + f 5^)x' 



0.45 



(5 = 0.38) 



(100) 



(101) 



In computing these class of corrections one may use the Feynman rules reported in appendix A. The 
result ( |101| ) will be used in the evaluation of the loop corrections to the matrix element of the weak 
current between tt and P {P = D, B), as discussed in section 5.2.2. 



3.3 Hyperfine splitting 

As an example of application of the chiral perturbation theory to the calculation of physical observables 
relative to heavy Qq mesons, in this section we work out in some detail the hyperfine mass splitting 
between 1~ and 0~ mesons. As a matter of fact, the spectroscopy of heavy mesons is probably the 
simplest framework where the ideas and the methods of heavy quark expansion can be quantitatively 
tested. As explained in section |]^, the splitting among the 1^ and 0~ heavy mesons masses is due, at 
the leading order, by the 1/M correction of eq. (2.48): 

A = {mp, ~mp) = -^ . (102) 

The exp erimental data, listed in table |[ supports quite well the approximate scaling law suggested by 
eq. ( |l02| ). These data can be used to estimate the parameter A2: 



Table 2: Experimental mass splittings between 1 and mesons. 



A (MeV) 


Md'+ 


- Md+ 


140.64 ±0.09 




- Mdo 


142.12 ± 0.07 






141.6 ± 1.8 


Mb- - 


-Mb 


46.0 ±0.6 


Mb'^ - 


-Mb^ 


47.0 ±2.6 



A2~ 0.10 -0.11 GeV'^ . (103) 

The second term in eq. (^), independent of the heavy quark flavour, is responsible for the mass splitting 
between strange and non-strange heavy mesons: 

= 2A1TO5 . (104) 

Experimentally one has psf : 

M^.± - Md± = 99.1 ± 0.6 MeV 

Mbo - Mb = 96 ± 6 MeV , (105) 
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leading to 

Ai ~ 0.33 . (106) 

Recently, attention has been focused on the combinations |6^. |7^, |7|, 36 : 

Ad - {Mot - Md, ) - iMD'+ - ) (107) 
As = (Mi3j - Ms, ) - [Mb.o - Mbo) (108) 

which are measured to be |p5|| : 

Ad ^1.0±1.8 MeV (109) 

As ~ 1.0 ±2.7 MeV . (110) 

This hyperfine splitting is free fr'om electromagnetic corrections and vanishes separately in the SU (3) 
chiral limit and in the heavy quark limit. In the combined chiral and heavy quark expansion, the leading 
contribution is of order ms/rnq and one would expect the relation p9t : 

Ab = —Ad . (Ill) 
rrib 

In our framework the lowest order operator contributing to Ad.b is: 

r,0, = ^Tr[H^a,^H,an^-^ . (112) 
The matrix is 

where m is the light quarks mass matrix and ^ the coset variable defined in eq. (2.33). By taking 
rUs/A^ ~ 0.15 and rj ~ Ag^^^ ~ 0.1 GeV^ one would estimate: 

Ad ^ 20 MeV (114) 

As ~ 6 MeV . (115) 

Given the present experimental accuracy, the above estimate is at most acceptable, as an order of mag- 
nitude, for Abi while it clearly fails to reproduce the data for Ad- If the contribution from O2 were the 
only one responsible for the hyperfine splittings, agreement with the data would clearly require a much 
smaller value for rj. 

In chiral perturbation theory, an independent contribution arises from one-loop corrections to the 
heavy meson self energies [ |n| , evaluated from an initial lagrangian containing, at the lowest order, both 
the chiral breaking and the spin breaking terms of eq. (^4|) and (|4^). These corrections can be computed 
by using the Feynman rules given in appendix A; they depend on an arbitrary renormalization point /i^ 
(e.g. the t'Hooft mass of dimensional regularization) . This depe nden ce is cancelled by the /i^ dependence 



of the counterterm r]{pL )02- Following the discussion in section |3.2| one can use /i = 1 GeV 



The possible sources of hyperfine splittings via chiral loops are the light pseudoscalar masses TOtt , uik 



and rrirf, the mass splittings As, A of eqs. (|102| ), (104) and, finally, the difference between the P*P*-k 
and the P*P7r couplings [P ^ D,B) induced by the last term of eq. (^). 
This splitting is of order 1/M, and, from eq. (IsJ), one obtains: 



Ag = 5p.p.^ - gp.PTT = 2 H . (116) 
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The second term in (^8|), proportional to gi, breaks only the heavy flavour symmetry, making the B*B^*\ 
and D* D^*^t: couplings different. The third term, proportional to 52, breaks also the spin symmetry and 
contributes differently to the P*Ptt and to the P*P*tt couplings. This is precisely the effect relevant to 
the hyperfine splitting. 

In terms of these quantities, one finds 0, 

f ^^2 i„„/ X.) 2TO„l0g(^ 



4m2,log(- 

'"'K 



A2 

6m2log(^) 



5 



6np g 



1 3 
2< 



(117) 



The dependence upon t he he avy flavour P = D, B is contained in the parameters A and A^. 

The first term in eq. (117) is the so called chiral logarithm |71 . In the ideal situation with pseudoscalar 
masses much smaller than A^, it would represent the dominant contribution to Ap. Calling A^ and A^ 
its value for the D and B mesons, respectively, one finds: 



A" 



-13 MeV, A^g ~ +4 MeV, 
where we are using the representative value g = 0.38 (see eq. (|9 



The second term in eq. (117) represents a non analytic contribution of order i 



3/2 



(118) 



I , which, although 



formally suppressed with respect to the leading one, is numerically more important, because of the large 
coefficient 2Att. The separate contributions to the D and B hyperfine splittings read: 



-30 MeV, 



-9 MeV 



(119) 



Finally, the last term in eq. (117) is also of order m'J'^ . Its evaluation requires the estimate of 
the difference Ag/g, which is not directly related to other experimental data. In ref. ||7^, this difference 
has been computed in the framework of QCD sum rules. Using mf, — 4.6 GeV and rric — 1-34 GeV one 
gets: 

fl> 9B-B-n = 0.0094 ±0.0018 06^2 

fl. go'D-n = 0.017 ± 0.004 GeF^ . (^^20) 
and for the gp^p-rc coupling 

fBfB'9B'B. - 0.0074 ± 0.0014 Gey2 

fo fo- go'DTT = 0.0112 ± 0.0030 Gey2 . (121) 

To derive the difference Ag at first order in 1/A/, one should expand the relevant sum rules in the 
parameter 1/Af, keeping the leading term and the first order corrections which are given by 



gi +92 
9 



2A' 



and 



2^ + {A' -A) 



-0.15 ±0.20 GeV 



0.99 ±0.02 GeV 



91 - 92 
9 



A' + A. 



-1.15 ±0.20 GeV 



(122) 



(123) 



The couplings gi and g2 have been defined in eq. ( pT] ) and the parameters A and A' are related to the 
1/M corrections to the leptonic decay constants, fp and /p.: 



F 



A 
M 



If 



F 

Jm 



(124) 
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Neglecting radiative corrections, A and A' are given by 



A = -L 



Gk 



3Gy 



A' = -- 



Gk 



-Gy 



(125) 



where lu represents the difference between the pseudoscalar meson and the heavy quark masses, at leading 
order in l/M. The sphtting of the couplings depends on the quantity 2 g2/g that contains only the 
difference A' — A given by: 



A' -A = -Lj- 4Gs 
3 



(126) 



There is disagreement in the literature on the values of the parameter G-^: at the b quark mass scale 
from ref.(73| one gets Gs = (0.042 ± 0.034 ± 0.023 ± 0.030) GeV, while in ref. [|2| the central value 
Ge — —(0.052) GeV is quoted. In view of this discrepancy, to provide an estimate of the difference (126), 
it is reasonable to approximate A' — A^i 2/3(jj « 0.4 GeV, obtaining 



2^ 
9 

From (111 
Ab 



0.6 GeV 



(127) and from the formula (117) of the hyperfine mass splitting one finds: 



^ (27.3 + 61.4 - 75.8) MeV = 12.9 MeV 



(127) 



(128) 



We notice that we have used in eq. (117) / = = 130 MeV for all the light pseudoscalar mesons of 
the octet. In eq. ( |128D we have detailed the contributions A°, and the one from Ag/g respectively. 



We have also taken A^^ = 1 GeV. It is evident that there is a large cancellation among the last term and 
the other ones. For the value g = 0.38 we obtain 



Ab ~ 1.9 MeV 



(129) 



The application of this result to the charm case is more doubtful, in view of the large values of the 1 /nic 
correction Ag/g. By scaling the result (129) to the charm case, one would obtain 



A 



nib 



D 



Ab =i 6.3 MeV 



(130) 



In conclusion we observe that the application of chiral perturbation theory to the calculation of the 
heavy meson hyperfine mass splitting is rather successful, even though, given the large cancellations in 



eq. (128), the results (129) and (130) should be considered as order of magnitude estimates only. 



3.4 Strong decays of positive parity states 

In this section we shall examine the applications of the effective lagrangian approach to the strong decays 
of the positive parity heavy meson states. We shall first review the experimental evidence for these states; 
next we shall give the formulas for the decay rates into final states with one pion. Finally we shall present 
some estimates of the couplings based on QCD sum rules and we shall apply them to the calculation of 
the strong decay rates. 

Strong transitions of the positive parity states, contained in the multiplets S and T introduced in 



m section 



section 2.1, are described in the present formalism by the lagrangian £3 of eq. (68), explicitly discussed 
2.5. The experimental data concerning these states are still at a preliminary stage. In the 



charm sector, the total widths of the 1)2 (2460) and _Di(2420) states, have been measured 

rfot(I?2(2460)) = 21 ±5 MeV 
TtotiD i{2A20)) = 18 ±5 MeV . 



(131) 
(132) 
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As for the B sector, evidence has been recently reported [[74| [T^] of a bunch of positive parity states i?**, 
with an average mass 

TO^.. = 5732±5±20 MeV (133) 
and an average width 

T{B**) ^ 145 ± 28 MeV . (134) 
The OPAL cohaboration of LEP has also reported evidence of a B** state with mass 

TOB*. = 5853 ± 15 MeV (135) 
and width 

T{b;*) = 47 ± 22 MeV . (136) 

The decay widths of the states (l''",0'''), belonging to multiplet sf = (1/2) + , here referred to as Pi 
and Pq, are expected to be saturated by the single pion channels ||l^, |l^: Pq Pn and Pi P*n. 
Therefore these transitions are controlled by the coupling constant h of eq. (2.71). In the mg — > oo limit 
one obtains 

1 /h^' 



T{Po^P+n-)^r{Pi^P*+n-) = ^i^jj {dm,)^ (137) 

where Snis is the mass splitting of the states S with respect to the ground state H. From estimates based 
on quark model ||7^, [T^ and QCD sum rules [^8[ ^ computations of the masses of these states, one has 
A = 500 ± 100 MeV. We notice that this mass splitting agrees rather well with the experimental result 
in the B sector, given in eq. (|133|). 



On the other hand, the formula (137) is of limited significance, especially for the case of charm, due to 



the large 1/M corrections coming from the kinematical factors. Keeping M finite, the formulas become: 
r(Po - P - ) = ^Gp..p, ^ . (138) 



r(Pi ^ P*+7r 



87r 2M|,^ ^ 

For the B system the coupling constants Gp<-*pt^ and Gp-^p^T^ are defined by the strong amplitudes: 

Gb..b„ =<7r+(g)P°(g2)|P**+((7i)> (140) 

Gb^B'^ = < n+iq) B°*{q2)\B+{qi) > (141) 

where B** and Bi denote the 0+ and the 1+ states in the = (1/2)+ doublet. Analogous definitions 
are understood for the D's. In the infinite mass limit, Gp-'p-n and Gp^p^jr coincide. The amplitude 
Gb'^Bit is related to the strong coupling constant h appearing in the heavy-light chiral lagrangian ( [7l| ) 
by the formula: 

r2 II f 2 



Mf,,. - Mi h 



Gb"b. = - VMbMb.' ^ — . (142) 

Mb" Jtt 
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In the limit mi, — > oo one has: 

■ 1 



Mb = mb + u} + 



rrib 



Mb- -Mb = Snis + O ( — ] (143) 

\mbj 

2h 

and one recovers eq. ( |l37[ ). 

Differently from the decays of the positive parity states having sf = (1/2)+, the single pion transitions 
of the sf = (3/2)+ particles, here denoted P2 and P{, occur with the final pion in D-wave. The decay 
rates for these transitions are given by ||l^, p^ : 



r(po - 


P+TT-) 


r(p° - 


^ P*+7r-) 


r(pf - 


^P* + TT-) 



1 Mp h!^ \p^f> 
ISttMp, A^^ 

1 Mp, h'^ \p^f 
IOtt Mp, A^^ 
1 Mp. \p^f 



(145) 
(146) 
(147) 



Gn Mp, A/ 

where the strong coupling h' is given by: 

h' = hi+h2 , (148) 

in terms of the parameters hi and /i2 of eq. (|7^). From the previous equations, one finds the following 
prediction in the D system: 

= 2.7 (149) 

in good agreement with the experimental result Ea, 2.4 ± 0.7. 



To get numerical results for the rates given in cqs. ( |138| ), (139) and (145-147), one should specify the 
relevant coupling constants. The parameter Gp.^p^ has been evaluated in the framework of QCD sum 
rules, by means of two independent methods [|oj. The first method, based on the single Borel transform 
of an appropriate correlator evaluated in the soft pion limit, gives the results: 

Gij. . = 13.3 ±4.8 GeT/ (150) 

Gi)..c^ = 11.5 ±4.0 GeT/ . (151) 

We observe substantial violations of the scaling law Gd-^Dtt/Gb'-'Btt ~ Tnc/mf,. 

In the limit M — > 00 we obtain, from the asymptotic {mi, — s- 00) limit of the sum rule: 

h= -0.52 ±0.17. (152) 

The second method is based on the light-cone sum rules |8^, |6^, |8^ . One obtains |Q : 

Gi3..B^ = 21 ± 7 GeV (153) 

Gp,..£,^ = 6.3± 1.2 GeV (154) 
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A two-parameter fit of the above results in tlie form 
h(rn) 



h (1 + -) 

TO 



gives for h (see eq. ( |142D the result: 
h = -0.56 ±0.28 



(155) 



(156) 



and for the parameter tr, ct = 0.4 ± 0.8 GeV. 

The values of h found by the two methods agree with each other. As for the finite mass results, the 
two methods sensibly differ (almost a factor of 2) in the case of the charm, while the deviation is less 
important for the case of beauty (around 40%) . These differences should be attributed to corrections to 
the soft pion limit that have been accounted for by the sum rule based on a light-cone expansion. 

Using Go'Dn = 6.3 ± 1.2 GeV, Gb-b^ = 21 ±7 GeV, and Ad = Ab = 500 MeV, from eq. (|l3|) 
one finds 



r(-Do 



Dtt) 
Bit) 



180 MeV 
360 MeV 



(157) 
(158) 



There is no direct information on the coupling Gp-^p'Tr- In the infinite-mass limit it coincides with 
GpfpTT and, in order to estimate the widths of the l"*" states, we assume that this equality holds for 
finite mass as well. From eq.(139) we obtain: 



r{Di D*tt) ~ 165 MeV 
T{Bi B*tt) ~ 360 MeV . 

Also in this case we have taken Mp^—Mp* 



(159) 
(160) 

500 MeV (P = B, D) as suggested by HQET considerations. 



To estimate the strong coupling constant h! j , one can make use of the total decay width of 
eq. (131), Ftot (1)2(2460)) = 21 ± 5 MeV. Assuming that only two body decays are relevant, one gets 
/I'/A^ K, 0.55 GeV~^. From this result and from eq. ( |l47|) one obtains for the state D\ the total width 
Ttot ~ 6 MeV to be compared with the experimental width of t he other narrow state observed in the 
charm sector, Ftot(£'i(2420)) = 18±5 MeV, also given in eq. (|l32| ). This discrepancy could be attributed 
to a mixing between the D\ and the D\ states ||84| . If a is the mixing angle, we have 



sin^(a) 



12 MeV 



0.08 



(161) 



and therefore one gets the estimate a ~ 16° [gO|. This determination agrees with the result of Kilian et 
al. in ref. p8[ . 

In ref. |l7|| the decay rates for the transitions of the (l+,2+) states with the emission of two pions, 
TTTT, have also been estimated. They appear to be suppressed with respect to the single pion 

rates. 

In the B sector, the recently observed positive parity states B** , whose average mass and width are 
given in eqs. (133) and (134), can be identified with the two doublets (2+, 1+) and (1+, 0+). To compare 
previous estimates with the data, we average the widths of the (l+,0+) multiplet, eqs. (158) and ( |160| ), 
with those of the 2+ and 1+ states, obtained from eq. (145- 147| ): 



rtot{B2) 

TtotiBi) 



12 MeV 
10 MeV 



(162) 



27 



It is difficult to perform a detailed comparison of these results with the yet incomplete experimental 
outcome. However, assuming that the result obtained by LEP collaborations in the B system represents 
an average of several states, and neglecting a possible mixing between the states Bi and Bi (a 1/M 



effect), the experimental width is compatible with the previous estimate. 

Finally, the total width in eq. ( |136| ) can be interpreted as connected to the decay B** BK, B*K. 
Assuming again that the width is saturated by two-particle final states, and using Mb" — 5853 MeV, 
we obtain: 

r{B**{0+)) ~ 280 MeV (163) 

T{B**^B*K) ~ 200 MeV (s^)^ = (1/2)+ (164) 

r{B**{l+)) ~ 0.45 MeF (s^)^ = (3/2)+ (165) 

T{B;*{2+)) ~ lAMeV . (166) 



Also in this case a detailed comparison with the experimental results cannot be performed without more 
precise measurements; we observe, however, that the computed widths of the different i?** states are 
generally smaller than the corresponding quantities of the B** particles, a feature which is reproduced 
by the experiment. 

4 B ^ D decays and chiral dynamics 

One of the most important applications of the heavy quark symmetry is the analysis of the exclusive 
semileptonic decays B Dlvi and B D*lvi. We shall here give a brief summary of this extensively 
studied subject: for more details see for instance ||^ and references therein. 

In the symmetry limit, i.e. infinite D and B masses, the six form factors generally needed to param- 
eterize the matrix elements < D^*\v')\J ^\B{v) > {v, v' velocities) reduce to a single function ^(w • w'), 
the Isgur-Wise function. One finds [Q: 

< D{v')\c-i,MB{v) > - ^MeMD^iv ■ v'){v + v')^ 
<D*iv',e)\cj^b\B{v)> - ^/M^M^i^ivv')e^,^pe*''v'"v^ 
<D*iv',e)\cj,,-f5b\B{v)> = ^MBMoiiv ■ v') [(! + «• v')el ~ {e* ■ v)v'^] . (167) 

Various calculations of the Isgur-Wise function exist in the literature; they use different non-perturbative 
approaches, such as QCD sum rules |8^ or lattice QCD A review of these results would be outside 
the scope of the present report, and we refer the interested reader to the literature. 

At the symmetry point, i.e. v = u', the normalization of the Isgur-Wise function is known: ^(1) — 1. 
This is a consequence of the conservation of the vector current = h'^^^hy = h'^v^h^ and allows a model 
independent determination of the CKM matrix element Vcb from semileptonic heavy to heavy decays by 
extrapolating the lepton spectrum to the endpoint v = v' . 

Of special interest for this determination is the decay B — > D*lv, since there are no l/vriQ corrections 
for the axial form factor Ai, dominating the decay rate, at the symmetry point. This is the content of the 
Luke's theorem A simple proof of this important result has been presented by Lebed and Suzuki 1^. 
Luke's theorem is an extension to the spin-flavour symmetry of the AdemoUo-Gatto theorem [^8| , which 
was originally stated for the SU{i) flavour symmetry of light quarks and refers to the matrix element of 
the vector current between states belonging to the same SU{i) multiplet at = 0. The statement is that 
matrix elements of a charge operator, i.e. a generator of the symmetry, can deviate from their symmetry 
values only for corrections of the second order in symmetry breaking. In the case of semileptonic decays 
B D{D*)li>, the only form factor protected by this theorem against l/mg corrections at the symmetry 
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point V = v' is Ai, dominating the decay B D*lv 'At v — v' . In practice, Luke's theorem reduces to 
the result: 

Mb + Md- , , 2 



(^max) ^VA+ <5l/m2 (168) 



where ija = 1 if strong radiative corrections are neglected. 

The l/rriQ corrections at the point v = v' have been estimated a combined analysis |^ gives 
a correction to ■^(l): Si/m^ = —(5.5 ± 2.5)%. Also leading and subleading QCD corrections arising from 
virtual gluon exchange have been computed: see for instance |8|. 

Experimental measurements close to v = v' suffer of large errors, due to the smallness of the phase 
space: high statistics is needed to reduce the uncertainty in the extrapolation of the lepton spectrum 
to this point: nevertheless, the exclusive semileptonic decay B D*lvi can provide a rather precise 
measurement of the element Vcb of the CKM matrix, complementary to the analysis of the inclusive 
semileptonic decay rate. 

The HQET has also been used to investigate the semileptonic decay of a B meson into an excited 
charm meson @, where S£ is the total angular momentum of the light degrees of freedom and I 

the corresponding orbital angular momentum of the charm meson (s£ = Z ± 1/2). At the leading order, 
the matrix element 

<D^'^^'\v')\J,\B{v)> (169) 

appearing in the semileptonic transition is described by a single form factor ^(**^')(ii . v'): the Isgur-Wise 
function for the B D^*' transitions is the function ^(^/^'O) = ^; for the P-wave heavy mesons we have 
^(3/24) ^ ^^^^ and ^(^/^^^) = T1/2; they have been computed by QCD sum rules in and by constituent 
quark models in |77|] . 

4.1 Chiral corrections 

Violations to SU (3) symmetry can be computed by means of the effective heavy meson chiral lagrangian. 

To estimate the size of the chiral corrections, it is common practice, as we have stressed already, to 
retain only the non-analytic terms arising from chiral loops. Moreover when the subtraction scale /i is of 
order of the chiral symmetry breaking scale A^ « 1 GeV , the coefficients of the higher order terms do 
not contain large logarithms, and therefore the numerical estimates are carried out at this scale. 

Chiral perturbation theory has been used to compute the leading corrections to the form factors for 
B D{D*) semileptonic decays, arising from the chiral loops of figure |l|. The dominant corrections at 
zero recoil, i.e. v — v' , are of special interest and have been computed in |9l| According to Luke's 
theorem, these corrections appear at the order l/mg. This class of corrections should not be confused 
with those coming from the I/^tt-q terms present in the effective lagrangian or in the current: the l/mg 
terms in the lagrangian, in particular the one responsible for the hyperfine mass splitting P* — P and 
the one giving the splitting between the couplings go and gB, generate at one-loop l/niQ non- analytic 
corrections. The effect of the gB and go splitting has been neglected in For instance, the 

B — > matrix element at the recoil point v = v' , as computed by the formulas of appendix A, is 

< D{v)\4''\Biv) > - 2v^[l + Ci^i)/ml+ 

.2 / A \ 2 



2 {AnfJ 



[/(A,/m.)-flog(MV™')] 
<D*{v,e)\4'\B{v)> = 2e;[l + C'{^,)/ml+ 

(170) 



9l{^\ [/(-A,/mO + log(/.Vm2)] 
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Figure 1: Diagrams for the calculation of one- loop chiral corrections to the B D^*^ transition matrix 
element. The box represents the b ^ c weak current, the dashed line a light pseudoscalar. 



where C and C stand for tree level counter-terms and 



dz- 



1 



(z2 + 1)3/2 \^ [(^2 + 1)1/2 



1 



(171) 



In ( |17C| ) only the dependence on Ac — Mo' — Mo has been kept, discarding the A;, terms and those 
proportional to the gs — go splitting. Numerically, for /x = 1 GeV and g = 0.38, the correction from the 
logarithmically enhanced term in (170) is —0.6%, and the correction from f{x) is 0.3%. 

A complete calculation of the l/mq and SU{3) breaking corrections to the Ba D^hvi process 
has been performed in [ p3| . This analysis includes non-analytic terms arising from chiral loops and the 
analytic counterterms, but it lacks predictive power due to the introduction of many unknown effective 
parameters. 

In the SU (3) limit, the Isgur-Wise function is independent of light quark flavor of the initial and final 
mesons, i.e. 

e« = ed = 6 , (172) 

where £,u,d,s is the Isgur-W ise f unction occurring respectively in Bu,d,s decays. In ^ th e leading 
corrections to the equality (172) have been computed in chiral perturbation theory, giving |9l[| 



g^njv ■ v') 



[rriK log (ml- ^2 



log (^2/^2) „ 3 2 (^2/^2 



(173) 



where 



2 + x 



2 



log 



+ 1 + y/jp 



4i/e2~ 



log 



\fx^ 



(174) 
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In (173) the analytic counterterms are neglected. Numerically, the nonanalytic chiral correction is a 
few percent. 

We mention here another calculation in the framework of the heavy meson chiral perturbation theory, 
the ratio of the parameters Bb^ and Bb, entering in the analysis of B^s) ~ ^(s) mixing and defined as: 



<Biv)\h^il~j5)dbr{l~j5)d\B{v)> = ^flBB 
<B,(«)|67^(l-75)s 67^(1 -75)s|i?.(«)> = Ifl^BB, 



(175) 
(176) 



In the chiral symmetry limit Bb^/Bb — 1. For non-zero strange quark mass, the ratio is no longer 
equal to 1, and the one loop chiral corrections, arising from the diagrams of figure 0, are ]95||: 



^ ^ l-|^K.o.(..4/M') + 

1 3 
+ 2^"' (w'/m') - log (ml/fi^) 



(177) 



Using = 1 GeV and g ~ 0.38 the previous formula gives BbJBb — 1.03. 



Figure 2: Diagrams appearing in the calculation of one-loop chiral corrections to the B — B mixing. The 
self-energy diagrams are not shown. The dot represents the /S.B — 2 operator 



4.2 The B D'-*hlui decay 

Another application of the chiral lagrangian can be found in the semileptonic decays of B into a charmed 
meson with the emission of a single soft pion, i.e. B D^*'>ttIvi. The phenomenological heavy-to-heavy 
leading current 1^ : 

Jf = -i{v ■ v') < Hi^h,{l - l^)Hi'^ > (178) 

does not depend on the pion field, and therefore the amplitude with emission of a single pion is dominated 
by pole diagrams, where the pion is emitted by the initial B or the final D^*\ and is proportional to 
the coupling g. The two diagrams are shown in figure 0. These decays might be used to determine the 




Figure 3: Pole diagrams for the B —f D^*^-!tIi/i decay. The square represents the b c current, the dashed 
line is the pion 
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value of or even to test the heavy quark flavour symmetry prediction for the D*Dtt and B*BTr vertices 
93 = go = 5- Moreover these processes may give indications on resonance effects. 

The chiral calculation is reliable only in the kinematical region of soft pions. In the decay B — > Dnlvi, 
the soft pion domain is a large fraction due to the inclusion of the cascade decay B D*lvi — s- Dirlvi. 
This process has been trea ted by va rious authors ^ and the analysis has been extended to 
B -f D*t:Ivi in |||, |^ |lO§: in flO^, in addition to the ground state mesons D, D* , B and _B*, also the 



contribution of the low-lying positive parity 0+ and 1+ resonances and some radially excited states is 
estimated. 

4.3 The heavy-to-light effective current 

The weak current for the transition from a heavy to a light quark, Q g^, is given at the quark level 
by 9a7/i(l — IbjQ'i when written in terms of a heavy meson and light pseudoscalars JTo| , it assumes the 
form, at the lowest order in the light meson derivatives, 

K = %<l'{^-l5)Hbila> ■ (179) 

This operator transforms as (3l,1_r) under SU{i)L x ^[/(S);^, i.e. analogously to the quark weak 
current, and is uniquely defined at this order in the chiral expansion. 
From the definition of decay constant of a heavy meson P 

< 0\qaY'l5Q\Pb{p) ip''fpjab (180) 

one gets 

fp^^^= . (181) 



We note that in the infinite quark mass limit, Mp^ — > mg, and there is no dependence on the light 
flavour. The previous formula shows the 1/ ^mq scaling of the heavy meson leptonic decay constant in 
the mq oo limit, and its light-flavour independence in the chiral limit (we neglect the small logarithmic 



dependence of F on mq). In section 3.1.2 we have already discussed the various determinations of F, see 
eqs. ©, ® . 

Higher derivative, spin breaking, and SU{3) breaking current operators are written explicitly i n |35|| : 
their introduction adds many unknown effective parameters, and they correct the leading behaviour ( |l8l| ). 
Lattice calculation [ |lOl| and QCD sum rules j6^, ^ indicate that the 1/toq corrections are sizeable at 
least for fjj. 

The current describing weak interactions between pseudoscalar Goldstone bosons and the positive 
parity S fields is introduced in a similar way: 

L^a^'-^<r{l-l5)Sbil> (182) 

The analysis done in j?^, based on QCD sum rules, gives for F'^: 

F+ ~ 0.46 Gey3/2 ^ ^^gg^ 

The current describing the interaction of the H fields with the light vector mesons, is, at the lowest 
order in the derivatives: 

= a,<jP{l-j5)H,Vp{p>'^VnUL> + 

+ a3<j''il-j5)HhV^ipp-Vp)Ul> . (184) 



32 



The current (184) is of the next order as compared to the currents (179) a nd (|18^ ), and does not 
contribute to the leptonic decay constant fp. As we shall see below, the term in (184) proportional to ai 
contributes in a leading way to the Ai form factor in the P p semileptonic matrix element, while the 
terms proportional to a2 and as contribute to the A2 form factors, but they are subleading with respect 
to the pole diagram contribution. 

We observe that there is no similar coupling between the fields T^, defined in ( p^ ) and ^. Indeed 
( |179D and (182) also describe the matrix element between the meson and the vacuum, and this coupling 
vanishes for the 1+ and 2+ states having si — 3/2. This can be proved explicitly by considering the 
current matrix element (A'' — qal^lzQ)'- 



<{)\A^^\b^>= fe^ , 

where Di is the 1"*" partner in the s; — 3/2 multiplet. Using the heavy quark spin symmetry, ( |18 
out to be proportional to the matrix element of the vector current between the vacuum and the 2 
which vanishes. 



(185) 

turns 
state, 



4.4 Chiral corrections for fpjfp 

In the chiral limit, the leptonic decay constant docs not depend on the light flavour, i.e. 
IP 



(186) 



As discussed in section (5.2), one can obtain an estimate of the SU{3) violations by computing the 
non-analytic terms arising from the chiral loops. 

The one-loop diagrams contributions to the leptonic decay constant fp are shown in fig. ^, and have 
been computed in keeping only the "log-enhanced" terms of the form m?' \og{rn? / p?') . For the 

ratio fus/fD one has: 



Id 



-?.log(^) + ^™^log(^ 



3 2 1 /™ 
-m^ log(^ 



(1 + 35^) 



(187) 



The corrections proportional to arise from the self-energy diagrams, fig. while the diagram ^c 
gives the ^-independent corrections. The diagram |^d, linear in g, vanishes at the leading order in I/toq. 
Using g ~ 0.38 in (187), one gets foJfD — l-H- The excited positive parity heavy mesons contribute 
to 5C/(3) violating effects as virtual intermediate states in chiral loops. In Ref. |102| the "log-enhanced" 
terms due to these excited-state loops have been computed: some of them are proportional to and 
others depend linearly on /i, h being the coupling of the vertex P**P'k. In |102| it has been pointed out 



that these terms could be numerically relevant and could invalidate the chiral estimate based only on the 
states D and D*] as discussed in sectio n 3.4 the coupling h is estimated by QCD sum rules in |80j, with 
the result h ~ —0.5, see eqs. (152, 156), where a more accurate chiral computation of the ratio foj fo 
is performed. We present here some details of the calculation: the vertices and the integrals needed for 
the loop integration can be found respectively in appendices A and B. 

The self-energy diagrams Wp give the following wave function rcnormalization factors: 



Zd = l- 



35' 



16^2/2 



3/2Ci(Ad'D, A_D.£), m^r) -I- Ci{Ad*d, Ad'^d^itik) + ^Ci(Ad.d, Ap-u, m,,) 



167r2/2 
1 



[3/2C(Ap„p, Ap„p, m^) + C(Ap„^p, Ap„^p, mA-)- 



6 



C(ApoP, ApoP,TO^) 



(188) 
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(a) 



(b) 



(c) 



(d) 



Figure 4: Diagrams for the calculation of one-loop chiral corrections to the heavy meson leptonic decay 
constant fp. The box represents the b q weak current, the dashed line a light pseudoscalar. 



3.9^ 



16^2/2 



167r2/2 



2C{ApgP^,ApgP^,mK) + -C(ApnP, ApoP,TO^) 



(189) 



where the mass splittings Ap.p — Mp» — Mp, Ap*p^ = Mp* — Mp^, and Ap*p — Mp* — Mp are 
©(I/toq), while the mass splittings Ap^p = A/p, — Mp, Apg^p = Mpg^ — Mp, and Ap,p^ = Mpg — Mp^ 
between excited and ground states are finite in the limit toq —^ oo. 

The functions Ci and C come from the loop integration and are defined in appendix B (here we use 
A = 0). 

The diagram ^ gives the same contribution as in ( 187 ) , while the diagram ||d is linear in h (the 
analogous term proportional to g vanishes), and proportional to F^: combining all the diagrams one 
obtains [pOl: 



ID 



F 



VMd 

327r2/.2 
F+ h 
"Fl6^ 



32^2/2 



log(Tf ) + log(^) + log(^) 



3 1 ' 

-Ci(A£,.£), Ap,'D,n^;r) + Ci{AD*D,AD'D,mK) + -Ci{AD'D,AD'D,m,^) 

3 1 
-C(ApoP, ApoP,m^) + C{Ap„^p,Apg^p,mK) + -C(Ap,p, Ap,p, to,,) 



3 1 

-C(Ap„p, 0, TO^) + C(Ap„^p, 0, tuk) + -C(Ap„p, 0, to„) 



(190) 



Id. = 



F 



/M 



D 



1 - 



1 



32^2/2 



2TO2,log(^) + iTO2log(4) 
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327r2/2 

327r2/2 

F+ h 



2C(Ap„p, , Ap„p^ , m/f ) + -C(Ap„p, Ap,p, m^) 



-t TrT72 2C(Ap„p^ , 0, mK) + -C{Ap,p,Q, m^) 
F Idtt^/^ 3 



(191) 



From the previous formulas, using Ap^p = 0.5 GeV, ^ 1, F+ ^ 0.46 GeV^^'^ and F = 0.30 GeV^/'^, 
one gets numerically: 



fD = 



/A7 



D 



[1 + 0.09 + 0.003g2 - 0.33/1^ - l.oo/i) 
1^1 + 0.17 + 0.59.g2 - 0.66/1^ - 1.15/i) . 



(192) 
(193) 



In the p revio us fo rmulas we have kept only the leading order in the l/mg, i.e. we have put Ad*d — 
in ( p^ and ( |l9l| ). 

It is found that the terms 0{h^) and 0{h), while important, tend to cancel o ut in (192, 193) and that 
the ratio of leptonic decay constants is numerically the same as obtained from (187): 



fD 



1.10 



(194) 



These values are obtained by using g — 0.38 and h = —0.5. 

The formula ( |187D is valid at the leading order in l/mg, and in this limit it is the same for B and D 
systems. In other terms, the double ratio Ri 



Ri = 



fs, I Jb 



(195) 



is equal to 1 in the chiral limit and in the heavy quark limit, separately. To see how Ri deviates from 
unity one has to take into account the 1/M terms in the chiral effective lagrangian and in the effective 
current. As discussed in four new parameters contribute at the order l/mg to the leptonic decay 
constants: two of them, pi and p2-, come from the l/mg terms in the current as 



f (l + |^)<7''(l-75)i?.4> + 

and they modify the leptonic decay constants as follows 

Pi + 2p2 ' 



(196) 



/Mp/p 
/Mp/p. 



Mp 

Pi - 2/92 

Mp 



(197) 



The t wo parameters pi and p2 can be related to the HQET matrix elements Gk and Gs defined in 
eq. (|125|), and estimated by QCD sum rules in [73|. 
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The other two couphngs, gi and 32, have been aheady introduced in (48), and they parameterize the 
1/M corrections to the couphngs gp^pj; and gp*p*.n'- 

gp-P'TT = 5+7^(51+52) (198) 



M 



p 



gp'PTT = 5 + Ti~C9i ^ 52) ■ (199) 



M 



p 



For the chiral correction to the double ratio neglecting as usual the analytic counterterms, only 
the quantity 171 — 52, i.e. the 1/A/p correction to gp.p^, is relevant, and one gets ||35|| 

i?i - 1 = -0.11(/2 -0.06 3(51 -52) GeT/'i . (200) 

5 Heavy-to-light semileptonic exclusive decays 

Most of the known CKM matrix elements have been determined using semileptonic decays. In particular, 
from semileptonic B decays one can extract Vet and Vut- 

The extraction of the value of Vcb from the exclusive process B — > D*lvi has been studied in the 
HQET context, and we mentioned it before. The situation for Vub is apparently more uncertain, both for 
the inclusive and the exclusive semileptonic rates. Its determination is one of the most important goals 
in B physics, but it involves great experimental and theoretical difficulties. At the moment there is a 
safe experimental evidence for b u transitions, and experimental data on the B —>■ nlvi and B — > plvi 



exclusive processes have been presented by the CLEO II collaboration |103|. 

The interpretation of the inclusive b u semileptonic rate is difficult because of the dominant b c 
background: to eliminate it, one works beyond the end-point region of the lepton momentum spectrum 
for 6 — > c processes. This is a very small fraction of the phase space, where theoretical inclusive models 
have relevant uncertainties. 

Predictions for the exclusive channels B — > X^lvi are also model dependent, and here HQET is much 
less useful than in the B ^ D process, because of the presence of a light meson in the final state. We 
shall show in the following sections how to relate B — > nlvi to D ttIvi, and B plvi to Z? — > plvi, in 
the b and c infinite mass limit. The main problem of this approach are the I/ttIc corrections, potentially 
relevant and not under control. 

The effective lagrangian approach can shed light on these semileptonic decays and can give indications 
on the values of the relevant form factors at the zero recoil point. In order to extract information from 
the experimental data the complete dependence of the form factors is required, which goes beyond the 
chiral lagrangian approach. For this reason external inputs, either phenomenological or purely theoretical, 
are required, and, in the next section, we shall discuss this issue in some details. 

5.1 Form factors 

We introduce now form factors that parameterize the hadronic matrix elements of the weak currents. 

In the case of semileptonic decays, P P'lvi {P, P' pseudoscalar mesons) there is no contribution 
from the axial-vector part of the current and the matrix element can be written as 

< P'ip'W'^lPip) >= [{p + pT + ^' , ''^ glFi(g^) - f^iPL^q^Foiq') (201) 

where = q'^f^Q and q = p — p' . There is no singular behaviour at = because i^i(O) = ^o(O). 

The form factor Fi(q'^) can be associated, in a dispersion relation approach, to intermediate states 
with quantum numbers = an d Fo(g^) to states with = 0+. In the limit of massless lepton. 



the terms proportional to in (|20l| ) do not contribute to the rate, so that only the form factor Fi{q^) 
is relevant. 
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For the pseudoscalar to vector matrix elements also the axial-vector current contributes and four form 
factors are required: 



<V{e,p')\ ( - A^)\P{p) >. 
+ i{Mp + Mv) 



Mp + 
e* • q 



Mv 



e-q 



+ it* ■ q 



{Mp + Mv) 
2Mv 

r,1 



q^Aoiq^ 



Al{q') 



M2 



A2{q') 



(202) 



where 



Ao(0) 



Mv - Mp 
2M^ 



^2(0) 



Mv + Mp 
2M^ 



A^{0) 



(203) 



Neglecting the lepton mass, only the form factors V{q'^), Ai{q'^) and ^2(5^) contribute to the decay 
rate. The form factors Ai and A2 can be associated to = 1+ intermediate states, and V to = 1~ 
states. 

The form factor dependence on q^ is still an open question, and, at the present time, there is no 
general theoretical agreement. Quark model calculations are based on meson wave functions, generally 
derived by some wave equation, and make use of them to compute hadronic matrix elements. These 
calculations are normally reliable only at some specific value of q^, and the dependence on the variable 
q^ has to be assumed as an additional hypothesis. The physical region for semileptonic decays covers the 
range < q^ < q^ax = i^'^P ~ ™7r)^ (in the hmit of massless leptons). Close to q'^^^x^ the form factors 
should be dominated by the nearest t-channel pole, located at the mass of the lightest heavy meson 
exchanged in that channel. 

With decreasing q^, the influence of the pole becomes weaker: in a dispersion relation the form factor 
can be written as a P* pole contribution {P* nearest pole for that channel) plus a continuum contribution, 
that, in the narrow width approximation, reduces to a sum over higher resonances P„. We can therefore 
write, e.g., for the form factor Fi: 



Fi{q') 



fP'SPP'Tl 

g2 - Ml. 



^ fp,i9PP, 



(204) 



where gpp*^ is the trilinear coupling among P, P* and tt. In the combined limit — > 0, Mp —f 00, 
the P* pole contribution goes like M^/"^ when q^ — > Mp {gpp'-n ^ Mp), while the higher resonances 

— 1/2 

contributions go only like Mp , as they do not become degenerate with the P in the heavy mass limit. 
But far away from the kinematical end point q^^^^, many resonances can in general contribute to the 
form factors. This observation leads to two-component models for the form factors ||l2|. In [104| it is 
shown that the form factors for i? — > tt semileptonic decay are dominated by the B* pole at all q^ in 
two-dimensional planar QCD; it remains to be seen if such dominance hold also in four dimensions. 
The nearest pole-dominance on the whole q^ range, i.e. 



F{q^) 



(205) 



should therefore be taken as an additional assumption in building models, as done for instance in the 



popular BS W model |105|. Other dependences can be found in the literature; for example in the ISGW 
model [106 1, which is expected to work well close to (Zmax; the extrapolation to lower q^ is done by an 
exponential dependence for the form factors. It is important to stress that in general predictions for the 
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widths are not sensitive to the assumed dependence only when the available range in is not large, as 
in D decays: for B decays into light mesons, as i? ^ T:{p)li'i, different form factor behaviours can lead 
to different predictions. 

Experimentally, only the decay — > K'e^Ve allows at the moment a study of the dependence of 
the form factor Fi. The data are compatible with the pole form (205), but the precision is still poor. The 
value of the pole mass, as fitted by the data, is compatible with the Z?* mass, and the intercept Fi{0) is 



F/'^ (0) = 0.75 ± 0.03 



(206) 



The Cabibbo-suppressed decay D nlvi suffers from poor statistics: MARK III and CLEO II data, 
extracted with the pole-dominance hypothesis, giveES: 



i^]"^(0) 



= l.Otol ± 0.04 MARK III 



1.3 ±0.2 ±0.1 CLEO II 



(207) 



Theoretically, QCD sum rules allow to compute the q^ dep e ndence o f the form factors, except when 
close to q^nax- The analyses performed in |l07| , 108 , 109 , IIC , 111 , 112 are generally compatible with 



nearest pole dominance for the vector current form factor, Fi a nd V : for the axial form factors, Ai and 
A2, there are discrepancies among the different calculations. In |113, 108 1, the form factor Af^'^ has an 
unexpected behaviour, decreasing from q^ — Q to q^ — 15 GeV^: for A^^'', a moderate increase in q^ is 
found, at least up to q^ ~ 15 GeV^. For higher values of q^ the estimate is unreliable. Such a behaviour is 
the resu lt of cancellations among large terms, and therefore can suffer from relevant uncertainties. Also, 
in [112| the form factor Ai decreases with q^, while A2 can be fitted by a pole formula. Light-cone sum 
rules I lie I show, on the contrary, an increasing Ai, with a dependence close to the pole behaviour, and 
a steeper increase for V. 

Current lattice QCD simulations cannot study directly the b quark, because its mass is above the 
UV cut-off. Quantities are computed around the charm scale, and then extrapolated up to the b mass 
using the Isgur-Wise scaling relations With this strategy, suggested in [114], one is forced to make 
assumptions on the g^-dependence at the b scale, because the extrapolation pushes the q^ value towards 
Qmax- Fo'" meson form factors, the determination of q^ dependence is still poor, but compatible with 



pole dominance |115, 116^ Preliminary lattice computations of the q dependence of the form factors 



Fi and Fq in B nlvi |117], and Ai in B plvi |11S], seem to favor a dipole/pole fit for Fi and Fq 
respectively, and a pole behaviour for Ai. the data have however large uncertainties, and will be improved 
by working with heavier quark masses and by using larger lattices. 

An additional constraint to the form factor q^ dependence is provided by nonleptonic heavy meson 
decays. It has been shown that the commonly used form factors, when used together with the additional 
hypothesis of factorization to evaluate non-leptonic decay amplitudes, do not agree with the data on _B — > 
J/ipK{K*) transitions |11£, 120|. These non-leptonic decays can be computed, using the factorization 
approximation, as functions of the leptonic decay constant /j/^ and the form factors F-^^ , Af^ , 

^ , and V^^ , at q^ = Mj^^. The problem is to fit simultaneously the rather small ratio of vector 
to pseudoscalar decay rates, T{B — > K*J/ip)/r{B KJ/t/j), and the large fraction of longitudinal 
polarization in _B — > J/^K*. 

In 1 121 1 , it has been shown that the discrepancy can be eliminated allowing for a non-polar behaviour 
of some form factors. Using the Isgur-Wise scaling laws, these authors compute the relevant form factors 
from the experimental data on the semileptonic transitions D K{K*) at q^ — 0. Subsequently 
they adopt for the form factors a generic dependence (1 — g^/A^)^", with n — —1,0,1,2. Simple-pole 
dominance corresponds to n = 1. Three scenarios survive to the phenomenological analysis: the form 
factors Fi, Ai, A2 and V can only have a dependence [np, ni, n2, ny] = [+1, —1, n2, +2] respectively, and 
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712 can be equal to 2, 1, (with a preference for n2 = 2). Notice that Ai{q ^) is hnearly decreasing in these 
scenarios: a similar behaviour is found also in a theoretical analysis |108|, at least for < 15 GeV^ . 

The hypothesis that Ai is linearly decreasing on the whole range in is certainly not valid: the 
form factor should h ave a pole at q^ = , which should affect the q^ dependence at least close to the 
zero-recoil point. In [120| it is argued that Ai should have a flatter q^ dependence than the one predicted 
by the pole dominance. 

In section 5.3 we shall present a phenomenological analysis of the B ^ V {V light vector meson) 
processes, showing that a constant Ai behaviour leads to discrepancies with the available data, when 
scaling laws are used to scale the form factors from D to B systems: the situation would be even 
worse for a decreasing Ai. There is thus some rough suggestions for an increasing Ai with q^, and in 

a constant term plus a pole term, is used with satisfactory 



5.3.2 a two-component model for Ai 



section 

phenomenological agreement. 

In any event, further theoretical and experimental studies are needed to clarify the problem of q^ 
dependence in the form factors. 



5.2 B ^ TT semileptonic decays 



In this section we shall analyze the semileptonic exclusive decays of a heavy meson P = B, D into a light 
particle belonging to the pseudoscalar octet, e.g. tt. 
The relevant hadronic matrix element is: 



(208) 



In absence of a knowledge of the hadronic current in terms of hadrons, various theoretical approaches 
to the evaluation of (208) have been developed: potential models, for instance in [105; 106|, lattice 
QCD lll|, |rg, |ig, |ll^ QCD sum rules |122[|, and the method based on the chiral and heavy quark 



symmetries, that we shall review here. 

In general two attitudes are possible. In the first one, that we shall call the scaling approach, one 
relates the different hadronic matrix elements (208) using the spin and flavour symmetries. In this way for 
instance it is possible to relate the matrix elements < 7r|j^|i? > and < K\j^\D >, tt and K belonging to 
the same chiral multiplet, and B and D being related by heavy flavour symmetry. In the second approach, 
one builds up an effective lagrangian incorporating chiral and heavy quark symmetries, and computes the 
form factors within such framework. The advantage of the second approach is the possibility to include 
in a rigorous way symmetry breaking corrections, at least formally. Both approaches lead to the same 
results at the leading order, and require some experimental or external inputs: in the scaling approach one 
starts with a known matrix element, while the effective lagrangian contains unknown couplings that are 
determined by the data. As we shall see, some of these couplings have also been estimated theoretically, 
and we shall review also these analyses. 



5.2.1 The scaling approach 

We shall discuss first the so-called Isgur-Wise scaling laws for the form factors. Let us parameterize the 
hadronic matrix element as follows: 



cPtt f„2\ 



(209) 



where q^ = {pp —p^r)'^. The Isgur-Wise relations, following from the SU{2) flavour symmetry between b 
and c quarks, give |5l|]: 



(./+ + /-) 
(/+-/-) 



+ 1/2 



(210) 
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We have neglected here the logarithms of mq arising from perturbative QCD corrections. These scaling 
laws are valid as long as v ■ p-^ does not scale with mq, i.e. in the kinematical regime of soft tt, close to 
Imax = {iTT-P — ™ir)^- In this region E-^ — rriTr << mq. q^-,^^^ is the no recoil point, where the final scalar 
and the dilepton system are at rest. 

Using the form factors Fi.Fq, the eqs. (|21C|) become, including QCD corrections. 



max.B / 



max . B ) 



-6/25 



-6/25 



1 vimax ) 



Ml 



max.D ) 



(211) 
(212) 



Eqs. ( pill) and (212) are valid in the mb, rric —> oo limit. The l/nic corrections can be large, as we 
will discuss later, and can be estimated in the effective lagrangian approach. 

The application of the chiral symmetry is straightforward, and, at the leading order gives: 



F^^(g2)=F^- (g2 



(213) 



Here n and tt' are two arbitrary light pseudoscalar mesons, for instance tt and K. We notice that (213) 
is valid for any value of g^. 

From the knowledge of the form factors at any value of q^ for a given decay mode one can compute 
a whole class of decays as follows: 



i) Using (213), one computes all the chiral-related decays. We notice that all the form factors related 



by light flavour symmetry have the same g^-behaviour, but Fq and Fi can have a different behaviour. 



ii) The Isgur-Wise scaling laws (211) and (213) allow to relate B and D form factors, at least close to 



Hmax 

iii) The strongest assumption concerns the evolution in of the scaled form factors. The B decay 
rates are quite sensitive to the explicit dependence, while for D decays the g^-range is much 
smaller. 

We proceed to the computation of the form factors, widths and branching ratios for the semileptonic 
decays of a heavy meson into a light scalar, following the str ategy we have described. Such an approach 
was followed in [123| for the semileptonic decays, and in |11£, |120| to compute the B K{K*) form 
factors. 

For the form factor Fi we shall assume a simple pole behaviour, which, as we have already discussed, 
agrees with present experimental and theoretical evidence. As an input, we use the decay D — > Klvi, 
i.e. the form factor at = (|206D . For chirally-related decays, we impose the same value to the form 
factors at = 0. For generic q^ we have 



Fi(0) 



l-g2/Af' 



(214) 



p* 



The 1~ pole P* can be a strange or a non-strange heavy meson, depending on the decay mode: their 
mass difference — Mp^ — Mp ~ 100 MeV is a chiral breaking effect, which we will not neglect. 

Let us make some comments on the Isgur-Wise scaling law, which follows from the observation that 
the matrix element < 'k{p-^)\J^\P{v) > behaves as ^mq^ in the limit v ■ << raq 00. These 
asymptotic scaling laws have 1/mq corrections 



< Tr{p^)\J^\P{v) >^ [1 + 



A 



QCD 

mq 



O 



V - Pit 

mq 



O 



myr 

mq 



(215) 
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and a number of different choices for the scahng relations are possible. They are all compatible with the 
asymptotic behaviour, but differ for corrections of the type indicated in ( ^15[ ). Eqs. (211) and (212) is a 
possible one, but other choices might be done. For instance we could have taken, for Fi, neglecting QCD 
corrections: 



pBrr 



Md Mb + 
Mb Mo+m^ 



(216) 



which reduces to ( pil| ) in the M^, Mjy —i- oo limit. Notice that the "soft scaling" ( |21q ) can lead to results 
numerically different when the mass of the light final meson is not so small, as for instance in the case of 
K or K* . In these cases the differences between (211) and (216) are of the order itik/Md, or rriK* /Md- 
Other possible scaling forms can be found in |119, 120|. QCD radiative corrections, written in (211 1212| ), 
are of the order of 10%: and therefore they can be neglected within our approximations. 

The chiral limit, 0, presents some subtleties. In order to examine them, let us find the scaling 

laws at q"^ 
finds: 



0, as arising from the simple-pole behaviour (214) and the asymptotic scaling at q^ax- One 



Ft" 



Md As Ff^ 



Mb Ad + Fj"'^ 



In the mt,, rric oo limit, Ap = Afp. — Mp — s- 0, and one gets: 



Bit 



Md 
Mb 



On the other hand, performing first the m^r 

pBn 



limit in ( piTl) , one finds: 



F- 



(MdY 

[mbJ 



(217) 



(218) 



(219) 



The contradiction between (218) and (219) means a breaking of the naive scaling laws at q'^^ix the 



chiral limit: as shown by pi] , in the limit 

pB. 



(211) becomes: 



F- 



( Mb 
\Md 



(220) 



This behaviour can b e exp lained , as we shall see, from the polar diagram with exch ange of the P*. 
Therefore, combining ( 217 ) and (|22C| ), we find also in the chiral limit the scaling law ( 218 ), which we 
shall use i n the subsequent analysis. 
Using (206) as an input, we get: 



Ff^(O) ~ 0.45 



(221) 



As previously discussed, the uncertainties due to scale corrections are expected to be of the order itik/Md, 
i.e. about 30 %. The uncertainties due to deviations from the polar behaviour of the form factor are 
hard to estimate and essentially unknown. 

We could use as an input the Cabibbo-suppressed decay D tt, which has however larger experimental 
errors. The scaling uncertainties, moreover, even if probably smaller than in the case of a final K , are 
expected to be of order A/Md (A = Mp — mq) and not of order rriTj/MD only; therefore they could be 
really significant. 



The prediction for widths and branching ratios following from (221) and chiral symmetry are reported 
m table |. 
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Table 3: Predictions for semileptonic D and B decays in a pseudoscalar meson, in the scaling approach. 
We have neglected the 77 — 77' mixing. The branching ratios and the widths for B must be multiplied for 
|T4;,/0.0032p. We assume tb, — tbo = r^-f = 1.55 ps. 



Decay 


Fi(0) 


BR 




exp. BR 




0.75 


3.4-10- 


-3 


(3.9;^ j) - 10-^ 


D+ ^T] 


0.31 


6.9-10- 


-4 






0.61 


3.2 - 10- 


-2 






0.75 


4.2 - 10- 


-3 




B^-^Ti- 


0.45 


2.2 - 10- 


-4 


(1.63 ±0.46 ±0.34) - 


Bs^K 


0.45 


2.2 - 10- 


-4 
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The available experimental data shown in table || have large uncertainties. Concerning the B^ 
■K^l^vi decay, the CLEO II collaboration |l03| quotes for the branching ratio two different values, de- 
pending on the model used for the detector efficiency: in the previous table ^ we have put the number 
corresponding to the BSW model |105| 



BR{B° 



'l+vi) = (1.63 ± 0.46 ± 0.34) - 10" 



while the value corresponding to the ISGW model |106] is 
BR{B^ ^ TT^l+n) = (1-34 ± 0.35 ± 0.28) - lO"'' . 



(222) 



(223) 



5.2.2 Effective lagrangian approach 

We now discuss the effective chiral lagrangian approach to the semileptonic heavy-light form factors. 

We have already presented the effective lagrangian that combines heavy quark and chiral symmetry 
and describes the low-momentum interactions of heavy mesons with light pseudoscalars, and the chiral 
representation in terms of meson fields of the weak current 97^(1 — 15)Q- In this framework one can 
compute the hadronic matrix elements < 7r| J^|P(?;) > in terms of the effective couplings of the lagrangian 
and of the weak current, at least in the soft-pion region, i.e. close to q^ax- Two diagrams contribute to 
the form factors, at least in the leading order: the P* pole diagram, proportional to the strong coupling 
constant g, and a direct diagram, as shown in fig. ||. At the leading order in l/mg they give |jl^, |l^ : 



B 



B 



B^ 



Figure 5: Tree diagrams for the B n transition matrix element. The box represents the h —> u weak 
current, the dashed line the pion. 



F 



(224) 
(225) 
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The same relations have been obtained assuming P* po le dominance in , and in 1 124 



PC AC and heavy quark spin symmetry. The formula (224) s atisfi es the asymptotic scaling (211): in the 
chiral limit v ■ ^ 0, however, the scaling is modified as in (22C), because Ap ^ 1/Mp. The scaling at 



q = is easily derived from (217) 



combining 



pDTT 



Mp 
Mb 



(226) 



as in (218). Therefore, if we use the input D - 
the same results as in table |[ as expected. 
Alternatively, one can take the estimates g ; 




■ to fix the effective coupling gF in (224), we obtain 
0.38 (|9|) and F = 0.30 ± 0.05 GeV^/^ (|^, which give 

(227) 
(228) 



The uncertainties on these results, as arising from QCD sum rule approximations, can be estimated around 
30%. We stress that at the charm scale the 1/toc corrections are potentially relevant: nevertheless the 
prediction (227) is in rather good agreement with data. In table ^we quote form factors and branching 
ratios obtained from (224), (|225|) and chiral symmetry. 



Table 4: Predictions for semileptonic D and B decays in a pseudoscalar meson, in the effective lagrangian 
approach. We have neglected the rj — rj' mixing. The branching ratios and the widths for B must be 
multiphcd for |Kib/0.0032p. We assume tb^ — tb« = tb+ = 1-55 ps. 
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Fi{0) 


BR 




Exp. data 


£)0 


0.63 


2.4- 10- 


-3 


(3.9+^;^) • 10-^ 


D+ ^t] 


0.26 


4.9 • 10- 


-4 




Ds^V 


0.51 


2.3 • 10- 


-2 






0.63 


3.0-10- 


-3 




B'^^TT- 


0.38 


1.6 • 10- 


-4 


(1.63 ±0.46 ±0.34) • 


Bs^K 


0.38 


1.6-10- 


-4 
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The effective lagrangian result (E2 



shows the dominance of the P* pole near the kinematic end- 
28|) have been presented in [p5|, where it is 



point. The l/mq corrections to the leading results (227), 
shown that at the order 1 /mg one has to introduce four new relevant couplings in the lagrangian and in 
the effective current. Two of the occurring parameters, pi and p2, parameterize the 1/M corrections to 
the leptonic decay constants fp and fp*, as shown by (197). The others, gi and 172, are related to the 
subleading corrections to the coupling gp.p,^^ and gp,p^^ as shown in (pT]). In terms of these couplings, 
the form factor Fi is |3q|: 



Fi = 



F 



Mp + vp^ 
V ■ Ptt + A 



Mp + v p^ 
yMp{vp^ + A) 



1 



Pi - 2p2 
Alp 



gp'PTrfp* — fp 



.91 - 92 

Mp 



- 1 



Pi + 2p2 
Mp 



(229) 



where, in the second expression, gp^p^,, fp, and fp include their own 1/M corrections. 

The chiral logarithmic corrections to the process P — > ttIvi in the effective theory have been computed 
in j6^: evaluating the one-loop chiral diagrams, it is found that t he ch iral-corrected form factors Fi and 
Fq at the leading order in l/ruQ have the same form as in (224), (225) but are expressed in terms of the 
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chirally renormalized leptonic decay constants /JJf" and f^'^", the heavy meson coupling g^ff to the axial 



vector Goldstone current (see also section 3.2) and the hyperfine mass splitting A'"'^", i.e 



^ gefffF"Mp 



ren 



Fo = ^ . (231) 

J TT 

A more detailed analysis is presented in , where all the non-analytic terms arising from chiral loops 
are kept. The 5'[/(3) violation in the pole term of the amplitude is of the order of 40%, but a rather large 
value of g is used in the numerical estimate, g ~ 0.7. For smaller values of e.g. g — 0.3, the chiral 
violation between P tt and P ^ K pole amplitudes reduces to 10%. 

5.3 B ^ V semileptonic decays 

We now discuss the semileptonic decays of a heavy meson P into a light vector meson V ~ p, K*, (f>. 

In the following, we will discuss how to relate the B and D form factors, following as before two 
different approaches: the scaling approach and the effective chiral lagrangian approach. Let us begin 
with a review of the available experimental data. In the D —^ K*lvi decay, the most extensively studied 
channel, the quality of the data does not yet allow to determine the dependence of the form factors. 
The analysis is performed assuming a simple pole formula for the form factors V{q^), Ai(g^), and ^2(9^), 
with pole masses given by the nearest resonance (i.e. 2.1 GeV for the vector form factor and 2.5 GeV for 
the two axial form factors). The average of three Fermilab experiments gives the results ]45[| : 

V{Q) = 1.1 ±0.2; ^i(O) = 0.56 ±0.04; ^2(0) = 0.40 ± 0.08. (232) 

Data have also been obtained for the decay Ds — + 4>lvu but the errors on the form factors are still 
large, and we shall not use them. 

In the case of semileptonic D decays, due to the limited range, the pole assumption does not 



sensibly affect the results (232); we have for instance extracted the form factors assuming Ai constant in 
q^ , but A2 and V pole-dominated, finding discrepancies of the order of 10%, which are within the quoted 
uncertainties in ( ^321 ). The q^ dependence of the form factors is on the contrary extremely important in 
B decays, as we discuss below. 

For B mesons the semileptonic rates are strongly Cabibbo-suppressed: the CLEO II Collaboration 



has only recently presented the new measurement giving |103| 
BR{B° ^ p-£+iyi) = (3.88 ±0.54 ±1.01)- 10"* WSB 

BR{B° ^ p-£+iyi) = (2.28 ± 0.36 ± 0.59) • 10"* ISGW (233) 



where the first value is obtained using the WSB model ]105| | in the Montecarlo code which evaluates the 
efficiencies, and the second one is based on the use of the ISGW model |l06| . 

Before discussing semileptonic B decays in more detail, let us stress that another source of information 
on the weak matrix elements between B and K or K* is represented by non leptonic B decays. As a 
matter of fact, the factorization hypothesis allows to relate non-leptonic to semileptonic rates. The 
color-suppressed decays B — > K{K*)J /tp give, in this approximation, indications on the form factors 
B K{K*) at q^ = Mj^^. There are two relevant experimental figures: the ratio of vector and 



pseudoscalar widths, measured by Argus |126| and CLEO II [125 , whose averaged value is 



T(B ^ J/t/jR'*) , , 
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and the fraction of longitudinal polarization 

Tl ^ T{B ^ J/^K*)l 
r T{B J/iPK*) 



0.74 ±0.07 



(235) 



Detailed phenomenological analyses have been performed in |119, 120, 121|, where it has been shown 
that most of the models fail to explain the previous data, in particular the fraction of longitudinal 



which corresponds to the average of the measurements of Argus: Fl/F — 0.97± 0.16± 0.15 |126|, CLEO 
IF Fl/F = 0.80 ± 0.08 ± 0.05 and CDF: Fl/F = 0.65 ± 0.10 ± 0.04 ^2?" 



polarization (235). We point out that all the current models use the hypothesis of factorization |128] 
which in general works satisfactorily in B decays |129| but could have corrections in specific channels 
like the colour-suppressed B — > J/ipK*. Possible non-factorizable contributions are introduced in il3i 
By the definition 



V 



BK* 



(M 



j/i,' 



one gets, assuming factorization [119 



R = 1.081 



Fl 
F 



aBK' \- 

^(M'j/i,)) [(«-M' + 2(l + cy2)] 
(a - bxf 



(236) 



(237) 



(238) 



(a-6a;)2 +2(l + cy2) ' 
The coefficients a, h and c are dimensionless combinations of masses; from the data one gets 

a = 3.16 6=1.31 c = 0.19 . (239) 



5.3.1 Scaling approach to B ^ V form factors 

The scaling approach, valid at the leading order in l/mg, is similar to the case P — > tt; the scaling laws 
for the form factors are derived from the asymptotic behaviour of the matrix element, i.e. 



<V{p')\J^\P{v) >^^[l + 



A 



QCD 



O 



V ■ p 



O 



(240) 



For D K* the violation to ( p40D can be important, namely of order m^* /Me,. This uncertainty is 
reflected in different choices of the scaling laws at ~ q^nax ■ ^'^'^ example one can follow the approach 



Mq. 
Mql 



' ) 

max / 



^ ) 

max I 



Mp 


+ Mv 




Ifp 




Ifp 


Mp 


+ Mv 


Mp 


+ Mv 



/Mp 



(241) 



The second choice we shall consider is "hard-scaling" : 
V{ql 
Aliql 



max / 



"imax / 



3 max ) 




(242) 
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Table 5: Predictions for form factors and widths for B'^ p~lvi- and Tl refer to p with respectively 
transverse and longitudinal polarization, r+ and r_ to p with positive and negative helicities. The 
branching ratios (BR) and the widths for B must be multiplied for |Ki6/0.0032p. We assume tb, = 
tbo — tb+ 1.55 ps. 



Extrapolation 


soft-pole 


hard-pole 


soft-cost. 


hard-cost. 


data 


Ai(0) 


0.21 


0.17 


0.42 


0.33 




^2(0) 


0.27 


0.34 


0.27 


0.34 




y(o) 


0.64 


0.81 


0.64 


0.81 




r+/r_ 


0.045 


0.13 


0.005 


0.06 






0.39 


0.15 


1.80 


0.30 




BR (10-") 


2.8 


2.8 


6.7 


3.5 


3.88 ± 0.54 ± 1.01 


TL/r{j/^K*) 


0.27 


0.02 


0.50 


0.18 


0.74 ±0.07 


r{K*)/r{K) 


1.74 


1.68 


3.35 


2.27 


1.68 ±0.33 



The two scenarios, (^4l|) and (242), differ by subleading terms of the order My/A/p, which can be 
nevertheless numerically important. Needless to say, some form factors might exhibit soft-scaling and 
others hard-scaling, in different combinations. 

The scaling laws allow to relate the D and B form factors near q^nax- as discussed above, the de- 
pendence on is practically unknown, and, as we stressed already, B transitions depend strongly on 
the extrapolation to — 0. The vector form factor V is generally believed to be pole-dominated, as 
discussed in section 5.1, while for Ai and A2 the theoretical situation is unclear. In [120| the soft scaling 
laws are justified by extending the heavy-to-heavy scaling relations down to the light final meson case. 
In the same limit one finds that A2/A1, V/Ai and Fi/Ai should have a polar behaviour in q^: assuming 
Fi as pole dominated, this implies a constant Ai and a pole behaviour for A2 and V. Nevertheless the 
extension of the heavy-to-heavy scaling laws to the heavy-to-light case remains arbitrary, and should 
be considered as an ansatz. As we will di scuss explicitly in the next section, the effective lagrangian 
approach leads to the soft-scaling solution (241): this follows from the factor {Mp + My) contained in 
the definition ( 202 ) of the form factors. 

To simplify the discussion, we assume that A2 is dominated by the nearest pole, while for Ai we 
consider two possibilities: the pole-dominance, and a flat Ai constant in q^. We have considered four 
different possible scenarios: soft scaling and Ai pole-dependent (called soft-pole), soft scaling and Ai 
constant (called soft-constant), hard scaling and Ai pole-dependent (called hard-pole), and finally hard 
scaling and Ai constant (ca lled hard-constant). For each of them we have computed, using as inputs the 
D K* form factors ( |232| ), the branching fraction and the ratios of decay widths r_|_/r_ and Tl/Tt 
for the process B^ —>■ p^l^vi. Here Tt and Tl refer to p with transverse and longitudinal polarization 
respectively, r+ and r_ to p with positive and negative helicities. We have computed the longitudinal 
fraction Fi/F, and the ratio of the vector to scalar BR's for B J/ipK{K*). The results are presented 
in table where we have extrapolated from D — > K* to B ^ K*, and then we have equated, by chiral 
symmetry, the form factors B K* and B ^ p. 

Let us comment on table |^. First of all, all the four scenarios give a rather low value for the ratio 
Fl /F ( J/ ipK*) : the soft-scaling, Ai constant (third column of table ||) , is the closest one to the experiment , 
but it produces a too high value for the branching ratio B p^lvi and for T{K*)/T{K). The value of the 
latter ratio depends mainly on the ratio Ai/Fi (at q^ — Mjy^), and could be smaller for a larger value 
of Fl ^ : however, a too large value would disagree with the measured branching ratio for the B — > ttIvi 
(excluding large SU{3) violation). This scenario is preferred in |12C|, where however the upper limit for 
the B ^ p and the _B — > tt data are not taken into account (soft-scaling is also applied to Fi, obtaining 
a larger value for it and a better agreement for r{K*)/r{K)). 
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The second comment is that a constant Ai gives a higher value than pole behaviour for the B ^ p 
branching ratio, when scaling is used. The situation improves assuming pole dominance for Ai for the 
semileptonic branching ratio B ^ p, and also for the ratio T{K*) /r{K), as it can be seen in the first 
column of the table. But, at the same time, the longitudinally polarized fraction rL/T{J/ip K*) decreases, 
because A2/A1 at = Mj^^ grows. 

Hard scaling decreases the value of Ai, and it raises A2 and V: this produces a smaller semileptonic 
BR{B p), but also lowers TL/T{J/ij K*) (because A2/A1 grows, see (p3^)) and T(K*)/T{K). Hard 
scaling, together with a pole-dominated Ai, as in the second column, leads to a value for V l/T{J /"ip K*) 
in disagreement with the data: its combination with a constant v4i, as in the fourth column of the table, 
improves the agreement with the data, even if TL/T{J/tp K*) remains still rather small, even smaller 
than in soft-pole scenario for Ai. 

Summarizing, the previous analysis indicates that a constant Ai requires strong scaling in order to get 
agreement with the data of the semileptonic BR{B ^ p). In ]120| | a different result is obtained, which, 
however, does not take into account the results for the semileptonic _B — > p transition. 

If Ai is single-pole dominated, soft scaling is required in order to get a reasonable (not too small) value 
for Tl/T{ J/^ K*): these data are however difficult to explain without spoiling other phenomenological 
requests. 

It should be stressed, however, that the figures of table ^ have large uncertainties. Leaving aside 
theoretical uncertainties that are however significant, the q uote d numbers have an uncertainty due to 
the experimental errors of the D K*lvi form factors (|232| ) used as inputs. We notice that, for 
instance, ^2(0)^ ^ is quoted with an error of about 20 %: this error alone implies 30 % uncertainty in 
TlIT{J/4! K*) and 15 % in T{K*)IT{K). Finally we have used flavour SU{?>) symmetry to relate B ^ p 
and B K* form factors and this is another source of theoretical error which in principle should be 
taken into account. 



5.3.2 Effective lagrangian approach 



Light vector resonances have been discussed in the framework of the effective heavy meson chiral la- 
2.4; applications to semileptonic decays have been developed in Chiral loop contributions 

. We now review this subject. 



grangian m 

to D ^ K*lvi have been partially taken into account in 



In the effective chiral lagrangian framework, five different diagrams contribute at the leading order in 
l/niQ to the matrix element < V{p' , e)\Jn\P{p) >. They are analogous to the diagrams of fig. ^ Four 
of them are polar diagrams: V{q'^) takes contribution from the 1~ pole diagram, proportional to the 
coupling A among P, P* and p introduced in (|2|). The C and p couplings in (^ give the vertex PP**p, 
where P** is a positive-parity states of the doublet (0+, 1+) and the corresponding polar diagrams, with 
exchange of a 1"*" meson, contribute to Ai and A2. The p and A terms have dimension higher than the C 
or g terms; nevertheless they are the lowest-order contributions to V and A2. The fourth polar diagram, 
with exchange of a 0~ meson, is proportional to the coupling /3 of the vertex PPp (see formula (|6^) and 
is relevant only for the form factor Aq. Finally there is a direct diagram, which arises from the effective 
current term proportional to ai in the heavy-to- light current ( 184 ). All these vertices can be found in 
appendix A. 

Imax ^1^*^ leading order in l/mq, one gets [p^ : 
V{qLJ = ^AF — \ (243) 



Computing the diagrams for q^ 
9v .f^Mp + Mv 

— —;=At 

V2 



' ) 

max / 



2gv 




V2 Mf 



M{q, 



) 

max/ 



F+ 



Mv 



ai 



F+{C/2-pMv) 
Mv + App,* 



Mp + Mv 



V2{Mv + App,,) 



(244) 



(245) 
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Mil 



9v 



9v an 



2\/2Mv(My + A') V2 Mv 



(246) 



where the A are appropriate mass spUttings. 

From (245)- ( [24^ ), we can extract information on the scaHng behaviour of the form factors as well as 
some indications on their behaviour. 

245|) imply the "soft-scaling" of eq. (^41|), since the coupling 



As to the scaling, we notice that (243)- 
of the effective lagrangian is flavour independent. Moreover, the form factors V and A2 have only the 
pole structure, which at ~ (Zmax signalled by the factor l/{Adv + A), while Ai contains a pole term 
but also a non-polar one, proportional to ai. This suggests a more complex q^ behaviour of Ai. 

The non-polar term could be a general polynomial in q^: the simplest way to take into account the 
indications coming from the effective lagrangian is to describe Ai as a sum of a constant term and a pole 
term, i.e. we write: 

while keeping for A2 and V the nearest-pole behaviour. The parameters a and b are flavour dependent, 
and scale differently at q^ = qf^ 
following scaling laws: 



Assuming soft scaling, as suggested by (p43f) - (246), we find the 



a{B 



V) 
V) 



Mb + Mv\I Md"'^ 



Md+M^ 
Mb + Mv 



V) 
V) 



(248) 
(249) 



The formulas (243)-(245) give the form factors at g^ax- the leading order in l/mg, the value in 



for pole-dominated terms is: 



F(0) 



2 (My + A) 
Mp 



Fiq 



' ) 

max/ 



(250) 



As for Ai, we identify the term proportional to ai in (244) with the constant term a of (247), and the 
other with the pole term, proportional to b. In this way one gets: 



ViO) = gv 



V2XF 



Mp + Mv 



M 



3/2 



(251) 



a = —qvV2 



-ai 



9v 



2V2F 



Mp + Mv 

C/2 - fxMv 



{Mp + Mv)VM^ 



^2(0) = -V2^igvF- 



_Mp + Mv 



M' 



3/2 



(252) 
(253) 
(254) 



From the experimental data ( p32[ ) on D — > K* we can fix the effective couplings appearing in the 
previous formulas: from (251), using F = 0.30 GeV^^^ and gv — 5.8 (see (p7[)), we obtain 



|A| = 0.41 GeV 



(255) 



We shall see, in the next section, how the sign of A can be fixed. It is intere sting to observe that this result 
agrees with the second, but not with the first determination obt aine d in 131 by a light cone sum rules 
calculation (the first determination gives a higher value; also in 132 a higher value of A is obtained). 
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Similarly, from (254) and F+ = 0.46 GeF^/z g get 



^i = -0.10GeV-' . (256) 

Having fixed A and /i, one can compute the form factors V and A2 for the B K* and B ^ p matrix 
elements, because these couplings are heavy-flavour independent at the leading order. The result is 

|ys-*^(0)| =0.50 Af^^(0) = 0.19 . (257) 



The mass difference between K* and p is numerically irrelevant in (257). 
Concerning Ai, from the input A^ ^ (0), we derive 

(a + 6)^^* = 0.56 . (258) 

The knowledge of the behaviour of ^ would allow to extract simultaneously a and 5, but at 
present we can only introduce an arbitrary parameter r, defined as 

" (259) 



a + 



where a and h are relative to the A^^ (r is not heavy flavour independent). 

When r varies from r = to r = 1, we have a smooth transition from a pure pole dominance (r = 0) 



to a constant Ai [r — 1). The analysis of section 5.3.1 has shown that a pure pole behaviour (soft-pole 



case) leads to a rather low value for Tl/T{ J/^ K*), while a constant Ai gives a semileptonic branching 



ratio B ^ p too high: therefore we expect that the two component form factor (247) can explain better, 
for some intermediate value of r, the large longitudinal polarization in _B ^ J/^pK* and at the same time 
can agree with the data for B plv. 

In figs. I, and I we plot respectively TL/T{J/i> K*), T{K*)/T{K) and BR{B" p-l+vi) as a 



function of r. We have assumed A2 and V pole-dominated, with the values (257) at q^ — 0. We see 
that all the three observables grow when the q^ dependence of Ai becomes flatter, i.e. r — 1. This is a 
satisfactory feature for Tl/T{J/iI] K*), but a large r gives too large values for the semileptonic BR and 
for the T{K*)/T{K) ratio. 

In table 1^ we quote the values of various observables at r = 0.5, where a good compromise is obtained: 
T l/T{J /ip K*) = 0.51. This value is still smaller than the data, but one should not forget that the 
factorization assumption could receive sizeable corrections. 



Table 6: Predictions for form factors and widths for B^ p^lvi and for B —> J/ipK'^*\ with r — 0.5. 
Tt and refer to p with, respectively, transverse and longitudinal polarization, r+ and r_ to p with 
positive and negative helicities. The branching ratios (BR) and the widths for B must be multiplied for 
|T4f)/0.0032p. We assume tb^ = t^o — tb+ = 1.55 ps. 



Observable 


r = 0.5 


data 


^i(O) 


0.28 




^2(0) 


0.19 




V{Q) 


0.50 




r+/r_ 


0.013 




Tl/Tt 


1.60 




BR (10-4) 


3.8 


3.88 ± 0.54 ± 1.01 


TL/nj/i'K*) 


0.51 


0.74 ±0.07 


T{K*)IT(K) 


2.01 


1.68 ±0.33 



From the value of r one can extract a and h separately: from the scaling relations (248) and ( 249 ) we 
have for r = 0.5: 



a{B ^V)^ 0.21 



b{B ^V)^ 0.07 



(260) 
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Figure 6: Ratio Fl/F for the decay B 
longitudinally polarized K*. 



J/^K* as a function of the parameter r. Tl is the width for 



i .e. Af '^(0) = 0.28, as quoted in table ^. The knowledge of a and b, together with the identifications 
(|252| ) and (|253| ), fixes the couplings ai and the linear combination C/2 — /iMy. The results depend on 
the value of r; for r — 0.5 we get: 



ai = -0.07 GeV^^^ 



c 



fiMv = 0.14 



(a and b are taken as positive). From ( |256| ) we can extract finally 
C = 0.10 . 



(261) 



(262) 



The previous phenomenological analysis has to be taken cautiously, due to the large uncertainties. 
Subleading corrections, dependence of the form factors, breaking of factorization and chiral violations 
could easily lead to substantial modifications of the chosen scenario. New experimental data will hopefully 
clarify the situation, and allow to distinguish among different models. We shall adopt in the following 
the effective lagrangian results of table |6| (in particular r ~ 0.5). 

In table [t] we present the values of the form factors of the b ^ u transitions in different models. 
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Figure 7: Ratio of the widths B J/tpK* and i? J/tpK as a function of the parameter r. 

6 Radiative decays 

6.1 Flavour conserving radiative decays: D* D'j 

In this section we shall consider the decay 

D: ^ i^a7 , (263) 

and the related processes for the B case: B* Ba'j- In ( |263| ) a = 1,2,3 is the light quark index 
corresponding to u, d, s. The matrix element for this radiative transition is as follows: 



In (264) is the photon polarization, whereas the coupling comprises two terms: 



^J■a = ^J.a + ^J■a > 

corresponding to the decomposition: 

M{D: ^ Dal) - e 6*^ < DMJ^-^mp,^) e e*^ < Da{p')\Jf, + J/^li^Ife ^y) > 
Here and are the light and the heavy quark parts of the electromagnetic current: 

3 



1 - 1_ ^ _ 



(264) 
(265) 
(266) 

(267) 



o=l 
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Figure 8: Branching ratio of the deeay B'^ ^ p Ivi as a function of the parameter r. We have taken 
Vub = 0.0032. 



and 



Correspondingly, eq. (265) becomes 



e , h , eg 

A^a = Aia + Ma = T- + , 

where Aq and Aq are mass parameters to be determined. 



(268) 



(269) 



Let us consider the two currents j/j and separately. The matrix element of j/j can be obtained 



from the Lagrange density: 



C" = — 



2mQ 



(270) 



which allows the transition Q —^ Qj and can be expressed in terms of the Isgur-Wise universal form 
factor £^{v ■ v') as follows: 



< Da{p')\j';\D:{p,ij) >= ee < Da{p')\crc\D:{p,Tj) >= i-^ Md^Md:^^^' K^cf^v" v^v'^ , (271) 
where p' = M^v' , p — Md»v and v ■ v' ~\ because: 



Q = q^ =m% + Ml,, - 2MdMd-v ■ v' . 



(272) 
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Table 7: Form factors at 


g2 ^ for 6 - 


^ u transitions 


in different models 




Reference 








yB^p 


This paper 


0.38 


0.28 


0.19 


0.50 


QCD sum rules 



DP 15 

CZ p. 

BBD^Oj 



Narison [|134 
Ball ll08| 
BKRJIO 
ABS m 



0.4 ±0.1 
0.36 
0.24 ±0.025 
0.23 ±0.02 
0.26 ±0.02 
0.24-0.29 



0.35 ±0.16 
0.5 ±0.1 

0.24 ±0.04 



0.42 ±0.12 
0.4 ±0.2 



0.47 ±0.14 
0.6 ±0.2 

0.28 ±0.06 



Quark models 
BSW [|105i 
ISGWjlOel 
FGM psf 



0.33 
0.09 
0.21 ±0.02 



0.28 
0.05 
0.26 ±0.03 



0.28 
0.0.02 
0.30 ±0.03 



0.33 
0.27 
0.29 ±0.03 



Lattice 



APE ]116| 
Abada et al. 

UKQCD Ins 



|114| 



0.35 ±0.08 
0.30 ±0.14 ±0.05 



0.24 ±0.12 
0.22 ±0.05 



0.27 ±0.80 
0.49 ±0.21 ±0.05 



0.53 ±0.31 
0.37 ±0.11 



Taking into account the normalization of ^(w • v'): ^(1) — 1, one gets, for the charm case, 

3Ae 

(^Jj = — 3^ for the b case), with 



(resp. Afj = y/Msa Mb^ ) ; in the leading order in 1/toc one finds: 

Ac = TOc 



(273) 



(274) 



(275) 



(resp. Ab = TOfe), independently of the light quark label a. 

Let us now consider the second term in (265), i. e. /x^, which cannot be computed within HQET 
since it involves light quarks; we shall now show that the chiral effective theory can be employed to get 
information on this quantity. We shall examine two approaches: the first one is based on the calculation of 
chiral loop corrections the second is based on the use of Vector Meson Dominance (VMD), together 
with the effective chiral lagrangian for light and heavy meso ns ||136| . Other approaches used to compute 
(|263|) are based on quark models E^, |2^, 137|; bag model |13S] and QCD sum rules |]139|, 140| (for a 



previous review of theoretical results see |141]). 

The first approach we consider is based on the chiral loop corrections to the tree diagram 
us start with the definition of /x^: 



in the limit of SU{i) symmetry the constants Aa's are equal, i. e. one gets Aq" 
unknown constant which can also contain effects suppressed by powers of l/rric. 



Let 



(276) 



/3, where /3 is an 
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The leading SU{3) violations to (276) are obtained by considering the loop diagrams of fig. 
the results [EtI: 



9 



2 m-K 



M2 



f 1 ? 



|, with 

(277) 
(278) 
(279) 



K 



Here g is the strong coupling constant of the vertex D*Dtt (in the jtIc oo limit) defined in (|40|). If one 
considers only the leading S'[/(3) violations, one should put Jk — fj^- This has to be the case if one uses 
the value g = 0.38, eq. (|99|). As a matter of fact, as discussed in j72|, the SU{2>) invariant coupling g is 
obtained assuming a unique value ~ HQMeV for all the light pseudoscalar meson decay constants in 
the sum rule. The analysis of |M provides a pattern for SU{i) violations in D* radiative decays, and 




y 



D 



D 




Figure 9: Chiral loops contributing to the radiative decays D* 

can be in principle used to determine g and /3, independently of the theoretical determinations based on 
the the QCD sum rule for g. One can use the two formulas: 



^•^^^^ 3M^ 



(k = photon momentum) 



*+ 



D^TT 



(280) 



(281) 



and the experimental results contained in table |l| ^ , together with the condition g < 1 (which is 
experimentally satisfied |Q). Because of the large experimental error (especially in the channel D*+ ~> 
-D~'"7), one gets, however, a rather broad range of values for g [p7[P|: 



0.3 < g < 0.7 



(282) 



The value of g obtained in this way is an effective coupling which takes into account part of the l/rric 
corrections, as it is obvious from the fact that /i^ has not been neglected in comparison with /i^ (/i^ is 
not negligible because rric is not sufficiently large: in this analysis one uses = rric = 1.7 GeV). It is 



^In the analysis of Ref. ^ fx 
^Similar results are obtained in 



.22/,r is used. 
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nevertheless interesting to observe that the small values for g in (282) are in broad agreement with the 



results of the QCD sum rules quoted in section |3.1.2 
The analysis of 



shows that smaller values of g favour values of (3 near the non relativistic quark 
model result 137 , where (3 — m^^, and mq ~ 300 — 500 MeV is a typical value of the light quark 
constituent mass. In particu lar, from g = 0.38 and rric — 1.7 GeV one gets the value f3 = 1.9 GeV^^. 
The pattern displayed by eq.(27£), fil — — — 3f;^mif /(47r/|.)) = — 1/(3A3), can be interpreted, in the 
quark model, as due to a constituent strange quark having a mass ~ A3 larger than rriq = Ai, or A2, 
which is what one would naively expect. 

Instead of considering loop effects one can therefore take into account SU{3) violations by choosing 
explicitly different Aq's. Quark model calculations that assume SU (3) violations arc considered in p9|, |27[|. 
In particular, the calculation of ref. |2^ in the semirelativistic quark model of ref. discussed in section 
as we stated already, makes use of the Salpeter equation |56l, i.e. a wave equation which takes into 
account relativistic kinematics, with an interquark potential modelled on the Richardson's potential | |55[ |. 
As we have observed, this model is able to explain the reduction of the value of the strong D*Dtt coupling 
constant from the (non relativistic) quark model prediction g = 1 down to g w 0.33, as a consequence 
of the relativistic kinematics relevant to the light quark in the heavy meso n; si nce this small value is 
favoured by the QCD sum rule analyses as well as by the data (eq.( ^82| )) one may take this as 

an indication that relativistic kinematics plays a role also in the case of the radiative decays. The results 
of the analysis in l27[| for the constants Aq are display ed in table ^, together with the results of the chiral 
loop calculation [[47[ , i.e. the results based on eqs. (277)- (|279| ). In the same table we also report the 
parameters of the model [136| based on Vector Meson Dominance, to be discussed below. 



Table 8: Theoretical inputs for mass parameters in radiative D decays, x-loop represents the chiral 
loop calculation, VMD is the model based on the effective lagrangian supplemented by the hypothesis of 
Vector Meson Dominance; RQM refers to the relativistic quark model; Ag and Aq are mass parameters 
(in GeV). 

Decay mode x-loop VMD RQM 







Aq 


Aa 


Aq 


Aa 


Aq 


Aa 


D*+ 




1.7 


0.61 


1.9 


0.50 


1.57 


0.48 


D*" - 




1.7 


0.79 


1.9 


0.50 


1.57 


0.48 


D*s - 




1.7 


1.11 


2.0 


0.59 


1.58 


0.497 


B*+ 


B+-f 


5.0 


0.79 


5.3 


0.51 


4.93 


0.59 


B*" - 




5.0 


0.61 


5.3 


0.51 


4.93 


0.59 


b:- 


^Bs-f 


5.0 


1.11 


5.4 


0.60 


4.98 


0.66 



In the case of the chiral loop calculation |47 , we have assumed as an input g = 0.38, which is 
the intermediate value among the different QCD sum rules results [^5[ |6^, Q; on the other hand (3 = 
1.9 GeV~^ is fitted from the experimental CLEO data of table |^, using the branching ratio of D*^ Z?°7 
as an input. As for Ac, following we take Ac = 1.7 GeV; on the other hand for Ah we take the value 
Ah — 5 GeV which, similarly to the Ac case, is slightly larger than the value derived by QCD sum rules. 

Let us now discuss the model based on Vector Meson Dominance (VMD) |136[| . In this model the 
calculation of fi^ is based on the results obtained by the effective chiral lagrangian approach. The idea is 
to use VMD to express < Da\Jf^\Dl > in terms of < DaV\D* > {V — light vector meson resonance) and 
then to employ information from heavy meson weak decays to compute < DaV\D* >. In other terms 
one writes: 



<Daip')\Jl\D:{p,r^)> 
3aY.<Da{p')V{q,e,{\))\D:{p,r,)> 



<0| \V{q,e,{X))> 



(283) 



v,x 
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where = and the sum is over the vector meson resonances V = uj, p^, 4> and over their hehcities. The 



vacuum-to-meson current matrix element appearing in (283) is given, assuming SU{3) flavour symmetry, 
by: 



< 0\qal^qa\V{q,ei) >= ei^fvTr{VT-) 



(284) 



where {T°')im = iai^am and, as usual, a = 1, 2, 3 for u, d, s respectively. From uj —>■ e^e^ and p'^ —^ e+e~ 
decays Q one has fuj ^ fp = fv = 0.17 GeV^; from (jj e+e" one obtains = fv + Sf, with 
5f = 0.08 GeV^, which implies a relevant SU{3) violation. Using (284) and the strong lagrangian 
containing the vertex D*DV (see eq.(p2)): 



C ^ iX < Hba^-" F^,{p)baHa > 



one can compute (283). The results in terms of the mass constants are as follows: 



Air 



-Pa 



(285) 



(286) 



where pi — p2 



fv/ml, p3 



fcji/m^, gv = 5.8. Equation (255) only gives the absolute value of A, 



but eq. (286) clearly shows that A < if has to be interpreted, as in the quark model, as a mass 



parameter. Therefore we take (see eq. (255)): 

-1 



A: 



-0.41 GeV" 



(287) 



The results of this approach are reported in table together with the chiral loop (x-loop) and the 
relativistic quark model predictions. 

From eqs. (280) and (281) and from table || we get the decay rates and branching ratios (BR) for 
both D* and B* decays; they are reported in table |^ for the three models examined so far. For the chiral 
loop calculation and the VMD approach we use the same value g = 0.38 for the strong BB*tt coupling 
constant, whereas for the third column we take g ~ 0.39 as predicted by the relativistic quark model pTtj. 



Table 9: Theoretical predictions for D* and B* widths and branching ratios. The radiative decay widths 
are computed by the parameters of the preceding table. 



Decay rate/ BR 


X-loop 


VMD 


RQM 


V{D*+) 




39.5 KeV 


40.0 KeV 


46.2 KeV 


BR{D*+ 


^ D+n") 


31.5% 


31.1% 


31.3% 


BR{D*+ 


^ D'<7r+) 


68.1% 


67.3% 


67.7 % 


BR{D*+ 




0.4% 


1.6% 


1.0% 








28.3 KeV 


37.1 KeV 


41.6 KeV 


BR{D*'^ 




63.6% 


51.5% 


50.0 % 


BR{D*^ 




36.4% (input) 


48.5% 


50.0 % 




not) = 


r{Dt ^ Z),7) 


0.06 KeV 


0.35 KeV 


0.38 KeV 




T{B*+) 


= T{B*+ B+7) 


0.14 KeV 


0.37 KeV 


0.24 KeV 


T{B*") = 


= T{B*" B"j) 


0.09 KeV 


0.12 KeV 


0.092 KeV 


r(B*) - 


r{B* ^ Bsj) 


0.03 KeV 


0.09 KeV 


0.08 KeV 



We see that the chiral loop calculation, which uses the data on -> D°7 to fix the light mass scale, 
reproduces quite well the D*+ decay branching ratios; also the quark model and the VMD predictions 
(that are parameter free) are in reasonable agreement with the data. 
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Let us compare these results with other approaches. The resuh of the calculation in [Q, based on 
the ideas of HQET, is as follows: T{D*° D°-f) = 8.8 ± 17.1 KeV and T{D*+ D+-f) = 8.3 ± 8.1 
KeV, if one uses nic = 1.7 GeV. For the B case, with mf, = 5.0 GeV, T{B*") = 0.13 ± 0.20 KeV, and 
T{B*^) = 0.66 ± 0.93 KeV are found. Clearly the results of this calculation are dominated by large 
experimental uncertainties. 



The result of a QCD sum rule calculation [140|, which updates previous analyses |13£], is as follows: 
T{D*° D°-f) = 2.43 ± 0.21 KeV, T{D*+ D+j) = 0.22 ± 0.06 KeV, and T{D; Dsj) = 0.25 ± 0.08 
KeV. D* decays have been also studied in the framework of the bag model |138| with the following results 
(for the value A = 1 of the relevant parameter in that paper): T{D*+) ~ 80 KeV and r{D*°) ~ 60 KeV, 
a factor 1.5 — 2 larger than the results contained in table ^; this model also predicts r(_D*+ D^l) — 1 
KeV and the correct ratio {D^ ^) / {D^ tt^) . Correspondingly, the computed B* radiative width is also 
larger than the entries in table |^. 

Let us finally observe that similar calculations can be performed for the radiative decays of positive 
parity charmed meson resonances ]136[ Qq, 142 1. For neutral resonances the computed branching ratios 



are of the of the order 10 — 10 , for the charged parity resonances the computed branching ratios are 
much smaller, due to an almost complete cancellation between the two contributions in the e.m. current 
|136|. 



6.2 Weak radiative decay: B tvy 



Another interesting process where the formalism of the effective chiral lagrangian can be applied is the 
radiative leptonic decay channel: 



B- 



(288) 



It has been suggested [143, 144, 145 1 that this decay channel can be use d to e xtract the B* decay constant 

F /\/Mb: the analysis of (288) can represent an alternative 



fs* ; since in the — > oo limit Jb* = fs 
way to measure F as compared to the purely leptonic decay channel 



B ^ . 

The branching ratio for the purely leptonic channel is given by: 



BR{B- 



2.6 



Vub 



0.003 



200Mey 



10" 



(289) 



(290) 



where one uses t^- = 1.55 ps. This result is two order of magnitudes smaller than the present experi- 
mental upper bound put by CLEO |l46t : BR{B~ /x~P^) < 2.1 • 10"^. For the channel B^ — > e^Pg 
one expects a much smaller BR (of the order 10~^^ — 10""'^^), because of the helicity suppression; on the 
other hand in the channel B tv^., the helicity suppression is absent and the expected BR is of the 
order 7 x 10"^, but the t identification represents a serious experimental problem. 

Because of the small value that is expected for the leptonic channel (p9G), the radiative decay (|288|) 



may be a serious competitor. Various estimates [145, 147 1 of its branching ratio indicate that the radiative 
decay rate is larger than the purely leptonic one by almost an order of magnitude, mainly because the 
radiative decay, differently from the leptonic channel, is not helicity-suppressed due to the photon in the 
final state. 

Let us now describe how one could extract the value F from (28S). First of all one should dis- 
tinguish between the two classes of diagrams describing the radiative process. The first class contains 
bremsstrahlung diagrams where the photon is emitted from the B~ or from the charged lepton leg. This 
contribution vanishes in the limit m„. and is negligible also for finite lepton mass. 



The relevant diagrams for this process are of the type depicted in fig. [Tc 1 14L ] . Other possible diagrams 
are chiral loop contributions, where, instead of the single particle intermediate state, one has a chiral loop 
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(a) 



(b) 



Figure 10: Diagrams dominating the B — > £ Vij decay mode in the limit mi 
meson state, B** is the 1+ axial vector meson state. 



0. B* is the vector 



with the B and a pseudoscalar particle (these contributions are discussed in |144|). One considers only 
the resonant pole diagrams, for which no problem of double counting arises. Let us assume, following 



[144 1 and [145|, that in the pole diagrams the intermediate state is a = 1 {B*) or a positive parity 
B** meson. The amplitude with intermediate P{— B*,B**) state is written as follows: 



M 



SD 



V2 



VubA{B ^ P7) 



{p — k)^ — rnj 



<0|u7^(l-75)6|P> 



(291) 



where p and k are the B and photon momenta respectively, = £{pi)^p,{l — 75)2^(^1^) is the le pton 
current, A{B P7) is the amplitude of the process B P7, and P indicates the pole. When in ( |29l| ) 
one takes P — B*, the matrix element becomes proportional to fs*', therefore, if the contribution of the 
higher mass resonances is negligible (and we shall see that this is indeed the case), and for light leptons 
in the final state, the radiative flavour changing B decay can be used to measure the decay constant 
/b*, provided the amplitude A{B* — > Bj) i s kn own. A direct measurement of the B* width would be 
extremely difficult (as mentioned in section 6.1, it is less than 1 KeV in all the models); however this 
amplitude can be indirectly obtained by using the heavy flavour symmetry already employed in section 
|6.l| to relate the radiative D* and B* decays. 

If the partial width T{D*^ D'^l) is measured (only the branching ratio is available so far, see table 
|l|, then one could extract, from the amplitude 



eg 



A{D:{v,7j)^Da{v'h{q,e))^ I 
the mass constant to be used in the formula giving A{B* 



(292) 



Bl): 



A{B:{v,r,)^Ba{v'Mq,e)) 



Sb ^ eg 

Afc Ag 



(293) 



The expression for the amphtude M a^j giving the contribution of the B* pole to the decay B' 



fi i'p7 is therefore as follows: 



M 



ClfE 



{vk + AY 
where A = Mb' — M^, and Ci is given by: 



efc 2 
Ah 3Ai 



(294) 



(295) 
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From (294) one can compute the contribution of the B* pole to BR{B —>■ fj, V/i'y) as a function of the 
parameter F. Before considering some numerical predictions, let us study the effect of the B** pole, i.e. 
the positive parity heavy meson having si — 1/2 (only the axial vector states 1+ can contribute as poles 
to the decay; moreover, the state Bi, having si — 3/2, has vanishing coupling to the weak current in the 
limit mt, oo (see the discussion after (|l85| ); therefore only the state Bi (with — 1/2) gives a 
contribution in the same limit. 

The contribution of this state is: 



M 



(Bi) 
SD 



C2fE 



{vk + A') 



(e • vk^ - V ■ kefj_)U' 



where A' = Mb^ - Mb ^ 500 MeV, 



< 0|67p75(?|Bi(p,?7) >= fB^MB^ilf, 



and 



C2 



^/2 2Mb 



2ebri/2(l) ea_ 
A' 



Mf. 



(296) 



(297) 



(298) 



The function rx/2(i' ■ v') is the universal form factor, analogous to the Isgur-Wisc function, that appears 
in the matrix element of the heavy quark current between a positive — 1/2 and a negative parity 
heavy meson state. ti/2(w • v') has been introduced in |148| and computed in |149| by QCD sum rules, 
with the result ti/2(1) ~ 0.24. 



Using this input and the ratio {f b^ / f B*){^a^'a) estimated in |145|, one obtains 



0.1, 



(299) 



which confirms the previous hypothesis, i.e. that the Bi pole represents only a small contribution to the 
final result. 

The results are sensitive to the experimental input T{D*^ D^l)- For reasonable values of this 
quantity, which has not been measured yet, one finds for BR{B^ /i^i/^7) a result in the range 
10~^ — 10~^, i.e. a radiative branching ratio larger than the leptonic BR{B~ IJ'~v)\ the enhancement 
is obviously still higher if the electron leptonic decay channel is considered, which contributes by a factor 
of 2. Therefore, in principle, the decay channel B^ j.i^Df^j can be used as a way to measure the 
leptonic B* decay constant. 



6.3 Weak radiative decays: B — > V'j 

This section is devoted to the analysis of the exclusive flavour changing radiative B decays with a vector 
meson in the final state. These channels, at short distances, are described by processes 



(300) 



or 



d'-f , 



(301) 



which are dominated by a peng uin d iagram with the top as intermediate quark; these decays have been 
intensi vely studied in the past |l50| and experimental data have been collected both for the inclusive 
decay fsil: 



BR{B X,"f) = (2.32 ± 0.57 ± 0.35) x 10" 



(302) 
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and for the exclusive decay process 



for which the foUowing result has been obtained |152|: 



BR{B K*-i) = (4.3:;:}-J ± O.e) X 10 



(303) 



(304) 



The process (303) has been studied within the framework of the heavy quark effective theory |153 
155 , |5^ , 15£] (for reviews see |157|), by QCD sum rules |15S, 1591 , 



160 



161 



154 



16^ or Lattice QCD |163|] . 

In this section we shall review the application of the heavy meson effective chiral lagrangian |164| to the 
decay process (303) as well as to the other exclusive decay channels 



Bs 
B 
B 



PI 



(305) 
(306) 
(307) 



Before doing this, let us consider however the simpler approach already considered in section 5.3. 1| , 
consisting in the application of the light and heavy flavour symmetries [155]. 

Let us begin with the 6 — > 57 transition. All the exclusive processes arising from this elementary 
decay are dominated at the quark level by the short distance 6 — > 57 hamiltonian 1 15C | given by 



CmbS(T''''(l + 75)6F^. 



h.c. 



(308) 



(neglecting terms of order ms/mf,). i^^^ is the electromagnetic tensor, ct^^ ~ \\~1ht1v\ and C is given by 



G = 



V2 167r2 



ml 



(309) 



where F2 is a factor including perturbative QCD corrections, and slightly dependent on the top quark 
mass. For mt — 175 GeV, it has the value: F2 — 0.63. Together with the short distance hamiltonian, one 
should take into account also long distance effects, such as transitions mediated by four quark operators; 
they will be discussed below. 

The short distance hadronic matrix element relevant to the transition B — > K*^ (B — B^ or S") can 
be expressed as follows: 



(i^*(p',e)|sa'^''(l + 75)&|B(p)) = ^[A{q^)l 



(310) 



where we have used the property ^e'' crxalb = 



Now, as first noted in |153| , as a consequence of the equations of motion of the heavy quark 

in the h rest frame one has 

7°6 - 6 , 
which means that 

9^0.(1 + 75)Q = -»g"7.(i- 75)0. 



(311) 
(312) 
(313) 
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Therefore the form factors (31C) can be related to those describing the weak semileptonic transition 
B —f K* (or B ^ p, using SU(3) symmetry). In computing the width for the decay B K*"f only 
the combination of form factors A + B is relevant, whereas H{q^) contributes to amplitudes with virtual 



photons. Using (313), the form factors A, B and H are related to the form factors V, Ai and A2 as 
follows: 



A{q') 
H{q') 



q 



Ml 



V{q' 



Mb 
2i 
Mb 



Mb 



Mb + rriK* 



Mb + niK* 
Ms 



^1(9') 



rriK' 



V{q^ 



Mf 



+ 



1 



q 



Mi 



'if* 



rriK' 2q'^ Mb + niK* 



A2{q^ 



(314) 
(315) 
(316) 



These relations are strictly valid for q^ 



However, following [153, 165, p8| (see however 156 ) 



one can assume their validity down to the value q^ = which is the kincmatical point relevant for decays 
with a real photon in the final state. 

In order to compute A{0) + B{0), one needs the values of the form factors for the transition B K*; 
they have been computed hy D K* semileptonic decays using heavy flavour symmetry (see table |^ of 



1.3.1 



neglecting SU{3) flavour breaking). Using the soft pole column result, we have V 



B->K' 



(0) 



section 

0.64 and A^^^ (0) = 0.21. We observe, however, that, on the basis of the scaling relations given in the 
previous section, eq. (241), since in this case we assume the same pole behaviour for V{q^) and Ai(g^), 
the term in Ai in eq.(314) is subleading (in l/mg) for any value of > and should be neglected in 
comparison to the term proportional to V{q^), if one works consistently at a given order in the l/mq 
expansion. In this way one gets the result: 



|A(0) +B(0)| = 0.53 



(317) 



Following the discussion in section 5.3.1, we may estimate for this result a theoretical uncertainty of at 
least ±30%. We consider now the radiative width, which is given by 



r{B 

One gets 



Ml - ml, 
2Mb 



^i^^|A(0) + B(0)|2 



BR{B ^ K*-i)= 2.4 X (|l^t3|/0.039) 



X 10" 



(318) 



(319) 



In the previous formula we have used \Vt 



tb\ 



1, mb — 4.7 GeV and tb+ — tb" — tb^ — 1.55 ps. The 



result based on (317) agrees within the errors with the experimental findi ng ( 304). It should be noted 
that eq.(31£) does not include long distance effects due the cc quark loop |16£] that raise the branching 
ratio by ~ 20%, thus improving the agreement; we shall discuss them in more detail below. A similar 
analysis, with obvious changes, applies to the decay Bs — > 07 and one obtains 



BR{Bs ^ cf>j) « BR{B ^ K*7) , 



(320) 



due to approximate SU{3) light flavour symmetry. 

Now we discuss how the 6 — *■ S7 exclusive decays can be described by the effective chiral lagrangian 
approach. 

At the lowest order in the derivatives of the pseudoscalar field, the weak tensor current between light 
pseudoscalar and negative parity heavy mesons is as follows: 



(321) 
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and it has the same transformation proper ties o f the quark current q'^a^'^{l + ^5)Q. Together with (321) 
we also consider the weak effective current (17S) corresponding to the quark V — A current g°7^(l — 75)Q: 



i: = 4(7,.(i-75)ff6eL) 



(322) 



We put the same coefficient iF/2 in both (321) and (322) because, due to (313), the relation 

Lo^ = -~iU (323) 
must be satisfied. 

We also introduce the weak tensor current containing the light vector meson p°' and reproducing the 
bilinear a^"" {I + -fz)Q 



(324) 



is related to the vector current L^^, eq. (|184|), introduced to describe at the meson level the quark 



current operator q°'^^(\ 



75)(5 taken between light vector particles and heavy mesons: 

- V^)6c4> 



(325) 



(we keep only the leading term in l/rng). We notice that in order to construct the tensor current we 
have imposed (323). 

To compute B K*"f we consider a pole diagram having as intermediate state between the current 
and the BK* system either a 1"*" or a 1~ heavy meson; moreover we add a direct term. The effective 
lagrangian and the effective tensor currents can reliably describe the process only for large values of , 
i.e. q^ « Imax — {Mb — rriK*)^- This is a general feature of the chiral lagrangian approach. Again, in 
order to extend our results to small values of q^ , we shall assume a polar dependence in q^ (with pole 
mass suggested by dispersion relations). By following a procedure similar to section 5.3.2, we obtain the 



results of table |10| that are valid for any q and in the limit mg — > oo. We notice that in writing the 
various contributions in table |l^ we have left the dependence of p- p' on q^^p-p' — {Mg + mj^, — q^)/2, 
in the term arising from the 1~ pole and we have assumed that the direct term has a pole dependence 
with mass given by the 1+ pole. These choices can be justified as follows. The results in table |l^ satisfy, 
for q"^ « QmaxJ the relations (314 - 1316| ) between form factors of vector and tensor currents. Eqs. (314) 

, the result of |15 



and (315) coincide with the relations found in ref. |153|; as for 
Mb I Mb + rriK- 2q^ 



H{q' 



Ml 



Me 



A2{q^) 



+ 2mK'AQ{q^)-{MB+mK')Ai{q^)^'^ 



(326) 



differs from ( |315|) for terms that are subleading in the limit mq oo and can be neglected. 

Following ]153| ] and |165| and our previous discussion, we assume that the results ( ^14| - |316D hold 
also for small values of g'^, which justifies the above mentioned choices in table ^ 

Before computing the entries of table ^ let us observe that the results in the columns " Direct" and 
1^ are subleading as compared to those in the column 1^. In other terms 



{A + B)d. 



{A + B),- ^\m,) ^^^^^ 

for any value of (positive) q^ . Therefore, consistently with the neglect of the 0{l/mi,) contributions, we 
do not keep them, which means that only the term arising from the exchange of the 1^ particle is taken 
into account. In this way one obtains: 



A{q') + B{q') 



2 FXgy 
Mb Ml, 



(g^ + Ms-ml-O 



(328) 



62 



Table 10: Terms contributing to the various form factors of the transition B K*^. mp is the pole 
mass {mp = 5.9 GeV for the direct and 1+ term; and 5.43 GeV for the 1^ contribution), p ■ p' = 
{Ml+m\,^q^)/2. 

Form Factor Direct 1 1+ 



A{q^) 



iV2g 



vOL\ q„ 



m. 



i2V2FXgv{p-p') -i^/2M^F+gv{C,-'i^J^mK') 



B{q') 



H{q') 



i2V2FXgi 



'i2^/2M^F+gvfi 
{ml - q^)MB 



Table 11: Theoretical values of the B K*j coupling \A{Q) + B{0)\ in different approaches: soft pole, 
chiral lagrangian, QCD sum rules calculations based on the evaluation of three point function and light 
cone sum rules respectively. 

Model \A{0) + B{0)\ 



soft pole (sec. 5.3.1 



X-lagrangian eg. ([328 ) 



QCD sum rules 
QCD sum rules 



159 



161 



0.53 

0.41 
0.70 ±0.10 
0.64 ±0.10 



which gives, for F = 0.30 GeV^/^ gv = 5.8 (see eq. ^) and A = -0.41 GeV~^ (eq.(|83)), the result: 
1^(0) +B(0)| = 0.41, (329) 



1.4 X {Vts/0M9Y 



10 ^. A similar analysis, with obvious changes, 
In this case one obtains \A{0) + 5(0) | = 0.42 and BR{Bs ^7) = 
we compare the analyses based on the scaling approach and 



In table 



11 



and, therefore, BR{B 
applies to the decay Bg - 
1.6 X (14^/0.039)^ X 10" 

on effective chi ral l agra ngian to the resu lts of QCD sum rule calculations |15£, 161]; other QCD sum 
rules analyses |16C, |162| agree with [159| and |l6l| . We note that the results based on the use of the 
heavy flavour symmetry (first and second row in the table) are generally smaller than the QCD sum rules 
outcome. 

Also lattice QCD |163| has been used to compute the transition B K*^; however in this approach 
the couplings for this decay are computed near the zero recoil point and for a heavy quark mass smaller 
than its physical value. Therefore a double extrapolation is needed to compute them and it is hard to 
compare these outcomes, that should be considered as still preliminary, with QCD sum rules or chiral 
lagrangian approaches. As for the comparison with the experimental data, as we have already mentioned, 
one should take into account also the so called long distance (LD) effects, that we now discuss. 

Let us begin with the decays ( [306D and (307), where these effects are larger. The decays B ^ and 
B — > ujj take contributions both from the short distance and the long distance mechanisms. The former 
is generated by a hamiltonian similar to eq. ( [308D , with obvious modifications {s ^ d and Vts — > Vtd)- For 
it an analysis similar to the one employed in i? — > K*^ applies, but it is obvious that this contribution 
is Cabibbo suppressed as compared to i? — > K*j, which can explain why the fr/ final state is more 



63 



difScult to measure (and indeed it has not been observed yet). Because of the smahness of the short 
distance contribution, LD effects are more important in these decays than in _B — > K*^ or Bg (j)^ . 
For _B+ — > p+7 decay, their contribution has been estimated by QCD sum rules |167, |168| . The ratio of 
the long distance to the short distance amplitudes, as expressed by 



Ald{B+ ^ p+7) 
AsDiB+ ^p+-f) 



R 



VubVu. 



VtdVtb 



(330) 



is estimated to be |f67| : R = — 0.30± 0.07, i.e. a significant contribution. LD effects mainly contribute 
to the weak annihilation diagrams and are therefore relevant for B^ P^li but less important for 
B — > p,w7. For B — > K*^ one does not expect significant contributions from the weak annihilation 
diagrams because of the CKM suppression (CKM non suppressed te rms can contribute by non factorizable 
diagrams, whose role, however, has been found to be very small |l69| ). Other LD contributions come 
from the four quark operator O2 ~ CL^^CLSL^^bL (or O2 ~ CLjfiCLd.Lj'^bL for B pj). It contributes 
via a charm quark loop (the up quark loop gives a negligible contribution), with the photon emitted by 
the charm quark line, and adds about 20% to the B K*j width [I66|; by duality this contribution 
may be seen as the result of a mechanism where K*J/'ijj are produced via O2, with photon conversion 
of J/tp [pTOl 171, 172, 173|; even though less reliable, the estimates of LD contributions based on this 



mechanism seem to agree [170| with the result of the charm loop calculation. Similar results hold also 
for the B p-f decays [166, 174, 175]. 



6.4 Weak radiative decays: B 

In this section we discuss the decays 



Ke+e-, B K*e+e 



B 
B 



K*e+e- 



They occur dominantly via a quark process b 



(331) 
(332) 

se^e^ (7* — virtual photon). In the effective 



lagrangian for b 



we have to include also the already mentioned long-distance contributions; they 



produce ip — j or ip' — j conversion, and are seen as peaks in the lepton pair invariant mass distribution. 
The effective lagrangian has been derived in [ [176| , |177| and we shall not report it here for the sake of 
simplicity. 

Let us begin with B Ke^e~. The transition B — > Kj can occur only by virtual photons and is 



described at short distance by the hamiltonian (308) and, therefore, by the hadronic matrix element 



{K{p')\sa'^-'il+^5)b\Bip))=^Siq' 



(333) 



The form factor S{q^) can be related, by using the heavy quark equation of motion and eq. (313), |152], 
to the form factors Fi{q^) and FQ{q^) for the weak transition B ^ K via the vector current. In this way 
one finds |153|: 



-F,{q') 



{F,{q')~Fo{q')) 



(334) 



Let us now consider the calculation of the transition _B — > by tensor current in the effective chiral 



lagrangian approach |164|. At the tree level the relevant hamiltonian is given by (321) together with the 
strong hamiltonian BB*tt of eq. (^). The short distance diagrams are similar to those of i? ^ K*j, 
that we have described in some detail in the preceeding section. In this case, however, at tree level, there 



is no direct coupling, and the only surviving term, in the limit 00, and for q 

contribution. 



is the pole 
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Assuming a dependence of a simple pole, with Mp — Mb* (as discussed before, this seems a 
reasonable assumption for Fx related form factors), one gets the result: 



5(0) 



1 - gV^i. 



with 



5(0) 



Fg 



UMl 



-.(Mb- + Mb - ruK) 



(335) 



(336) 



When expressed in terms of the form factors of the B 
2Fi(g2) 



Mb 



K transition by vector current, this result is: 

(337) 



which coincides with the Isgur and Wise relation (334) only at q ~ qfnax ^^'^ — *■ oo, namely in the 
range of validity of the effective hamiltonian. Exactly as in the case oi B ^ K*j transition, we find that 
some form factors (in this case Fq) are subleading when mg oo, which is expected because the 0"*" 
state, contributing to Fq, cannot couple to the antisymmetric tensor current sap^ui^ + 75)^- 



The numerical result of this analysis is as follows: using g : 
5(0) ~ -0.13 GeV'^ . 



0.38 and F = 0.30 GeV^/^, one obtains 

(338) 



This result already contains some SU(3) corrections, namely in the meson masses, but the bulk of the 
chiral corrections should come from loops containing pseudoscalar bosons (tt, K and rj). They have been 
computed in . Three classes of corrections are found. First one has correction to the pole amplitude 
(|336|) ; second there are corrections to the direct term (called point contribution in |68|); last there is 
the renormalization of the B meson wavefunction. Taking into account only the nonanalytic corrections 
that arise from the loop diagrams, and having no uncertainty related to unknown analytic higher order 
terms in the phenomenological lagrangian, one finds [6^ a correction of —51% to the dominant pole 
contribution (the correction to the direct term is much smaller). 

This analysis, as stressed in ]6^ , is not conclusive, since the analytic corrections could be significant; 
nevertheless it is interesting to observe that, with the nonanalytic correction alone, the outcome of the 
effective chiral lagrangian becomes 5(0) w —0.06 GeV~^, which is compatible with the result of a QCD 
sum rules analysis based on three point functions |178|: 5(0) = —0.05 ± 0.01 GeV^^ . 



The results we have reported, together with the long distance contributions [176, 177|, can be used 
to obtain the distribution dT{B — > Ke^ e^) / dq^ in the invariant mass squared of the lepton pair . 
This distribution is dominated by the contribution of the resonances J /ip and ; however, as discussed 
for instance in ]178{ , for far from the resonance masses, one could still obtain from the data, when 
available, useful information on the short distance dynamics. 

The same analysis can be performed for dr{B K* e~) / dQ'^ . The long distance contribution 
can be evaluated starting from experimental data on the nonleptonic decay modes B — > K^J/tp and 
B K*tp'. F or th e short distance part one needs the form factors A{q'^), B{q^) and H(q^) we have 
defined in eq. (|310| ). In the effective lagrangian approach one finds, from table |l^ 



Aiq') 
B(.q') 
Hiq') 



iV^FXgvimg 



q') 



(M|. - g2)VMB 
-i2V2FXgvMl^^ 



Ml, 



-i2V2gv 



FX 



F+n 



Ml, 



Mf. 



(339) 
(340) 
(341) 
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wher e B** is the l"*" sb resona nce. ^ (0) and B{0) have been discussed in section 6.3; as for H{0), using 
(p87|) and the resuh of section |5.3.2| (obtained with r = 0.5): /Lt = -0.10 GeV"\ together with F = 0.30 
GeV3/2 and F+ = 0.46 GeV^/^, eq. (|l83|), one gets: 



|iJ(0)| = 0.04 GeV^ 



(342) 



to be compared with the QCD sum rule resuh [|78| \H{0) \ « 0.10 GeV . Analogously to the A{0) + B{0) 
case (see table we see that the chiral lagrangian approach at tree level gives significantly smaller 
results than QCD sum rules, which might indicate either a relevant l/mg correction or, most probably, 
a relevant contribution from chiral loop, similarly to the case of >5'(0) discussed above. 

The distribution in the invariant mass of the lepton pair is largely dominated by the long distance 
contributions | 179| , 172 1, exactly as for the Ke'^e~ final state; nevertheless an accurate measurement of 
the lepton pair spectrum below cc resonances would display the effects of the short distance hamiltonian. 
This measurement would therefore complement the analysis of the B K*j decay process, providing 
further information on the fundamental parameters appearing in the short distance hamiltonian as well 
as on the validity of the effective lagrangian approach. 
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7 Symmetries for heavy quarkonium states 



Quarkonium, a heavy quark-antiquark bound state, is one of the most interesting systems for the study of 
quantum chromodynamics (QCD). The physical idea is that quarks with mass larger than the QCD scale 



^QCD would form bound states resembling positronium |18C|. Many properties of quarkonium can be 
predicted by the use of non-relativistic potential models. The overall description one obtains in this way 
for the charmonium and bottomonium spectroscopies is quite satisfactory provided corrections originated 
by leading relativistic terms are included and the possible multichannel structure of the phenomenon is 



taken into account for certain expected effects |18C] 



The heavy quark and anti-quark are bound in these models by an instantaneous potential, meaning 
that gluons have typical interaction times much shorter than the time scale associated with the motion of 
the heavy quarks. We indicate with k the relative momentum and with Vr — k/mg the relative velocity 
between the two heavy quarks of mass mg . It is interesting to examine the dependence of these quantities 
on the quark mass. 



For instance, Buchmiiller and Tye |18l[] have studied a QCD-motivated potential reproducing the 



behaviour 1/r for small r, and behaving as r at large distances (this model is similar to the model in 



|35| ). Analo gous results can be obtained using other models, such as Quigg and Rosncr |182| or Grant 
and Rosner |l83| , indicating that, by increasing the quark mass, the kinetic energy and the residual 
momentum increase, whereas the relative velocity decreases. Going further up with the mass mg, the 
heavy quarks separation becomes smaller, and eventually the Coulomb term of the potential energy, 
proportional to ag/r will dominate (a^ is the strong coupling constant evaluated at the momentum scale 
1/r). Taking r of the order of the size of the bound state 1/ (jnQVr) in the potential energy and equating 
its average value to the kinetic term we get < Vr >~ Ofs, with evaluated at the momentum scale mQVr, 
going to zero in the limit mg — > oo. We notice that in such a coulombic regime the relative velocity 
decreases logarithmically. 

Concerning the spin symmetry, the coupling of the gluon to the spin of the heavy quark is expected 
to be of the order pg/niQ, with the gluon momentum pg k k. Therefore the quantity k/mg = Vr gives 
information on the degree of spin decoupling, and strictly in the limit mg oo one has an exact spin 
symmetry. 

The heavy quark flavour symmetry, on the contrary, is badly broken in quarkonium systems. In general 
the gluon radiation exchanged between static quarks gives rise to infrared divergences. In a bound state. 



potential and kinetic energy play a delicate balance against each other [184|. The regularization of the 
infrared divergences then implies a large breaking of flavour symmetry because of the explicit appearance 
of the heavy quark mass in the kinetic energy. 

For charmonium, potential models give < Vr >~ 0.5, for bottomonium < Vr >~ 0.25; one expects 
corrections, even important, to the "leading order" velocity and spin symmetry description, especially 
for charmonium. 

7.1 Non-relativistic QCD description 

An effective approach to quarkonium is given by the non-relativistic heavy quark QCD description, 
which provides a general factorization formula for annihilation decay rates of heavy quarkonium |185| . It 
consists in exploiting the fact that in quarkonium the heavy quark moves with a small relative velocity and 
nonrelativistic quantum chromodynamics (NRQCD) is a good approximation. The lowest order dynamics 
is given by the Schrodinger equation for the heavy quarks. The resulting effective theory |186 consists in 
fact of a nonrelativistic Schrodinger field theory for the heavy quarks coupled to the relativistic theory 
for gluons and light quarks. Relativistic corrections can be included systematically into this picture at 
any given order in the heavy quark velocity v. 

In this framework the scales entering the problem are written in terms of the heavy quark velocity v 
and mass mg. As shown before, the typical velocity of the heavy quark decreases as the mass increases. 



67 



When mq is sufficiently large, the heavy quark and antiquark are non-relativistic {v « 1) and the scales 
'T^Qi fnqv (typical three-momentum of the heavy quark in the meson rest frame), and rngv^ (typical 
kinetic energy) are well separated: 

Aqcd < {rriQV^f < {rriQvf < m| (343) 

In NRQCD the effects at the scale mg are taken into account through the coupling constants of 4- 
fermion operators, while the effects of the lower momentum scales rnqv, mgw^, and Aqcd are included 
into matrix elements organized in terms of their dependence on v'^. The lagrangian is obtained from 
QCD by introducing an ultraviolet cut-off of the order mg, which excludes relativistic heavy quarks from 
the theory. It also excludes light quarks and gluons with momenta of order mg. Then the heavy-quark 
and heavy-antiquark degrees of freedom are decoupled by a Foldy-Wouthuysen transformation. The full 
NRQCD lagrangian consists of the part describing the heavy quarks and anti-quarks in terms of a non- 
relativistic Schrodinger theory with separate two-component fields for the quarks and anti-quarks Cheavy , 
plus a fully relativistic part for the light quarks and gluons Cught , plus a correction term 5C reproducing 
the relativistic effects of full QCD in terms of new local interactions: 

^NRQCD — flight + Cheavy + 6C (344) 

with 

Cheavy = M^o + ) V' + N^O + ) X , (345) 

where '0 and x are the two component fields for quarks and anti-quarks and and D are the time 
and space components of the covariant derivative. The term bC contains all possible gauge invariant 
counterterms, whose coefficients must be matched with QCD in order to avoid ultraviolet divergences 
in the calculation of long distance quantities and to reproduce the results of full QCD. In principle the 
NRQCD lagrangian consists of infinitely many terms. However they can be classified in powers of the 
heavy quark velocity v and their relative importance can be established. 

The annihilation of quarkonium can be reproduced in this framework only indirectly, through its 
effects on QQ scattering amplitudes. At long distance, these amplitudes can be described adding to 
the lagrangian four fermion operators that annihilate and create a heavy quark-antiquark pair. Due to 
the optical theorem the imaginary parts of the coefficients of the four fermion operators are related 
to the annihilation of heavy quarkonium. It should be noted that the annihilation decay rates for 
heavy quarkonium are small perturbations of the energy levels. In this approach the contributions to 
annihilation widths from the dimension-6 four fermion operators contain extra suppression factors, due to 
the coefficients of the operators. The widths are of order Q!^(mQ)u or smaller, while the splitting between 
radial excitations is of order rnqv^ . 



7.2 Heavy quarkonium effective theory 



Quarkonium, in the heavy quarkonium effective theory [187, 188, |189| is described as a bound state in 
the particle- antiparticle sector of the HQET ^, ^. In quarkonium systems the internal motion of the 
heavy quarks cannot be neglected due to the delicate balance between the potential and kinetic energy 
in the bound state. This suggests to go beyond the static limit to describe quarkonia states. One must 
therefore keep the kinetic energy operator even when working at the lowest order. The kinetic energy 
operator is spin symmetric, but it includes a factor of l/mg. Therefore heavy flavour symmetry is lost 
while spin symmetry is still present. The leading order lagrangian is ]18^: 



£0 = U^\^vDW) + /.W^^) - hi-H^n.D)hl-^ + UJ^^hi-^ , (346) 
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where as usual a sum over heavy quark velocities is understood. The heavy field is obtained from the 
field Q of QCD by removing the dominant part of the heavy quark momentum: 

Q = eM-^nlQVx) {l + i±^^ + h^+\x) (347) 

and h\^^ describes a quark with velocity v. In a similar way hw ' describes an anti-quark moving with 
velocity w. The lagrangian in (346) is the starting point of an effective lagrangian description of heavy 
quarkonium decays. The two velocities of the heavy quarks differ only by a quantity of the order of 
Aqc d /™q so that it is convenient to work in the limit in which the two heavy quark velocities become 



equal. This limit can be taken consistently starting from the effective lagrangian (346). In any case a mass 
dependence in the lowest order dynamics is unavoidable so that heavy fiavour symmetry is destroyed. 
In this picture spin symmetry, as we anticipated, still holds since the kinetic energy operator is spin 
symmetric. 

As to relativistic corrections (proportional to the relative velocity of the heavy quarks) and non- 
perturbative corrections, they turn out to have the same origin in this approach, namely they come 
from higher order terms of the l/mg expansion. Short and long distance contributions for the inclusive 
annihilation decays are separated by means of the operator product expansion. The distance scale is 
given by the Compton wavelength of the heavy quark. The annihilation rates are written in an expansion 
in (A/toq) where A is the inverse Bohr radius of the system. The coefficients of this expansion can be 
calculated perturbatively. 

This approach is similar to the one of non-relativistic QCD and in fact in that case the lowest order 
dynamics is basically the one obtained by adding the kinetic energy operator to the static HQET part. 
The two approaches are not completely equivalent though: for example if we consider two operators like 
the gluon field strength [iD, iD] and (iD)^, they have the same dimension, but they are not equivalent if 
an expansion in relative velocity v/c is considered, as in the case of NRQCD (for a detailed comparison 
of the two approaches see |l89| ). In any event if the same set of assumptions is applied in the two cases, 
both the approaches yield the same results up to the order (A/mg)^. 

7.3 Heavy-meson effective theory 

This approach consists in constructing a heavy meson multiplet field and writing a lagrangian including 



the exact and approximate symmetries of the problem |190|. The procedure is analogous to the one 
introduced in the preceeding chapters of this review treating of heavy-light mesons. Symmetry breaking 
terms can be easily added to the formalism as we shall show in the following. The velocity description 
and spin symmetry are still useful, but flavour symmetry is broken. As in the single heavy quark case 
]lO| , [Tl| , |l2| , an effective lagrangian describing the low- momentum interactions of heavy quarkonia with 
light mesons can be written down. The heavy quarkonium multiplets are described by a simple trace 
formalism ^ , which can be also applied to the description of the Be system. 

7.3.1 Heavy quarkonium states 

A heavy quark-antiquark bound state, characterized by the radial number to, the orbital angular mo- 
mentum I, the spin s, and the total angular momentum J, is denoted by: 

m^'+'^lj (348) 

Parity P and charge conjugation C, which determine selection rules for electromagnetic and hadronic 
transitions, are given by: 

P =(-!)'+! (349) 
C =(-!)'+'' (350) 
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and are exactly conserved quantum numbers for quarkonium, together with ,/. If spin dependent inter- 
actions are neglected, it is natural to describe the spin singlet m^lj and the spin triplet m^lj by means 
of a single multiplet J(m, I). For the case I — 0, when the triplet s = 1 collapses into a single state with 
total angular momentum J = 1, this is readily realized: 

,,.(i±i!)K..-„l<l^ . ,351) 

Here denotes the four velocity associated to the multiplet J; and rj are the spin 1 and spin 
components respectively; the radial quantum number has been omitted. The expressions for the general 
wave jMi-Mi aj-g given in the appendix C. In the sequel K^^'"^^ represents the spin singlet component 
^Ij, the spin triplet ^Ij in the wave Jt^^-t^'. 



From eqs. ( p49[ ) and (350) one has the following transformation properties of Hi and Ki under parity 
and charge conjugation: 



^l+l 










ttI 

fj.l...fj.l 






ttI-1 

^fii...fii-i 






-K' 



(352) 



i-1 V "-1+1,1,1-1 



c; (353) 

As one can easily verify, the previous transformations laws are reproduced by assuming that the multiplet 
jfj.i...fj.i transforms as follows: 

Jill.. .11, £^ ^0 J 

•J ^ ! 'Jfj.i...fj.i I 

^ (354) 

where C = ij'^j^ is the usual charge conjugation matrix. 
Under heavy quark spin transformation one has 

jpi-w 5'jMi-W5't (355) 

with 5, S" e SU{2) and [S,^ — [S',)^ = 0. As long as one can neglect spin depe ndent effects, one will 
require invariance of the allowed interaction terms under the transformation ( |356| ) . 
Finally under a Lorentz transformation A we have: 

^ A%...A'';,^D{A)J'''-'"D{A)-^ (357) 

where -D(A) is the usual spinor representation of A. 

7.3.2 Be meson states 

The study of Be meson decays gives important information about the QCD dynamics and the weak 
interactions; moreover the Be system allows to use theoretical insight and phenomenological information 
obtained from charmonium and bottomonium. One important difference is that the total widths of 
excited Be levels are about two orders of magnitude smaller than the total widths of charmonium and 
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bottomonium excited levels, as the excited Be system does not have strong or electromagnetic annihilation 
decay channels and can only decay weakly. For the Be there is a large probability for the decay modes with 
a heavy meson in the final states as it can be seen in a simple constituent quark picture. The Be mesons 
spectra a nd t heir decay modes ha ve b een studied using potential models (see for example [191, 192, 
193, 194, 195|), lattice calculations |19(:| and QCD sum rules [197|. QCD perturbative calculations |19S| 



and fragmentation functions |199| were used to study its production. The approximate spin symmetry 
independence of the system can be implemented in an effective meson lagrangian and the corresponding 
symmetry relations impose restrictions on the form factors of the exclusive weak semileptonic decays of 
Be [ ^00[ . In the following we shall consider only this effective theory approach, as it is related to the 
material of the previous section. The consequences of spin symmetry for hadronic matrix elements may 
be derived using a trace formalism ^, 20C | analogous to the one used for quarkonium: 



(1 + ;^) 



(358) 



where H^'^^'^ is the 4x4 matrix representing the lowest-lying pseudoscalar and vector meson cb bound 
states. Under spin symmetries on the heavy quark and antiquark, the heavy meson field transforms as 



Se H^"''^ Sl 



(359) 



The definitions are analogous to the ones given for quarkonium. We shall not examine these systems any 
longer and we shall refer the interested reader to the existing literature, since no experimental data are 
available yet, and, therefore, the analysis would be only speculative. To give only an example, it has been 
suggested that the above formalism can be applied to the study of the semileptonic decay Be Dlv, 
which could provide a way of extracting the mixing angle \Vub\ 120C 



8 Heavy quarkonium decays 

The heavy quark spin symmetry leads to general relations for the differential decay rates in hadronic 
transitions among quarkonium states that essentially reproduce the results of a QCD double multipole 



expansion |201] for gluonic emission. Further use of chiral symmetry leads to differential pion decay 



distributions valid in the soft regime |19C], [202| (see also [203|). At the lowest order in the chiral 
expansion for the emitted pseudoscalars we find a selection rule allowing only for even (odd) number of 
emitted pseudoscalars for transitions between quarkonium states of orbital angular momenta different by 
even (odd) units. Such a rule can be violated by higher chiral terms, by chiral breaking, and by terms 
breaking the heavy quark spin symmetry. Specialization to a number of hadronic transitions reproduces 
by elementary tensor construction the known results from the expansion in gluon multipoles, giving a 
simple explanation for the vanishing of certain coefficients which would otherwise be allowed in the chiral 
expansion. In certain cases, such as for instance ^P^ '^P2TnT, "^Pi '^P2TnT, or D — S transitions 
via 27r, the final angular and mass distributions are uniquely predicted from heavy quark spin and lowest 
order chiral expansion. 

The effective heavy-meson description of quarkonium does not seem to present special advantages to 
describe heavy quarkonium annihilation. 

At the heavy quark level such annihilations can be described introducing four-fermion operators. The 
optical theorem then relates heavy quarkonium annihilation rates to the imaginary part of QQ — s- QQ 
scattering amplitudes. 

Heavy quarkonium annihilations were among the first tests of perturbative QCD on the assumption 
that one could factor out the non-perturbative bound state features and use asymptotic freedom to 
calculate the quark-antiquark short distance annihilation process. Such an approach has met a general 
phenomenological success and indeed it has provided a basic support to the quarkonium picture and 
to the asymptotic freedom. Calculations for P state annihilation were however disturbed by infrared 
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divergences |204] appearing within the simplest perturbative-nonperturbative separation scheme and 
requiring a suitable prescription to obtain physical predictions. The separation between short distance and 
long distance effects in P-states has been recently reexamined 1 205 , leading to a factorization prescription 
which introduces an additional long-distance parameter setting the problem of infrared regularization. 

Very recently the theoretical and experimental analysis of quarkonium production has given interesting 
results. In the following we shall concentrate on quarkonium decays. It is worth mentioning that a new 
set of data [ ^06| has encouraged to deepen the theoretical understanding of quarkonium production from 
suitable models [207| to more sophisticated calculations in the framework of QCD [208|. 



8.1 Radiative decays 

Here we discuss radiative transitions QQ — > QQ"f, where the recoiling system QQ has a mass close to 
that of the radiating system. In such a case radiative transitions are usually studied in the framework 
of multipole radiation. Radiation can occur through electric or magnetic multipole transitions, when 
allowed by the conservation rules of spin, parity and charge conjugation. Radiative decays provide a 
simple test of the formalism. We expect that this approach should reproduce the well established results 
of QCD motivated potentials |209|. The application demonstrates the power of the formalism in the 
evaluation of radiative decay amplitudes between the S— and P— wave states, both for charmonium and 
bottomonium. If absolute predictions are to be made, the formalism requires data to fix the unknown 
parameters in the effective lagrangian. 

The analysis of radiative decays in quarkonium can also be carried out directly in terms of reduced 
matrix elements of the appropriate interaction Hamiltonian, using the usual angular momentum proce- 
dures 1 210 . The two procedures are equivalent, as in this approach spin and angular momentum are 
described directly within the multiplet field. 

We write the lagrangian for radiative decays as follows: 



C = Y^ S{m, n) < J{m) J^,{n) > VyF^"' + h.c 



(360) 



where a sum over velocities is understood, F'^'^ is the electromagnetic tensor, the indices m and n 
represent the radial quantum numbers, J(m) stands for the multiplet with radial number m and S{m, n) 
is a dimensional parameter (the inverse of a mas s), to be fixed from experimental data and which also 
depends on the heavy flavour. The lagrangian (36C) conserves parity and charge conjugation and is 
invariant under the spin transformation of eq. (356). It reproduces the electric dipole selection rules 
A£ = ±1 and As = 0. It is straightforward to obtain the corresponding radiative widths: 



'5i7) = 



• Ms, 



-P 



(2J+ 1)^2.3 Ms, 



97r 



Pj 



(361) 



(362) 



r(iPi ^ 1^07) 



^ 3^ 

37r^ Mp 



(363) 



where p is the photon momentum. Once the radial numbers n and m have been fixed, the lagrangian (360) 



describes four no spin-fiip transitions with a single parameter; this allows three independent predictions. 
For the triplet states they are reported in table 12, where we give the ratio of the width for the state 
with J = 1 to the state with J = and that for the state with J = 2 to the state with J = 0, within a 



given multiplet. The theoretical numbers refer to eqs. ( [361^3 6 3D and to the results obtained by Cho and 
Wise ilT 
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For the state hc{lP) ( ^Pi state) no data on radiative widths are available yet, but extracting the 
value of S from the data of the corresponding triplet states, we can predict the width of hc{lP) rjcj- 
This transition is an electric dipole (El) transition, which is expected to be the dominant decay mode 
of the hc{lP) state (with a branching ratio of order 80% [211|). Using for the hc{lP) mass the value 
3526.2 ± 0.2 GeV of E760 experiment fll, we obtain: 



T{h^{lP) f]c-f) = 0.45 ± 0.02 MeV 



(364) 



We note that QCD-motivated potentials using first order relativistically corrected wave functions give 
a prediction of the width T{hc{lP) rjc'y) of 0.39 MeV |213|, while QCD predictions based on the 
factorization formulas of [ p05| give 0.45 ± 0.05(±20%) MeV for the partial radiative width and 0.98 ± 
0.09(±22%) MeV for the total decay rate of he, to be compared with the experimental upper bound 
on the total width of rr(/ic(l-P)) < 1-1 MeV at 90% C.L. [|l2|. A similar prediction for the Xb{ ^Pi) 
state can be easily extracted from eq. ( [363[ ), once the mass of this state and the S parameter of the 
corresponding multiplet are known. 



Table 12: Results for Ti/Tq and r2/ro, where Tj stays for the radiative width of the process involving 
■^Pj. The number in parentheses have been obtained by Cho and Wise (see text) 



Process 




Ti/To{th.) 


ri/ro(ea;p.) 


T2/To{th.) 


r2/ro(ea;p.) 


i;{2S) ^ 


Xc{l'Pj)l 


0.82(0.85) 


0.94 ±0.12 


0.55(0.58) 


0.84±0.11 


Xcil'Pj) 


^ J/V'(15)7 


2.05(2.11) 


2.61 ± 1.24 


2.74(2.84) 


2.93 ± 1.34 


T{2S) 


Xb{l'Pj)l 


1.58(1.56) 


1.56 ± 0.41 


1.56(1.54) 


1.53 ±0.41 


T(3S') ^ 


Xb{2'Pj)l 


1.61(1.61) 


2.09 ±0.26 


1.76(1.76) 


2.11 ±0.28 


Xbil'Pj) 


^ T(15)7 


1.25 




1.45 




Xb{2'Pj) 


- T(15)7 


1.07 




1.12 





Xb(2^Pj) ^ T(25)7 1:29 - 1:52 



8.2 Hadronic transitions in heavy quarkonia 

An important class of hadronic transitions between heavy-quarkonium states is provided by the decays 
with emission of two pions, for example: 

■0' — > "0 TTTT (365) 

To describe these processes we use the chiral symmetry for the pions and the heavy-quark spin 
symmetry for the heavy states. The first one is expected to hold when the pions have small energies. We 
notice that the velocity superselection rule applies at — q^axi when the energy transfer to the pion is 
maximal. Therefore we expect these approximations to be valid in the whole energy range only if (Z^^q^, 
is small. 

Nonetheless a number of interesting properties of these transitions can be derived on the basis of 
the heavy quark symmetry alone. Therefore, before specializing the pion couplings by means of chiral 
symmetry, we discuss the implications of the heavy quark spin symmetry in hadronic transitions. 

As an example, we consider transitions of the type '^Si —> '^Si + h and ^Sq — > ^Sq + h, where h 
can be light hadrons, photons, etc. By imposing the heavy quark spin symmetry, one is lead to describe 
these processes by an interaction lagrangian: 

£ss' J' J > Hss' + h.c. , (366) 

where the dependence upon the pion field is contained in the yet unspecified operator IIss'. It is imme- 
diate to derive from Css' the averaged modulus square matrix elements for the transitions ^Si -^'^ Si + h 
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and ^5*0 — s-^ Sq + h with an arbitrary fixed number of pions in the hght final state h. We obtain: 

= iMsMs'lUss'^h]^ (367) 

where Ms and M'g are the average masses of the two 5'- wave multiplets; Iiss',h is the appropriate tensor 
for the emission of the light particles h, to be calculated from the operator IIss' . By denoting with dT 
the generic differential decay rate, we have: 

dr{^Si~^^ Si + h)=dTCSo^'^ So + h) . (368) 

This is the prototype of a series of relations which can be derived for hadronic transitions as a 
consequence of the spin independence of the interaction terms. In all the known cases they coincide with 
those calculated in the context of a QCD double multipole expansion. We notice however that we do 
not even need to specify the nature of the operator 11, which may depend on light fields different from 
the pseudoscalar mesons (e.g. the photon, or a light hadron, etc), provided that the interaction term 
we are building is invariant under parity, charge conj ugation, and the other symmetries relevant to the 
transition considered. Indeed the label h in eq. ( p67| ) stands for an arbitrary combination of light final 
state particles. In this sense, this approach provides a generalization of the results obtained in the context 
of the QCD multipole expansion. 

By assuming a spin independent interaction, we can easily extend the previous considerations to other 
transitions |202| . In general, as a consequence of the heavy quark spin symmetry, the allowed transitions 
between two multiplets I and V will be related by a set of equations, independently of the nature of the 
light final state h. 



8.2.1 Chiral invariant hadronic transitions 

A useful symmetry that can be used in processes involving light quarks is the chiral symmetry. It is 
possible to build up an effective lagrangian which allows to study transitions among quarkonium states 
with emissions of soft light pseudoscalars, considered as the Goldstone bosons of the spontaneously broken 
chiral symmetry. The chiral symmetry is explicitly broken through light quark mass terms, which allow 
for rarer processes that could be, in some circumstances, kinematically favored. 

The relations ( |368| ) among the differential decay rates are direct consequences of the assumed domi- 
nance of spin independent terms for the operators describing the corresponding transitions. More detailed 
predictions can be obtained by specifying the form of the operators II's appearing in the expressions of 



the interaction terms (366). We restrict here to hadronic transitions with emission of light pseudoscalar 
mesons. 

The light mesons are described as pseudo-Goldstone bosons, included in the matrix E = (see 



formulas (33), ( p4| ) and the discussion in the section | 2 .2 [) . Frequently occurring quantities are the functions 
of ^ and its derivatives Afj, and given in eqs. (p8D and (|36|). 

For the subsequent analyses we are interested in the transformations under parity and charge conju- 
gation: 



(369) 



E 




Et 










p^ 




E 


c 
— > 


E^ 






A^ 









(370) 
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Fin ally we recall that light vector mesons can be introduced as gauge particles as discussed in section 
2.4. Under parity and charge conjugation, one has: 



P, ^ 



P, ^ -pI ■ (371) 

By imposing the heavy quark spin symmetry, parity and charge conjugation invariance, and by as- 
suming that the pseudoscalar meson coupling are described by the lowest order (at most two derivatives) 
chiral invariant operators, we can establish the following selection rules for hadronic transitions: 

even number of emitted pseudoscalars Al ~ 0,2,4, ... 

odd number of emitted pseudoscalars ^ A/ = l,3,5,... 

(372) 

In fact the spin independent operator describing Al = 0,2,4,... transitions has charge conjugation 
C = +1 (see eq. (355)). On the other hand, the lowest order, chiral invariant terms with positive charge 
conjugation are: 

< (V^-Pp)(V,-p,) > (373) 

whose expansion contains an even number of pseudoscalar mesons. Spin independence of the interaction, 
on the other hand, requires that the Al = 1,3,5, ... transitions are described by C = — 1 operators. At 
the lowest order we can form just one chiral invariant term with C = — 1: 

<A,XVu-Pu)> (374) 

whose expansion contains an odd number (> 3)of pseudoscalar mesons. 

This selection rule is violated at higher orders of the chiral expansion or by allowing for terms which 
explicitly break the heavy quark or the chiral symmetries. 

To further characterize the hadronic transitions respecting chiral symmetry, we consider below explicit 
expressions for the most general operators H-w . For simplicity, we limit ourselves to those contributing 
to two or three pion emissions: 

Hss' = Ass' < ApA" > +Bss' < {v ■ Af > 

n^s = Dps e''''''^v,<Ap{V„-pa)> 

U'p'p, = App, < ApA" > gf"" + Bpp, < [v Af > g'''' + Cpp, < A^A" > 

n^"^ = Cos < A^A" > (375) 

The constants Aw, Bw, Cw and Dw are arbitrary parameters of dimension (mass)^^ , to be fixed from 
experiments. One can easily derive amplitudes, decay rates and distributions for the corresponding 
hadronic transitions. 



For instance, the amplitude for the decay (365) is given by: 



M ("Si ^3 S, + tttt) = ^'^^sMs' ^, (^Ass'Pi ■ P2 + Bss'V ■ piv ■ p^) (376) 

where e and e' are the polarization vectors of quarkonium states; pi, p2 are the momenta of the two pions. 
It is well known that the use of chiral symmetry arguments |214] leads to a general amplitude for the 
process in question which contains a third independent term given by: 

Css' 75 (e -Pie ■P2 + e •p2e -Pi) ■ (377) 
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By combining the soft pion technique with a QCD double multipole expansion, Yan [201| finds Css' = 0. 
It is interesting to note that, within the present formahsm, this result is an immediate consequence of 
the chiral and heavy quark spin symmetries. 

Experimentally the amplitude (37£) describes well the observed pion spectra in the transitions ip' 
ipTTTT |215| and T(2s) T(ls)7r7r |21£]. Th e spectrum for the transitio n T( 3s) T(ls)7r7r seems to 
exhibit an unusual double-peaked shape |217] and cannot be fitted using (376). We observe that in this 
case, due to the large available phase space, probably the soft-pion approximation is however not reliable. 



8.2.2 Chiral breaking hadronic transitions 

In this section we discuss possible chiral breaking but spin conserving terms |202|] , which are important 
for transitions forbidden in the SU{3) x SU{3) symmetry limit. Examples of such kind of transitions are 



The transitions 



(378) 
(379) 



require terms which violate also the spin symmetry and will be discussed in the next section. 

We first discuss the masses and mixings of the octet and singlet r]' pseudoscalar light meson states. 
The term which gives mass to the pseudoscalar octet, massless in the chiral limit, is 



Cm - Ao < m(E + St) > 
Here m is the current mass matrix: 



(380) 

















rris 




ives. 


e-q 





(381) 



The lagrangian (|380| ) gives, in addition, a mixing tt^ — ij: the physical states tt^ , fj turn out to be: 



f] = 1] — CTT 

where the mixing angle e is 



md 



(382) 



(383) 



The 77', which is a chiral singlet, mixes with tt", rj. Such a mixing can be described by the term 

Z:^^, = '-I^~X < to(S - Et) > r]' (384) 



where A is a parameter with dimension of a mass. At first order in the mixing angles the physical states 
are: 



tt" + er/ + e'f]' 



fj = rj — eir'^ + 9rj' 
77' = Tj' -Orj- e'n° 



(385) 
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where 



X{md - rriu) 



n 

A ( m 



md 



(386) 



and e as given in (383). 

We consider chiral violating, spin-conserving hadronic transitions between charmonium states at first 
order in the chiral breaking mass matrix. We therefore consider the quantities: 



> 



< rfi(S + 

< m(S- E^) > 



(387) 



The first one is even under parity, the second is odd, and both have C = +1. 

The only term spin-conserving and of leading order in the current quark masses contributing to the 
transition (378) is 



a'A < ^(s _ st) > +13 fW 



(388) 



where a and jS are coupling constants of dimensions {mass)~'^ . The direct coupling to rj' contributes 
through the mixing (385). The spin symmetry of the heavy sector gives relations among the modulus 
square matrix elements of the transitions between the two p-wave states. In particular we find that 



\M\'{'Po Pott) = |M|^(^P2 Pqt^) = 
and that all non-vanishing matrix elements can be expressed in terms of ^Pq 



(389) 



Pitt: 



IMW^'Pi 

IMl'CPi 

\M\'CP2 
\M\'CPi 



Pitt) 

P2Tt) 

>^ P2n) 
Pit:) 



\MfCPo 



Y^lMfCPo 
llMl'i'Po- 



Pitt) 
Pit:) 
>^ Pitt) 



= \M\'i'Po^' Pitt) 



(390) 



where tt stays for 7r° or rj. The relations (390) can be generalized for any spin conserving transition between 
I — 1 multiplets, leading to the same results of a QCD double multipole expansion [201|. Predictions for 
widths can be easily obtained from (p83). 



8.2.3 Spin breaking hadronic transitions 



We study here transitions which violate spin symmetry 1 202 . For heavy mesons there are only two types 



of operators that can break spin symmetry. The reason is that on the quark (antiquark) indices of the 
quarkonium wave function act projection operators (1 +^)/2 and (1 — ^)/2 which reduce the original 
4 X 4-dimensional space to a 2 x 2-dimensional one. Obviously, in the rest frame, the most general spin 
symmetry breaking term is of the form a ■ a, where a are the Pauli matrices. In an arbitrary frame one 
observes that any F-matrix sandwiched between two projectors (l-|-;^^)/2, or (1 — /)/2, can be reexpressed 
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in terms of cr^^ sandwiched between the same projectors: 



1 

1 1 
2 2 

1 +/ 1 

1+/ 1+/ 
1+/ 



1 



= 



1 



'2^""'' 2 2 



and analogous relations with (1 +]^)/2 (1 —^)/2. We use here 60123 = +1- Let us define 



l±ji 1±; 
— :^ — cr,,,, — 



(391) 
(392) 
(393) 
(394) 
(395) 

(396) 



In the rest frame, a^jf^ reduce to Pauli matrices. From the previous identities it follows that the most 
general spin symmetry breaking terms in the quarkonium space are of the form Gi'^crj^\ or Gj'^cr^j/^ 
with G'"' two arbitrary antisymmetric tensors. One expects that any insertion of the operator cr^t^ gives 
a suppression factor l/mg. 

A relevant example of spin breaking is the splittings of the levels in a multiplet; one can easily write 
down the spin-spin, spin-orbit and tensor terms | |202{ . 

We apply as an example the formalism to the transitions ip' — > J/i/iTr" and ip' — > J/iprj. Of particular 
interest is the ratio 



R 



which provides for a measure of the light-quark mass ratio 
rrid - mu 



rrid 



(397) 



(398) 



Using partial conservation of axial- vector current, loffe and Shifman [218| give the prediction 



R 



27 
16 



Pn 



(399) 



The calculation of R is straightforward with the heavy quark formalism. Eq. ( |399| ) will be recovered 
when neglecting the mixings tt" — 7y and rj — rj' (or a possible direct coupling of 77'). 
The most general spin breaking lagrangian for the processes ip' J/tpn^, rj is 



..pX [< J'fT^'J >-< J<J'"'J' >] VP 



iA 

— < m{Y. - S^) > +Bri' 



h.c. 



(400) 



The couplings A and B have dimension {mass)~^; the B term contributes to the ratio ( ^97| ) via the 
mixing tt'^ — 77' and rj — rj' , in the same way as the (3 coupling in ( |38§| ). There are no terms with the 
insertion of two ct; the two P and C conserving candidates 



e^^pA [< J'crP^Ja/ > + < Ja^'^J'a/ >] yPd^ < m{Y. - St) >• 
e^i^pA [< J'aP^JaP^ > + < Ja^"'J'aP^ >] < m(E - St) > 



(401) 
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are both v anis hing. Using the lagrangian (400) and taking into account the mixings (385) we can calculate 
the ratio ( |397| ) 



R = 



27 
16 



pi 


3 


rud - m„ 


2 


.Pv. 




jUs - 1/2(to„ + m-d) 





2B XU 

3A m^, - 7 



A' 



(402) 



If we neglect the mixings tt*^ — rj' and i] — rj' (A = 0) or the direct coupling of 77' (B—Q) (402) reduces to 



Eq. (402) can receive corrections from electromagnetic contributions to the transition ip' J/ipir^. It 
has been shown that such corrections are suppressed [219, 22C|. A second type of corrections is associated 



with higher order terms in the light-quark mass expansion (the lagrangian (400) is the first order of such 



an expansion); a discussion can be found in ref. [221|. 

9 Appendix A 

We list here the Feynman rules for the vertices appearing in the heavy meson chiral lagrangian and used 
in the text. Dashed lines refer to light mesons, solid lines to heavy mesons of fixed masses (Mp or M_|_) 
and . The heavy meson propagators, for a state with velocity v and residual momentum k, are 



2(i;-fc + |A) 



j^ = o- 



and 



2(^^-fc-iA) 



r = r 



(403) 



(404) 



where A = Mp* — Mp. For the 0^ and 1+ states one has similar formulas with the appropriate mass 
difference A. 







2Mf 



■9 (e(l")-9) 
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h{vq) 



Y . 

I 







Y 



^/2 



80 



I 
I 



p 







+ 











iV2gvMp Piv ■ e{p)) 



10 Appendix B 

In this appendix we list some integrals that are encountered in computing loop corrections in the effective 
chiral theory for heavy mesons [||, |l|, |lO§. In the sequel we put A = 2/e - 7 + /n(47r) + 1. 



1 







(27r)4-e q-v - A 
d^-'q 1 1 



h{m) 



(27r)4-^ g2 _ „j2 IQ^2 

d'^-'q 1 

(27r)4-^ (g2 -m?){q-v- A) ~ IQtt^ A 



1 1 



him, A) 



(405) 
(406) 
(407) 



J^''(m,A) 



where 



d-^-'q 



(27r)4-^ {q^ - m'^){q ■ V - A) 
A[Ji(m, A)g^'' + J2(m, A)^'^^^ 



/i(m) = 
him, A) = 



167r2 



TO In — 5 — TO'^A 



-2A2ln^ - 4A2f(^) + 2A2(1 + A) 



(408) 



(409) 
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and 

F{x) = Vl tanh-^yi""^ |a;| < 1 

= -y/x"^ - 1 tan-Va;2 - 1 \x\ > 1 (410) 

Ji(m,A) = (-m^ + ^A>^ + l(A^-m2)^(J) 

- ^A2(l + A) + im2(2 + 3A) 

Mm, A) = (2m2-^A>^-i(4A2-m^)F(^) 

+ ^A2(l + A)-^m^(l+3A) . (411) 
Moreover, if we put: 

then it follows that 

Ci(A,A',m) = ^-I^[AJi(m,A)-A'Ji(m,A')] (413) 
C2(A,A',m) = ^-^[AJ2(m,A)-A'J2(m,A')] . (414) 

It can be useful to write down explicitly 

C(A, A', m) = (7i(A, A', m) + C2(A, A', m) , (415) 



which is given by 
C(A,A',m) 

with 



2m^ 



9(A - A') 



H( — ,m)-H( — ,m) 
m m 



(416) 



-xA + (9^3 - -a;) log(-^ 
+ 18x^F{-) . (417) 



if(a;,m) = -9x^(1 + A) + -xA + (Qx^ - ^a;) log(^) + 
.1 



a; 

For A = A' the previous formulas reduce to 

Ci(A,A,m)= Ji(m,A)+A^:^i^^ (418) 

C{A,A,m) = ^H'{-,m) (419) 
9 m 

where H'{x,m) = ^^^^p^- 
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11 Appendix C 



In the I 7^ case, the multiplet J for quarkonium is generaHzed to Jmi - m;^ with a decomposition 

I 



(1+/) 



1 \^ pl^ia0f,, o-Pl---Mi-lMi+l---/'i 



i=l 



(1-/) 



(420) 



In the above equation, K^^'"^^ represents the spin singlet "^Ij. Since J — I, K^^'"'^' is a completely 
symmetric, traceless tensor, satisfying the transversality condition: 







(421) 



The spin triplet ^Ij is represented by iJ^/^i-^'+i for J = I + 1, Hj"^-'" ioi J = I and Hj'^i'"''^ for 
J = Z — 1. These three tensors are completely symmetric, traceless and satisfy transversality conditions 
analogous to eq. (421). Moreover, in order to avoid orbital momenta other than I, we require that J'^i - '^' 
itself is completely symmetric, traceless and orthogonal to the velocity: 







(422) 



This allows to identify the states in (420) with the physical states. The normalisation for J^i- -^' has 
been chosen so that: 



^ Jill. ..Ill T ^_ 
^ 'J 'Jill. ..Ill 



TTi-ii...m+i ttH+i 

"l + l ^ni...ni + i 

fii...fii-i 



I fll...fll 
i fll...fll 



(423) 



where J = 7° J'''7° and < . . . > means the trace over the Dirac matrices. 



For example the expression for the P-wave multiplet J'^ that can be obtained from eq. (420) is: 



1 + :^ 



ttUQ. 

H2 la 



'V2 



=''"^^^;„7/3^i7 



V3 



(7^ - v^Ho + i^f 75 



(1-/) 



(424) 
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